Piezoelectric Ultrasonic Micromotors 


Anita M. Flynn 

MIT Artificial Intelligence Laboratory 
December, 1997 




Abstract 


This thesis describes development of a new actuator technology for integrated machines 
of the future: piezoelectric ultrasonic micromotors. Ultrasonic motors offer the advantages 
of low speed, high torque operation without the need for gears. They can be made compact 
and lightweight and provide a holding torque in the absence of applied power due to the 
traveling wave frictional coupling mechanism between the rotor and the stator. 

Whereas ultrasonic motors would typically be made from a bulk ferroelectric ceramic 
such as lead zirconate titanate, or PZT, this thesis describes the implementation of a new 
idea - that of using PZT in a sol-gel form deposited directly onto silicon to create high-torque 
motors compatible with silicon integration technologies. Due to large dielectric constants 
and increased breakdown strengths of thin-hlm PZT, ultrasonic micromotors offer a factor 
of 1000 improvement in energy density over electrostatic micromotors. In a joint project 
with the Penn State Materials Research Laboratory and MIT Lincoln Laboratory, 2 mm and 
5 mm diameter stator structures were fabricated on 1 fj,m thick silicon nitride membranes. 
Small glass lenses placed down on top spun at 100-300rpm with 4 V excitation at 90 kHz. 

While generation of appropriate traveling bending waves in the stator is fairly well 
understood, less is known about how the frictional coupling and surface properties at the 
rotor-stator interface affect mechanical power output performance. This thesis proposes 
models for production of torque in rotary ultrasonic motors. Models of line contact, Hertzian 
contact and linear spring contact for Coulomb and viscous friction have been derived and 
simulations are presented which predict speed-torque curves, efficiencies and overall output 
power for various conditions of operating voltage and normal force. 

To validate these models, a set of 8 mm diameter X 3 mm tall motors, in a designed 
experiment, has been fabricated. These devices have demonstrated maximum stall torques 
of 10 -3 Nm, maximum no-load speeds of 1710 rpm and peak power outputs of 27 mW. The 
resulting peak power density is 108 more than double that of human muscle. 

This thesis further describes a laser-etching process which has been developed to fab¬ 
ricate more practical piezoelectric ultrasonic micromotors. This laser-based process pro¬ 
duces thin-hlm PZT-on-silicon stators without the need for mask alignment, wet-etching 
or fragile membranes. The large power densities and stall torques of these piezoelectric 
ultrasonic motors offer tremendous promise for integrated machines: complete intelligent, 
electro-mechanical autonomous systems mass-produced in a single fabrication process. 
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Chapter 1 


Introduction 


Ferroelectric thin films incorporated into new ultrasonic micromotors can create actuators 
with orders of magnitude more energy density than electrostatic micromotors. Because the 
coupling of vibrational energy from the stator into rotational motion of the rotor involves 
frictional contact however, these motors are not as well understood as conventional elec¬ 
tromagnetic or electrostatic motors. The aim of the research described herein has been 
two-fold: to develop a new high-torque, low-speed actuator technology useful in micro- 
robotic applications and to increase our understanding of the frictional coupling mechanism 
at the rotor-stator interface of an ultrasonic motor. 

The research effort undertaken ranges from theoretical modeling of torque production 
in ultrasonic motors to development of a completely new type of microactuator using sol- 
gel lead zirconate titanate, PZT. These devices were fabricated in a joint project between 
the MIT Mobile Robot Group at the Artificial Intelligence Laboratory (the author’s home 
laboratory), the Pennsylvania State University’s Material Research Laboratory and MIT 
Lincoln Laboratory’s Solid State Division. 

The ultimate goal has been to integrate all the subsystems of a mobile robot, the intelli¬ 
gence system, the sensors, the actuators and the power supply, onto a single silicon substrate 
and mass-produce robots much the way we mass-produce integrated circuits. This idea of 
such a gnat robot was conceived after initial proposals were made for fabricating electrostatic 
micromotors on silicon. However, the low-torque, high-speed nature of variable-capacitance 
electrostatic micromotors made them ill-suited for robot propulsion, hence the search for a 
better motor. 
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While ultrasonic motors inherently deliver high torque at low-speed, piezoelectric ultra¬ 
sonic motors were typically built by bonding a bulk ceramic ferroelectric material with very 
large piezoelectric coefficients, such as PZT, onto a non-piezoelectric material to create a 
ring-shaped bimorph structure capable of sustaining traveling flexure waves. Due to the 
bulk ceramic formulation of the PZT material, ultrasonic motors were not amenable to 
fabrication in a silicon process. Our idea was to create a sol-gel formulation of PZT, spin it 
onto a silicon wafer in much the same manner as photoresist and use silicon micromachining 
techniques to etch motor structures. This report details our initial fabrication attempts and 
relates subsequent efforts aimed at modeling and improving the devices, in both thin-hlm 
and bulk ceramic form, for increased output performance. 

Specifically, the contributions of this dissertation include: 

• Fabrication of the world’s first ferroelectric thin-hlm motor 

• Analysis of models of interface contact mechanics 

• Characterization of 8 mm bulk motors based on design-of-experiments techniques 

• Development of a laser-machining process for fabricating new thin-hlm stators which 
requires no masks and no wet-etching 

Subsequent to a brief overview of ultrasonic motors later in this chapter, Chapter 2 
details the development of microfabricated ultrasonic motors. Much of this work was origi¬ 
nally published in [Flynn et al. 92]. After initial experiments were performed with thin-hlm 
PZT ultrasonic micromotors to determine feasibility, the project was split in order to sep¬ 
arate further development of the new thin-hlm PZT materials from development of new 
devices. Research was undertaken at the Penn State Materials Research Laboratory to 
fabricate thicker hlms with larger-area electrode coverage, while a parallel effort was be¬ 
gun at the MIT Artihcial Intelligence Laboratory to model, fabricate and test bulk-ceramic 
8 mm piezoelectric ultrasonic motors to enhance our understanding of the frictional coupling 
mechanism. 

Chapters 3 through 6 describe these models of ultrasonic motors, from stator vibration 
analysis through contact mechanics. Models of line contact, Hertzian contact and linear 
spring contact for Coulomb and viscous friction have been derived and simulations are pre¬ 
sented which predict speed-torque curves, efficiencies and overall output power for various 
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Figure 1-1: There are two levels of energy conversion in a piezoelectric ultrasonic motor. Electrical 
energy is first converted to strain energy through piezoelectric elements, where time changing elec¬ 
trical fields induce vibratory motion. The high-frequency, small-amplitude vibration is then rectified 
to lower-frequency unidirectional movement through vibro-impact frictional coupling. 

conditions of operating voltage and normal force. Chapter 7 discusses fabrication and test¬ 
ing of the prototype 8 mm diameter bulk motors. It was found that these devices could 
produce maximum stall torques of 10 -3 Nm (lOgf-cm), maximum no-load speeds of 1710 
rpm and peak power outputs of 28 mW. The resulting peak power density is 108^-, for 
these motors weighing on the order of one third of a gram. Chapter 8 discusses a new 
process developed for creating tliin-film PZT-on-silicon motors on the same scale as the 
8 mm bulk motors, but free from the wafer and from other constraints of a typical silicon 
microfabrication process. Chapter 9 gives a summary and conclusions and Appendices A 
and B relate additional background on mechanics and materials parameters. 

1.1 Overview of Ultrasonic Motors 

While most electric machines convert electrical energy to mechanical energy through the 
interaction of currents and magnetic fields, a new class of actuators has arisen which uti¬ 
lizes the effect of electrically induced vibratory motion converted to unidirectional motion 
through frictional coupling. These actuators, or vibration converters, come in a variety of 
physical embodiments and can be realized via different modes of vibration. They can be 
excited through piezoelectric, electrostrictive or magnetostrictive transduction mechanisms, 
but typically, the exciting transducers are piezoelectric elements and driving frequencies arq 
in the range of 20 kHz to 150 kHz. This subset of vibromotors, piezoelectric ultrasonic mo¬ 
tors, was first invented by the Soviets [Vishnevsky et al. 77], [Ragulskis et ah 88], first 
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Figure 1-2: Traveling waves can be induced in a waveguide structure such as the annular disk 
shown here. Points on the surface of the ring move in retrograde elliptical motions. A rotor pressed 
against the stator is propelled along in the reverse direction from the propagating wave. Figure from 
[Inaba et al. 87]. 

commercialized by the Japanese [Inaba et al. 87], [Sasliida 82], [Sashida 85], [Shinsei 89], 
[Okumura and Mukolijima 87], [Hosoe 89], [Ivasuga et al. 92], [Sasliida and Ivenjo 93], 
[Tomikawa and llelia 93] and first microfabricated by the Americans [Moroney et al. 89], 
[Flynn et al. 92], 

As shown in Figure 1-1, piezoelectric ultrasonic motors have a two-stage energy con¬ 
version process. In the first stage, piezoelectric elements convert electrical energy into 
oscillatory bending motion. Depending on the geometry of the device and the form of 
the excitation, longitudinal, torsional or flexural modes of bending can be induced in the 
structure to produce either standing or traveling waves of deformation. A traveling wave 
ring-type ultrasonic motor is depicted in Figure 1-2. 

Whatever the cause of the motion, all ultrasonic motors have a common form of second- 
stage energy conversion, wherein high frequency oscillatory vibration of a stator is rectified 
into macroscopic, unidirectional rotary or linear motion of a rotor or carriage. The mech¬ 
anism for energy conversion is a frictional impact between the rotor and stator surfaces. 
While free vibration of the stator in the first stage of energy conversion is a linear phe¬ 
nomenon and the equations of motion can be formulated as an eigenvalue problem, the 
second-stage conversion of stator to rotor motion is a vibro-impact system and inherently 
displays non-linear dynamics because vibration cycles of the stator surface cease to be 
symmetric due to impact with the rotor. 
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This thesis work concentrates on the second-stage energy conversion process, the fric¬ 
tional interaction between rotor and stator. Specifically, we focus on traveling wave motors 
which have time-continuous forms of contact between stator and rotor. 

We are interested in studying piezoelectric traveling wave ultrasonic motors for a number 
of reasons. One is that they have been shown to exhibit high-torque, low-speed character¬ 
istics without the requirements for gears. The other is that, because the stator structures 
are planar, ultrasonic motor technology is symbiotic with microfabrication techniques, and 
early results [Udayakumar et al. 91] have shown that new thin-hlm forms of piezoceramics 
can yield energy densities three orders of magnitude larger than energy densities in mi- 
crofabricated electrostatic motors. High energy densities and high torques are important 
in robotics applications, especially for autonomous machines which must carry their own 
power supplies. 

While small piezoelectric motors appear promising, there are many phenomena which 
are not well understood. The form of the friction law describing rotor-stator interaction, 
impact processes, effects of surface roughness, material hardness, wear, normal pressure, 
boundary layer regimes and bearings and mounts need to be studied in order to optimize 
and design efficient, compact ultrasonic motors. 
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Chapter 2 


Ultrasonic Micromotors for 
Microrobots 


Today’s robots are large, expensive and not too clever. Robots of the future may be small 
and cheap (and perhaps still not too clever). But if we could achieve even insect level 
intelligence while scaling down sizes and costs, there may be tremendous opportunities for 
creating useful robots. From autonomous sensors to robots cheap enough to throw away 
when they have completed their task - microrobots provide a new way of thinking about 
robotics. 

Our goal of building gnat-sized robots has been driven by recent successes in developing 
intelligence architectures for mobile robots which can be compiled efficiently into parallel 
networks on silicon. Brooks’ subsumption-style architectures [Brooks 86] provide a way 
of combining distributed real time control with sensor-triggered behaviors to produce a 
variety of robots exhibiting “insect level” intelligence [Brooks 89], [Angle 89], [Connell 90], 
[Mataric 90], [Flynn, Brooks, Wells and Barrett 89], [Maes and Brooks 90]. This assemblage 
of artificial creatures includes soda can collecting robots, sonar-guided explorers, six-legged 
arthropods that learn to walk, and a “bug” that hides in the dark and moves towards noises. 

One of the most interesting aspects of the subsumption architecture has to do with 
the way it handles sensor fusion, the issue of combining information from various, possibly 
conflicting, sensors. Typically, sensor data is fused into a global data structure and robot 
actions are planned accordingly. A subsumption architecture however, instead of making 
explicit judgments about sensor validity, encapsulates a strategy that might be termed 
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sensor fission, whereby sensors are only dealt with implicitly in that they activate behaviors. 
Behaviors are just layers of control systems that all run in parallel whenever appropriate 
sensors fire. The problem of conflicting sensor data then is handed off to the problem 
of conflicting behaviors. “Fusion” consequently is performed at the output of behaviors 
(behavior fusion) rather than the output of sensors. A prioritized arbitration scheme then 
selects the dominant behavior for a given scenario. 

The ramification of this distributed approach to handling vast quantities of sensor data 
is that it takes far less computational hardware. Since there is no need to handle the 
complexities of maintaining and updating a map of the environment, the resulting control 
system becomes very lean and elegant. 

The original idea for gnat robots [Flynn 87] came about when this realization that sub¬ 
sumption architectures could compile straightforwardly to gates coincided with a proposal 
[Bart et al. 88] to fabricate an electrostatic motor on a chip (approximately 100 /urn in di¬ 
ameter). Early calculations for this silicon micromotor forecast small but useful amounts of 
power. Already, many types of sensors (i.e. imaging sensors, infrared sensors, force sensors) 
microfabricated on silicon are commercially available. If a suitable power supply could be 
obtained (solar cells are silicon and thin him batteries are beginning to appear in research 
laboratories), the pieces might begin to fit. 

The driving vision is to develop a technology where complete machines can be fabricated 
in a single process, alleviating the need for connectors and wiring harnesses and the necessity 
for acquiring components from a variety of vendors as would be found in a traditional large- 
scale robot. The microrobots would be designed in software through a “robot compiler” 
and a foundry would convert the hies to masks and then print the robots en masse. One 
critical technology presently missing is a batch-fabricatable micromotor which can couple 
useful power out to a load. 

Various types of intriguing microactuators have recently appeared. One example is the 
variable capacitance silicon electrostatic motor (which is based on the force created due to 
charge attraction as two plates move past each other) [Tai, Fan and Muller 89], [Fujita, 
Omodaka, Sakata and Hatazawa 89], [Mehregany, Bart, Tavrow, Lang and Senturia 90]. 
Figure 2-1 illustrates one such electrostatic micromotor. Another type of micromotor is 
a “wobble” motor, where one cylinder precesses inside another, again due to electrostatic 
forces [Jacobsen et al. 89], [Trimmer and Jebens 89]. Figure 2-2 illustrates a wobble motor. 
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Figure 2-1: A variable capacitance motor has a 100 pm diameter rotor which revolves around a 
bearing as oppositely placed stators are sequentially stepped with the applied drive voltages. Figure 
from [Tavrow 91]. 

In general, electrostatic motors are preferred over magnetostatic motors in the microworld 
because electrostatic forces scale favorably as dimensions shrink and because dielectric ma¬ 
terials are more easily patterned and processed than magnetic materials, especially in the 
realm of silicon processing. The three-dimensional windings required for magnetostatic 
motors would be very hard to fabricate in silicon, but the small gap sizes that allow electro¬ 
static motors to take advantage of the ability to withstand increased electric fields before 
breakdown are easily fabricated using photolithographic techniques. Electrostatic micro¬ 
motors have demonstrated successes but also uncovered limitations. Problems with these 
types of motors arise in the areas of friction, fabrication aspect ratio constraints, and low 
torque-to-speed characteristics. 

[Flynn, Brooks and Tavrow 89] provides a detailed summary of these problems and 
proposes a piezoelectric ultrasonic micromotor as an alternative approach. This structure, 
fabricated from tliin-hlm lead zirconate titanate, PZT, circumvents many of the drawbacks 
of electrostatic micromotors. 

Our idea is based on the underlying principles of commercially available ultrasonic mo¬ 
tors presently popular in Japan [Inaba et al. 87], [Akiyama 87], [Sliinsei 89], [Kumada 
90], [Sasliida and Ivenjo 93] and [llelia and Tomikawa 93], which essentially convert elec¬ 
trical power to mechanical power through a piezoelectric interaction. Mechanical power is 
then coupled to a load through a frictional phenomenon induced by a traveling wave de- 
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Figure 2-2: The wobble motor contains a rotor which is attracted to active electrodes as the drive 
voltages are sequenced around the perimeter, similar to a variable capacitance motor. Since the 
rotor is the bearing, it tends to “wobble”. Figure from [Jacobsen et. al. 89]. 

formation of the material. Piezoelectric motors display distinct advantages over traditional 
electromagnetic motors such as small size, low noise, and high torque-to-speed ratios. These 
commercially available motors however, use PZT in its bulk ceramic form, which must be 
cut and milled. 

Our contribution has been to realize that if PZT can be deposited in a thin-film form 
compatible with silicon processing, then motors can be manufactured in a batch printing 
process instead of being individually machined. 

Additionally, these motors should show significant improvements in performance over 
bulk PZT motors. That is, because the films are very thin, it is possible to apply much 
higher electric fields than in thicker bulk devices. This leads to higher energy densities. 

2.1 Advantages of Piezoelectric Motors 

Energy Density — The argument for pursuing piezoelectric ultrasonic micromotors is 
based on energy density considerations. The maximum energy density storeable in the air 
gap of an electrostatic micromotor is 

r jp 2 

2 tair^bd 

where E ^ is the maximum electric held before breakdown (approximately 10 8 ^ for 1 /mi 
gaps) and where e ai > is the permittivity of air (equal to that of free space). 

For a piezoelectric motor made from a ferroelectric material such as PZT, the energy 
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density becomes 


I-r j? 2 

2 e pzt-t^bd 

Thin film PZT can similarly withstand high electric fields (E^d — 10 8 ^-), but the dielectric 
constant is three orders of magnitude larger (e pzt — 1300eo) than air. Other types of thin 
film piezoelectric (but not ferroelectric) ultrasonic actuators have been produced [Moroney 
et al. 89] from zinc oxide, but the dielectric constant is only one order of magnitude larger 
(e zo = lOeo) than air. The greater the energy density stored in the gap, the greater the 
potential for converting to larger torques, or useful work out. 

Low Voltages — Piezoelectric motors are not required to support an air gap. Mechani¬ 
cal forces instead, are generated by applying a voltage directly across the piezoelectric film. 
Because these ferroelectric films are very thin (ours are typically 0.3p), intense electric 
fields can be established with fairly low voltages. Consequently, we drive our thin film PZT 
motors with two to three volts as opposed to the hundred or so volts needed in air-gap 
electrostatic actuators. 

Geardown — Energy density comparisons may be the primary motivators in pursuing 
PZT micromotors, but there are other advantages as well. Because this strong dielectric 
material also bends with apphed voltage, mechanical power can be coupled out in unique 
ways. Figure 1-2 illustrated an ultrasonic traveling wave motor marketed by Matsushita 
(Panasonic). Two bulk ceramic layers of PZT are placed atop one another. Each layer 
is segmented such that neighboring segments are alternately poled. That is, for a given 
polarity of applied voltage, one segment contracts while its neighbor expands. These two 
layers are placed atop one another but olfset so that they are spatially out of phase. When 
also driven temporally out of phase, the two piezoelectric layers induce a traveling wave 
of bending in the elastic body. Any point on the surface of the stator then moves in an 
ellipse and by contacting a rotor onto the stator, the rotor is pulled along through frictional 
coupling. Fast vibratory vertical motions are transformed into a slower macroscopic motion 
tangential to the surface where peak performance is attained at resonance. This geardown 
means that we can fabricate motors without the need for gearboxes. This is especially 
important when we compare to electrostatic variable capacitance micromotors which spin 
at tens of thousands of revolutions per minute [Bart, Mehregany, Tavrow, Lang and Senturia 
90]. Gearing down to a useful speed for a robot from such a motor would require a gear 
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several feet in diameter. While electrostatic wobble micromotors are also able to produce an 
inherent gear reduction, they do not incorporate the advantage of high dielectric materials 
which the piezoelectric motors possess. 

No Levitation — Friction is another major player in problems besetting micromotors. 
In a variable capacitance electrostatic micromotor, frictionless bearings are something to 
strive for, as the rotor needs to slide around the bearing. Piezoelectric traveling wave motors 
on the other hand, are based on friction - it is sliding that we need to prevent. Consequently, 
there is no need to levitate the rotor, a fact which makes a piezoelectric motor much more 
amenable to designs for transmitting power to a load. Furthermore because the rotor in an 
electrostatic variable capacitance micromotor flies above the stator, it needs to be very flat. 
Electrostatic micromotors are small, on the order of 100 /im in diameter, because of the 
difficulties in fabricating large rotors without warpage. In a piezoelectric ultrasonic motor, 
the rotor is in physical contact with the stator, so the actuator can scale to much larger 
sizes for resulting higher torques. 

Axial Coupling — The consequences of the effects of friction and stability in various 
types of micromotors force specific geometries on these actuators. Variable capacitance 
motors require a radial gap design due to stability considerations. That is, the capacitor 
plates sliding past each other are radially distributed about the bearing. Since silicon pro¬ 
cessing techniques cannot create large structures in the vertical dimension, that leaves very 
little area for energy transduction. Similarly, the physics of wobble motors constrains them 
to have cylindrical coaxial geometries. Ultrasonic motors however, due to this frictional 
coupling, can be formed into either linear or rotary motors and in addition, have the advan¬ 
tage that the rotor can sit atop the stator, creating more area over which to couple power 
out. The large planar area of ultrasonic motors is also symbiotic with planar lithographic 
techniques. 

Rotor Material — Friction coupling (as opposed to charge attraction) leads to another 
trait characteristic of piezoelectric ultrasonic motors - the rotor can be any material. That 
is, the rotor need not be a good conductor as in variable capacitance or wobble motors. 
Rotor conductivity is unimportant, in contrast to an electrostatic induction micromotor. 
Most importantly, a piezoelectric ultrasonic motor could actuate a pump, and the fluid can 
then be any solution, without regard to conductivity. 

Holding Torque — Finally, in terms of complete systems such as autonomous robots 
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or battery-operated machines, total energy consumption over the lifespan of the system 
is of critical concern. Piezoelectric ultrasonic motors, again due to friction coupling, can 
maintain holding torque even in the absence of applied power. This is a unique trait for an 
actuator that does not contain a gearbox, much less a transmission system or a brake. 

2.2 From Materials to Devices 

Bulk ceramic PZT has been widely used for decades but thin him ferroelectrics are new ar¬ 
rivals, having only recently been developed for non-volatile memory applications [Ramtron 
88], [Udayakumar, Chen, Krupanidhi and Cross 90]. One problem with these new ferro¬ 
electric memories is fatigue, as the chips actually flex when memory cells switch. But that 
is exactly the effect we seek to exploit! 

We would like films that maximize the piezoelectric effect in order to design useful 
high torque, low speed micromotors but the leap from materials to devices is a large one. 
Standing on the shoulders of previous technology is, in general, a good idea (and one which 
has been the approach in electrostatic micromotor research to this point - some even going 
so far as to label them “IC-processed micromotors” [Tai, Fan and Muller 89]). Stepping 
away from known silicon processing techniques and incorporating a new material can be a 
large undertaking, especially when the aim is to develop a new device. Consequently, the 
design for the device has to be as simple as possible in terms of materials processing to 
ensure a reasonable chance of success. 1 

2.2.1 Keeping Things Simple 

Figure 2-3 and Figure 2-4 illustrate our initial designs for the stators of linear and rotary 
motors (carriages and rotors have not been microfabricated - at the moment we simply 
place small glass lenses or other materials down on the stators). A silicon-rich nitride layer 
is deposited on a silicon wafer and is then patterned on the backside to create a membrane. 
120 stator structures are patterned per two-inch wafer. After the membranes are etched, 
piezoelectric capacitor structures are built. These structures consist first, of a bottom 
electrode formed from titanium and platinum. The PZT dielectric is then added and finally 

1 Experienced designers usually note however, that the first way you design something is always the most 
complex way [Angle 90]. 
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Figure 2-3: Linear motors utilizing thin-film PZT are illustrated I litre. By etching a membrane 
into a silicon wafer and patterning the stator on the membrane, the stator will be able to deflect 
more than if it were trying to bend the entire thickness of the wafer. Silicon-rich nitride is used for 
the membrane. The stator consists of a bottom electrode of titanium and platinum (ground), the 
PZT Him and the patterned gold top electrodes. A carriage would have to be placed down by hand. 
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Figure 2-4: A rotary motor is made in the same way except that the top electrodes are patterned 
in a circle. We typically place down a small glass lens for a rotor. 


the patterned gold top electrodes are deposited. The bottom electrode and thin film PZT 
are laid uniformly over the entire wafer, while gold top electrodes are positioned only above 
membranes. 

A close-up of the membrane cross section is shown in Figure 2-5. Note that the silicon 
wafer and the silicon nitride membrane provide only structural support for the stator. No 
electrical properties or charge attraction effects of silicon are presently used in this motor. 
Future iterations might find other manufacturing technologies more attractive, but for the 
present we use silicon for its accompanying tools and lithographic techniques. 

These stators were designed in this fashion because the materials requirements here 
are much simpler than in, for instance, a Matsushita motor (Figure 1-2), which would 
require two layers of ceramics with alternately poled segments throughout each piece. Our 
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Figure 2-5: This scanning electron micrograph shows a cross section of the nitride membrane 
structure with PZT and the gold top electrode. The titanium-platinum bottom electrode is too thin 
to see here. 

microfabricated stators require only one layer of PZT and that layer is poled uniformly 
everywhere. The tradeoff is that our motors now require a four-pliase drive to induce a 
traveling wave, whereas the Matsushita motor requires only two phases. Patterning and 
wiring is straightforward with photolithography. However, while multilayer materials with 
various geometries of poling are easily realized in macro-scale assembled motors, these steps 
would be cumbersome from a microfabrication point of view. 

A more recent design is even simpler. We start with thinned wafers which are 75 /mi 
thick (as opposed to the usual 250/mi) and omit membranes entirely. Since stress-free 
nitride is no longer required as an etch stop to create the membranes, the entire layer 
sequence simply becomes silicon, oxide, Ti-Pt, PZT and then gold (titanium is required for 
adhesion reasons and oxide is necessary to separate silicon from reacting with the PZT). 

2.2.2 Stators and Rotors 

In either case, whether utilizing membraned wafers or thinned wafers, there are a variety of 
possible geometries for patterning the top electrodes for inducing traveling waves. Figure 2- 
6 shows the simplest layout for a rotary motor. Eight electrodes are patterned radially 
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Figure 2-6: This 8-pole stator has an inner diameter of 1.2 mm and an outer diameter of 2 mm 
placed over a 2.2mm by 2.2mm square membrane. The eight pads are driven in a four phase 
sequence (sin, cos, -sin and -cos), repeated twice. The extra four pads at the north, east, south and 
west positions are undriven pads which can be used as passive piezoelectric sensors. That is, a signal 
will be generated as the wave passes through the pad. 


around a center point and driven four-pliase over two wavelengths. Eight probes would be 
needed to drive the motor in this particular example. However, other patterns on our test 
wafer use an interconnect scheme between pads to reduce the requirement to four probes. 
Note in Figure 2-6 that there are four extra pads. These can be used as sensors, since the 
piezoelectric him is reciprocal, where a bending moment can induce a voltage. 

The simplest way to observe electrical to mechanical energy conversion is to place a 
small object down on the stator as portrayed in Figure 2-7. We have found that glass lenses 
spin the best, although dust particles dance wildly, signaling the onset of resonance as drive 
frequencies are swept from 50 kHz through several hundred kilohertz. Typical rotational 
velocities of the glass lens are on the order of 100-300 rpm. One interesting point to note is 
that four phases are not necessary to induce spinning. In fact, the lenses rotate with only 
single pad excitation. This is likely due to wave reflections off the edges of the membranes 
setting up parasitic traveling waves. 

In addition to rotary stators, we have fabricated linear stators as shown in Figure 2-8. 
These structures can also be used as surface acoustic wave devices for measuring various 
properties of the ferroelectric him, such as acoustic velocity, etc. 
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Figure 2-7: Here a plano-convex 1.5mm diameter glass lens is placed convex surface down upon a 
rotary stator which has the same dimensions as Figure 2-6. Although there is no bearing, the lens 
spins at 100-300rpm when the stator is driven at 90kHz. 



Figure 2-8: Linear stators have also been fabricated. Here, the probe pads to the right are 200 pm 
square and connect to every fourth electrode for setting up four-phase traveling waves. These 
structures are patterned over similarly shaped linear membranes. 
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2.2.3 The Process 


Nitride membranes are first fabricated using standard bulk micromachining techniques into 
two-inch silicon <100> wafers. A 1 fj,m thick low-stress chemical vapor deposition silicon- 
rich (nonstoichiometric) nitride film acts both as the membrane and the mask for the tetra 
methyl ammonium hydroxide (TMAH) anisotropic silicon etch. The electroded PZT film 
(for the stators) is then formed on the membranes. The reduced stiffness of the membranes 
permits larger stator deflections than would be possible on a full thickness wafer. 

Electrode material selection is critical to fully utilize the PZT properties. We have used 
a 0.46 n m thick platinum layer for our bottom electrode which is deposited on top of a 20 nm 
titanium adhesion layer. The nitride layer together with the titanium-platinum layers act 
as a separation barrier preventing the silicon from reacting with the PZT. 

Sol-gel PZT film is deposited by a spin-on technique in a series of steps. These films have 
been characterized as reported in [Udayakumar et al. 91] and show significant improvements 
over bulk PZT, including greater breakdown strength and dielectric constant. Although 
thin-film PZT was first developed for memory devices, much of that work has focused on 
sputtering and chemical vapor deposition methods, even though it is very difficult to get 
the correct PZT makeup with these techniques. Sputtering from three separate elementary 
targets (or even from a single ceramic PZT target) to get lead, zirconium and titanium 
atoms all in their proper atomic positions in the crystal lattice is significantly harder than 
preparing a solution of the proper composition and spinning it onto a wafer as in a sol-gel 
process. These sol-gel fabricated films do in fact, exhibit the proper perovskite structure 
and show strong ferroelectric characteristics. For memory applications, piezoelectric hexing 
is not of interest, but the conformal coating properties of vapor deposited films are. Sol-gel 
films on the other hand, are planarizing, which can be a problem where uniform thicknesses 
even over undulating surfaces are required. One critical requirement for preparing quality 
sol-gel films is cleanliness, as wet spin-on techniques are more susceptible to the particle 
contamination than vacuum-based methods. 

The sol mixture is prepared from the lead precursor, lead acetate trihydrate, together 
with alkoxides of Ti and Zr in 2-methoxyethanol as the solvent. Films are spun-on in 
approximately 50 nm layers. The film is pyrolyzed under quartz lamps after each step to 
remove organics. After 6-8 layers, the film is annealed to crystallize into the perovskite 
phase, which is the type of crystalline form which brings out the strong ferroelectric and 
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piezoelectric traits. Annealing is carried out above 500 °C. These PZT films are of the 52-48 
mole ratio of zirconium-titanium which places them on the morphotropic phase boundary. 
The morphotropic phase boundary composition is that composition for which the crystallites 
have the maximum number of possible domain states because the composition lies at the 
boundary of the tetragonal phase (6 possible domain states) and the rhombohedral phase 
(8 possible domain states). This position among the possible spectrum of compositions in 
the lead zirconate - lead titanate solid solution system is the most amenable for attaining 
the distinctive ferroelectric properties. One interesting point about these thin films of PZT, 
in contrast to its bulk form, is that poling, the process of aligning domains in order to 
bring out these strong ferroelectric characteristics, need no longer be performed at elevated 
temperatures. 

Characteristic measurements, described in more detail in [Udayakumar et al. 91] are 
summarized in Table I. Similar measurements reported in the literature for bulk PZT [Jaffe, 
Cook and Jaffe 71] yield some interesting comparisons. The breakdown field strength of 
1 is significantly improved over bulk PZT which is often on the order of 40 —. Our sol- 
gel PZT films also exhibit almost twice the relative dielectric constant, 1300, of (similarly 
undoped) bulk PZT which is 730. 

Table I. Sol-gel PZT Film Characteristics 


Ebd 

i mv 

cm 

Breakdown field 

tpzt 

1300e o 

Dielectric constant 

tan 6 

0.03 

Loss tangent 

CO 

CO 

220^ 

Longitudinal piezoelectric coefficient 

^31 

—88.7 

Transverse piezoelectric coefficient 

Pr 

36^ 

cm z 

Remanent polarization 

co 

CO 

0.49 

Longitudinal coupling factor 

hi 

0.22 

Transverse coupling factor 

kp 

0.32 

Planar coupling factor 


Once the PZT film has been annealed, 0.5 fj,m thick gold electrodes are deposited and 
patterned by lift-off. A variety of eight- and twelve-pole rotary stators on three sizes of 
square membranes (0.8 mm, 2.2 mm and 5 mm per side) and various configurations of linear 
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stators are patterned. 

The structures built on thinned wafers are fabricated in an analogous manner except 
that the membrane etch is skipped and 0.5 /im oxide is used in place of nitride. The oxide 
layer together with the titanium-platinum layer act as a separation barrier preventing the 
silicon from reacting with the PZT. 

2.3 Results 

Initial experiments with these thin him PZT actuators have raised some intriguing questions. 
On the one hand, we have observed phenomena we expected such as high energy densities, 
gear down and low voltage operation. A 4 V peak-to-peak drive signal at 90 kHz competently 
spins a fairly large rotor, a glass lens 1.5mm in diameter, at 100-300rpm. On the other 
hand, the lens spins competently with only one phase excitation, and does not spin any 
better with four, something we did not expect! Furthermore, changing directions when 
applying four-phase drive does not cause the rotor to reverse, although in one instance, it 
did cause the lens to stop. Essentially, we are not inducing traveling waves in the manner 
we would like, but evidently there is enough energy density that the lenses spin anyway. 

We have observed other indications of high energy density. Not only do dust and par¬ 
ticles vibrate across the stators upon resonance, but in certain instances in which a pad’s 
impedance is very low, applying a voltage on the order of 10 V causes a static deflection 
dependent on voltage that can be seen through the microscope as a darkened area where 
the surface is deformed and reflecting light away from the eyepiece. An example with a 
unique stress pattern is shown in Figure 2-9. Note that even at 10 V, the electric fields that 
we can apply across our 0.3/im thick films contribute to energy densities will beyond those 
achievable with bulk ceramic PZT motors. 

The single-phase drive is intriguing and brings into question the effects of the boundary 
conditions on the waves imposed by the edges of the membranes. At high enough frequen¬ 
cies (several hundred kHz) standing waves become visible on both square and rectangular 
membranes as shown in Figure 2-10. Rotors continue to spin however, even though the 
motors are not working in the manner in which they were designed. The plano-convex glass 
lenses seem to spin because the contact is a point. We have observed glass lenses, convex 
side down, jiggling across the stator until brushing up against the inner radius of the gold 
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Figure 2-9: Static deflection of this partially shorted stator pad can be seen through the microscope. 
The darkened portion is deformed, deflecting light away from the eyepiece. 10 V is applied from the 
electrode at the right across the PZT, to the ground plane beneath it. 


electrodes which are approximately 0.5/mi tall, whereupon they sit and spin. We have 
also observed that plano-convex lenses flat side down do not spin, nor do annularly shaped 
objects such as jeweled bearings. 

In this first fabrication sequence, we made no attempt to microfabricate a bearing or 
etch a rotor in place. Consequently, the amount of frictional coupling is only determined 
by the weight of the rotor. In fact, it is possible there is no frictional coupling and the 
lens is merely sliding on air as the surface vibrates. Nevertheless, a mass is spinning and 
it is possible to calculate a torque by measuring the inertia of the lens and its acceleration 
when starting. Approximating our lens as a disk 1 111111 thick and estimating from visual 
observation that it reaches a final velocity of 3 Hz in one quarter of a second, the resulting 
torque is dlpNni. We can compare to variable-capacitance electrostatic micromotors by 
normalizing over voltage squared, which gives a figure of merit over a class of experiments. 
This normalized net torque for such electrostatic micromotors (typically run at 100V) is 
1.4x 10 -15 The figure of merit for our piezoelectric motors then becomes, for 5 V 

excitation, 1.6 X 10 -12 ■^ i , about three orders of magnitude larger. Wliat this number 
actually signifies for a piezoelectric motor with no bearing and no traveling wave drive is 
debatable. Mostly, it serves to underscore that the films are indeed very active, encapsulate 
high energy densities and can move fairly large objects with low voltages. 

True motor action will depend on future attempts to fabricate a bearing and to measures 
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What went wrong... pin holes and standing waves. 



Figure 2-10: Standing waves on a membrane are visible through an optical microscope when a 
single-phase drive is applied. 


torques across a spectrum of normal forces. Further experiments are needed to determine 
better structures for guiding traveling waves. Instruments need to be developed for visual¬ 
izing the waves throughout a spectrum of frequencies and for ascertaining the amplitudes 
of these dynamic deflections. Determining proper rotor-stator interface coatings for high 
friction contact would also be helpful. 

Electrostatic motors are essentially the duals of magnetostatic motors which have been 
around for years and are well understood. Piezoelectric ultrasonic motors on the other hand 
are fairly new and ferroelectric thin films are newer still. Many factors conspire to produce 
complexities and difficulties in analyzing these structures: non-linear materials, coupled 
electrical and mechanical fields, resonance drives, clamped and unclamped boundary con¬ 
ditions and friction-based interactions between rotors and stators, to name a few. 
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Chapter 3 


Ultrasonic Motors 


Due to the constraints of the silicon micromachining process, our first attempts at thin-hlm 
PZT microactuators were built on l^m thick membranes. We chose membranes in order 
to achieve a thin substrate against which the thin him could bend. Due to constraints in 
fabricating sol-gel PZT, the PZT could not be much thicker than 0.3/um or cracking would 
occur during pyrolization. 

Consequently, because the membranes acted as drumheads rather than mechanical 
waveguides, and because the films sometimes shorted due to pinholes, it was difficult to 
achieve true traveling wave motion. At this junction, we split the project into two separate 
studies: one to develop thicker, pinhole-free films at the Penn State Materials Research 
Laboratory, and one, at the MIT Artificial Intelligence Laboratory to develop models for 
predicting performance of traveling wave motors on ring-type waveguides. 

Here, we discuss a succession of models and relate early experiments on prototype motors 
made at a larger, but still small (and useful for our mini-robot applications) scale, made from 
commercially available bulk ceramic PZT. These motors are 8 mm in outer diameter with a 
5 mm inner diameter. Two generations of these motors were actually prototyped during the 
course of the project, with significant improvements incorporated into the second batch. 

Later, we will describe a new laser-based process for bringing the two projects back 
together and achieving thin-hlm PZT-on-silicon micromotors where membranes are circum¬ 
vented and the motors are cut free from the wafer. 
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3.1 Bulk 8 mm Motors 


Traveling wave piezoelectric ultrasonic motors can be made with a very simple structure. 
Figure 3-1 illustrates an 8 mm ring-type motor. The stator is composed of a steel ring 
with piezoelectric elements bonded onto the underside for exciting vibration modes in the 
ring. When traveling waves are developed, points on the stator surface move in retrograde 
elliptical motions, creating a tangential component of velocity which propels the rotor. The 
rotor can be any material, but a normal force must be provided to press the rotor against the 
stator in order to frictionally couple the vibratory motion of the stator into the rotational 
motion of the rotor. 

The most commonly used material for piezoelectric actuators, and that used in this mo¬ 
tor, is lead zirconate titanate, PZT. PZT is a ceramic material that is ferroelectric, meaning 
that it displays a hysteretic effect between polarization and electric held where the polar¬ 
ization direction can be reversed with opposite polarity drive held - but it also happens 
to have very large piezoelectric coefficients. Bulk forms of PZT, on the order of several 
hundreds of microns in thickness, are commercially available and found in bimorph and 
multi-layer stack actuators. Thin-hlm PZT, under one micron in thickness, has recently 
been developed for non-volatile memories, taking advantage of the ferroelectric switch¬ 
ing characteristics [Udayakumar 93]. Fundamentals of piezoelectric notation and material 
properties can be found in Appendix B. Here we will model ultrasonic motors utilizing the 
piezoelectric properties of both bulk and thin-hlm PZT. 


3.2 Vibration of Rings 


An ultrasonic motor is essentially a vibrating annular ring. Vibration of strings, membranes, 
plates and disks was studied by Lord Rayleigh and Kirchoff in the nineteenth century 
[Rayleigh 1894] and the equation of motion for free transverse vibration of a circular plate 
was formulated as: 


Eh 2 

3(1 - v 2 ) 


V 2 V 2 w + p 


d 2 w 


= 0 


where 

2 d 2 Id Id 2 

^ dr 2 + r dr + r 2 dO 2 

for the coordinate system as shown in Figure 3-2. 
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Figure 3-1: This 8 mm diameter piezoelectric ultrasonic motor is composed of two pieces, the stator 
and the rotor. The stator is shown in the figure at top and is a steel ring with piezoceramics bonded 
onto the bottom side. A brass rotor is shown atop the stator in the bottom photograph. 


Here, id is the transverse deflection of the plate in the ^-direction, E is Young’s modulus, 
h is the lialf-tliickness of the plate, v is Poisson’s ratio and p is the mass density of the 
material. Vogel and Skinner [1965] give a detailed analysis with numerical evaluation of 
various boundary conditions. Solutions to this free vibration problem are of the form: 


w (r,9,t) = A(r) cos n9e Jujt 


where the mode shapes are transcendental functions in the ^-direction with scaling factor 
A(r). Plugging this solution into the equation of motion leads to Bessel function forms for 
A(r): 

A(r) = C\J n {fir) + {fir) + Cfil n {fir) + G±I\. n {fir) 

where 

4 _ 'ip (1 - v 2 ) 

1 Eh 2 

and .J n and Y n are Bessel functions, and I n and K n are modified Bessel functions. 
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33 kHz 


34 kHz 


39 kHz 



Umbrella mode 


Radial mode 





Figure 3-2: Finite element simulations of the modes of vibration for an 8 mm outer diameter, 5 
mm inner diameter, composite ring of steel on PZT, show the mode shapes and natural frequencies. 
The steel layer here is 820 pm thick and the PZT layer is 195 pm thick. 
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Figure 3-2 shows finite element simulations (ANSYS finite element package, [Ostergaard 
89]) which help to visualize the free vibration of an 8 mm outer diameter, 5 mm inner 
diameter composite ring, similar to the stator in Figure 3-1. The boundary conditions for 
these structures are free-free along the inner and outer circumferences. 

We can see that simulations show that not all modes of vibration are composed of trans¬ 
verse displacements. Such transverse deformations are known as flexure or bending modes, 
and the 2-, 3-, and 4-wavelength flexure modes can be seen to have natural frequencies 
of 12 kHz, 33 kHz, and 63 kHz, respectively. Due to axial symmetry, these modes are 
degenerate, having eigenvalues of multiplicity two. By exciting both solutions as standing 
waves, but phased 90 degrees apart in time, traveling waves result due to superposition. 
For an ultrasonic motor that is top-drive, where the rotor sits atop the stator, we want 
to induce traveling waves from these flexure modes. Other modes, such as the umbrella 
mode, will not sustain traveling waves of bending. Note, however, that the 39 kHz radial 
mode could be used for an ultrasonic motor if the rotor was placed circumferentially around 
the stator. Kumada has used this technique to produce very thin clock motors [Kumada 
91]. Here, we will study top-drive flexural traveling wave motors because we are interested 
in investigating piezoelectric thin films for excitation, a technology more compatible with 
top-driven devices. 

Because the stator of our motor has free-free boundary conditions at the inner and outer 
diameters, and the radial dimension of contact with the rotor is very short, we can model 
the radial variation of vibration amplitude as constant, and unfold the annular ring into a 
infinite Bernoulli-Euler beam capable of sustaining traveling waves. For a Bernoulli-Euler 
model and small deflections, rotary inertia and shear forces are neglected and normal cross- 
sections are assumed to remain normal after bending. While a finite beam can only support 
standing waves due to reflections at each end, we can model one standing wave component 
of the traveling wave as a finite beam of the same length as an open, unfolded ring as shown 
in Figure 3-3. 

Figure 3-3(a) represents the four-wavelength flexure mode at which the motors of Fig¬ 
ure 3-1 are designed to run and Figure 3-3(b) illustrates the Bernoulli-Euler beam model 
of one of the standing wave components. Using this model and superposition, eigenfre- 
quencies, traveling-wave speeds and surface-point trajectories can be determined. To pre¬ 
dict deflection in the out-of-plane direction, we may take one half-wavelength section of 
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Figure 3-3: (a) The fourth flexure mode vibration pattern of an annular ring, (b) Bernoulli-Euler 
beam model of the ring in its fourth flexure mode. 

this standing-wave beam and model it as a simply-supported piezoceramic-metal-composite 
beam of length L , where L = the supports being located at the nodes of the standing 
wave, as illustrated in Figure 3-7 and discussed later. 

3.3 Eigenfrequencies and Wavespeeds 

With the simple beam model of Figure 3-3(b), we can calculate the natural frequencies of 
an ultrasonic motor and the corresponding wavespeeds and rotor speeds. 

For a beam of transverse displacement w, mass density per unit length p, Young’s 
modulus E , cross-sectional moment of inertia /, and cross-sectional area A, the equation of 
motion of flexural free vibration of a beam is: 

d 4 w(x,t.) pAd 2 w(x,t.) 
dE 1 + Tl W 2 = 

Solutions are of the form: 

w (x,t) = W ( x ) 

with pinned boundary conditions at the nodes of a standing wave where L = 

d 2 W d 2 W 

W (x = 0) = W (x = L) = —-7- (x = 0) = — (x = L) = 0 

a,x* a,x* 
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yielding mode shapes with a sinusoidal form: 


W (x) = C\sinkx 


Here k is the wavenumber: 



w 

Cy 


A is the wavelength and cy is the speed of the transverse bending waves. For a pinned-pinned 
beam, 

, mr 


Plugging this solution back into the equation of motion, we find the dispersion relation: 


l 4 _ , ,2 PA 




El 


and the natural frequencies: 



We can also note that: 



where c; is the longitudinal speed of sound in the material. Furthermore, we can take: 



where k is the radius of gyration. This shows that the speed of flexure waves is dependent 
on excitation frequency: 


LO 

CT = T 




3.4 Traveling Waves and Elliptic Motion 

The wavespeed we have just described is the phase velocity of traveling flexure waves along 
the neutral axis of a beam. However, for a beam of half-thickness h, the points on the 
surface move in elliptical trajectories. Coupling of this motion to the rotor produces a rotor 
speed different than the wavespeed. 


35 



d 4 w(x,t ) pAd 2 w(x,t) 

Flexing Beam - + —-^- = 0 

OX 4 £7 C/£“ 

Traveling Wave Solution w(x,t) — W 0 COs(kx — out) 



(b) 

Figure 3-4: (a) For a beam in flexure, a traveling wave is one solution to the beam equation. A 
vibration of amplitude w 0 will cause a point on the surface to displace from P to P a . (b) Traveling 
wave excitation causes point P to undergo retrograde elliptical motion. A carriage pressed against 
this infinite beam is propelled along through friction in the direction opposite to the direction of the 
traveling wave. 
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Figure 3-5: (a) A thin glue layer is used to attach an elect.roded ferroelectric to the underside of a 
steel ring, (b) The segments are then cut. and poled alternately. The poling process aligns domains 
in the piezoceramic to create a remanent polarization, (c) After poling, two groups of electrodes are 
connected as shown at the bottom with silver paint., and wires are soldered in place. 


Generation of traveling waves and equations of motion for surface points on beams and 
rings have been reported in the literature [Sa.sliida. 86], [Inaba. et a.l. 87], [Zemella. 90], 
[Hagedorn and Wa.lla.scliek 92], [Hagedorn et a.l. 92]. 

It can be shown (see Appendix A) that for a. traveling wave solution of: 


w(x, t) = w 0 cos(kx — ujt) 


a. point, P, on the surface of a. beam as shown in Figure 3-4(a), moves with horizontal 
displacement, (, arid vertical displacement after bending, C where: 


£ = w 0 cos(kx — ujt) 


( = hkw 0 sii \(kx — ujt) = 


'hthws 

A 


sii \(kx — ujt) 
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The velocity of the horizontal displacement of a point on the surface of a beam is: 


d( 




dt 


2TTUjhw 0 

A 


cos (kx — 


Lot) 


where it reaches its maximum value at the peak of the deflection contacting the rotor. 
Assuming no slip between the stator and the rotor, the rotor is propelled in the opposite 
direction of the traveling wave at a speed of: 

2lTLohw 0 

max — 7 


This maximum horizontal velocity is simply the no-load speed of the motor. Note that 
this is very different than the phase velocity, cy, of the traveling wave: 



Figure 3-4(b) shows the retrograde elliptical motion on the surface of the stator as the 
traveling wave moves to the right. The angular velocity of the rotor is 


^rotor,max 


^t,max 

r 


where r is the radius of the rotor at the circumference of contact with the stator. Comparing 
this rotor speed to the frequency of vibration, we find: 


^rotor,max 

LO 


2ir hw 0 
A r 


3.5 Mechanical Modeling of the Stator 


We can model our stator of Figure 3-1 as a composite beam in the manner illustrated in 
Figure 3-5(a). Teeth are cut in the steel ring to reduce stiffness while maintaining maximum 
height for mechanical amplification. That is, for calculation of the natural frequencies: 


^n 



we model the teeth as not entering into the stiffness term, El, but as contributing to added 
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mass per unit length: 


(M) tot (pA)beam (P^)teet/i 

where n = 1 in our model of each half-wavelength piezo-segment acting as a pinned-pinned 
beam inducing a bending moment in the structure. 

The beam is considered as the composite structure of PZT and metal, a momomorph, 
where the height of the beam is the sum of the PZT thickness and the base portion of the 
metal thickness, h p + the pitch of each tooth is p t , the length of the beam is L = 
and the depth of the beam is b. The cross-sectional moment of inertia, /, is measured with 
respect to the neutral axis. The location of the neutral axis is the distance from the bottom 
of the piezoceramic layer, g , calculated as a fraction of PZT thickness, and weighted by the 
elasticities (stiffnesses) of each material: 



Here, c m is the elastic modulus of the metal layer and c p is the elastic modulus of the 
PZT layer (equal to cf t in the piezoelectric notation outlined in Appendix B). 

The piezoelectric elements are attached to the steel ring in the manner shown in Figure 3- 
5(b) and Figure 3-5(c). The piezoceramic is purchased with a thin electrode layer of nickel 
coating on each side. After bonding, the bottommost electrode is cut into sectors. Once 
the backside electrode is segmented, neighboring segments are poled alternately. 

The poling process utilizes a high electric held to initialize a polarization direction in 
a ferroelectric material to induce its piezoelectric properties. Poling is accomplished by 
application of 350 V at room temperature for two to three seconds. The convention for 
coordinate systems describing piezoelectric phenomena is to call the axes 1,2, and 3, where 
the poling direction is taken as the 3-direction and the 1 and 2 directions are mutually 
orthogonal. After poling, the electrodes are silver-painted together and wires are attached 
to top and bottom electrodes. 

Later, when a drive voltage is applied, neighboring segments which are poled up and 
down, expand and contract respectively in the 1-direction (using d^i motor action), causing 
the structure to buckle. Oscillatory drive voltages near the resonant frequency force the 
structure to vibrate with maximum amplitude in standing waves. 


39 




Figure 3-6: A four wavelength traveling wave motor is electroded with six half-wavelength 
alternately-poled segments and four quarter-wavelength spacer segments. 

The piezoelectric elements are designed in a mariner so as to induce the fourth flexure 
mode as was shown in Figure 3-3. To create a four-wavelength traveling wave ultrasonic 
motor, the electrodes are patterned as illustrated in Figure 3-6. As two neighboring elec¬ 
trodes create one wavelength of bending, eight electrodes are required for a four-wavelength 
motor. One pair is left passive and cut in quarter-wavelength and tliree-quarter-wavelengtli 
segments respectively, in order to space the two electrode groups by ninety degrees. The 
two opposing electrode groups are driven out of phase in time, also by ninety degrees, in 
order to induce traveling waves. 

We can investigate the viability of this model by evaluating our four-wavelength motor of 
Figure 3-1, whose dimensions are noted in Figure 3-5(a). For an average radius of r = 3.25 
mm, we find that A = 5.1 mm and L = = 2550 /mi. The elastic modulus of steel, c m , 

is 200 GPa and the elastic modulus of PZT (PTS-1195) is very close to that of aluminum.,, 
c p = 70 GPa [Piezo Systems 85]. The beam is modeled as not including the teeth for 
purposes of calculating the stiffness and location of the neutral axis. The base portion of 
the metal layer is then /?.& = 300 /mi, the PZT layer is h p = 188 /mi, and the depth of the 
beam is b = 3000 /mi. This gives a neutral axis location at g = 294 /mi from the bottom of 
the PZT layer. 

The equivalent stiffness of the composite beam is: 

El — Cplp T “'at G 
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El — c p b j z 2 dz + c m b j z 2 dz 

pzt base 

where x here is measured from the neutral axis. The equivalent stiffness is calculated to be 
El = 3.39 X 10 -3 Nm 2 . 

The total mass per length, with the mass density of steel taken as p me t = 7860 and 
the mass density of PZT as p pzt = 7600 is: 


(P A )tot — (P A )pzt A Pmet (Abase + At 00 th) — 2.32 X 10 


-2 kg 


m 


The resonant frequency of the motor vibrating in the fourth flexure mode is the first 
natural frequency in our model of this pinned-pinned beam of length L: 


El 


xq = 


L 2 V (pA) 


= 5.80 X 10 5 


tot 


rads 


sec 


= 92.3 kHz 


In actuality, the motor runs at several resonant frequencies, one of them being 98.5 kHz. 
Interferometric measurements verified that this was indeed the fourth flexure mode. The 
calculated natural frequency then is within: 


6.2 kHz 
98.5 kHz 


= 6 % 


Our model was based on the assumption of a half-wavelength monomorph acting as a 
Bernoulli-Euler beam where the length of the beam is much greater than the height (usually 
10 or 100 times larger). For our 8 mm motor: 


hp T /q 


= 5.2 


and we find that our beam is far from ideal, so 6% error is not unreasonable. 

We can also check the wavespeed and the rotor speed. The phase velocity of the traveling 
wave along a beam was given as 


c T = y/uZ 


, El m 

4/ = 471 — 


(PA) 


tot 
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and for our motor of radius r = 3.25 mm, 


cr k rads R 

u wave = — = 1.45 X 10 5 - = 1.39 X 10 6 rpm 

r sec 

The rotor no-load speed, however is: 


^rotor,max 


2irijJihw 0 

A r 


where the minus sign signifies that the rotor moves in the reverse direction from the travehng 
wave. 

Interferometric measurements determined that with no rotor, w 0 was 0.1 fj,m when 
excitation was 10 V peak at 98.5 kHz. With h equal to half the total thickness of the beam, 
| (h p + hb), the no-load rotor speed is calculated to be: 

rads 

^rotor,max — 5.4 — 51 rpm 

sec 


Experiments with a small glass rotor, similar in dimension to the brass rotor in the lower 
photograph of Figure 3-1, resulted in rotational velocities of 36 rpm, which is fitting, as the 
rotor has finite inertia. 

It is interesting to compare the final rotor rotational speed to the original excitation 
frequency: 

, i ^ A rads 1 

^rotor s ^ 1 

uq “ 5.8 X 10 5 ^ ~ 107,000 

s 

or to the speed of the traveling wave around the ring: 

/ i ^ A rads 1 

^rotor s ^ J- 

Wwave ~ 1.45 X 10 5 ^ ~ 26,900 

These large numbers signify the magnitude of the rectification of the high-frequency 
excitation, enabling resonant vibration and maximum amplitude deflection of the stator to 
be converted into low-speed rotational velocity of the rotor. Because of this property, ultra¬ 
sonic motors can run many applications (depending on the torque requirements) without 
the need for gears. 

Before we address the issue of torque production and the mechanisms behind the fric¬ 
tional coupling of vibrational energy of the stator into rotational kinetic energy of the rotor, 
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Figure 3-7: (a) One half-wavelength section of the stator is modeled as a simply-supported beam, 
(b) With polarization and electric Helds both in the 3-direction, the monomorph will bow upwards 
as shown. 

let us first look at the conversion process of electrical energy to strain energy of the flexing 
stator. This will enable us to predict the amount of deformation, w 0 , for a given drive 
voltage. 

3.6 Deflection of the Stator 

In order to investigate the dynamic deflection of the beam, we begin with an analysis of the 
static deflection of the stator due to piezoelectric excitation. We examine a half-wavelength 
portion of the stator modeled as a simply-supported beam as illustrated in Figure 3-7(a). 
For now, we ignore the teeth on the stator and focus on the monomorpli which consists of a 
bottom layer of piezoceramic of thickness h p , bonded onto the base metal layer of thickness 
/?.&. Because we model the teeth as contributing only added mass to the structure (and 
lienee mechanical amplification) and not stiffness, only the piezoceramic and base metal 
layers affect the final bending configuration. 

After determining the static deflection of the beam, the deflection magnitude at reso¬ 
nance can be calculated by multiplying the static value by a dynamic amplification factor. 
The dynamic amplification factor of the composite monomorph structure is the quality fac¬ 
tor of the piezoceramic material weighted by the strain energies of the two material layers 
[Nasliif et al. 85] as will be discussed later. First however, we find the static deflection. 

For a piezoceramic poled in the positive ^-direction and driven with an electric field 
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also in the ^-direction, the piezoceramic will expand in the ^-direction and contract in 
the ^-direction causing the structure to bow upwards as shown in Figure 3-7(b). Since the 
common convention defining a positive bending moment is for the beam to bow downwards, 
the effect of this polarity drive is to induce a stress held which creates a negative bending 
moment throughout the structure. 

In static equilibrium, force and moment balance equations lead to determination of the 
curvature in terms of stresses and strains, which can then be related to deflection of a beam, 
given appropriate boundary conditions. [Lucas 75] describes curvature and deflection of a 
monomorph cantilever beam and we follow the notation here, but apply our analysis to a 
simply-supported beam where we are interested in the magnitude of the deflection at the 
midpoint. 

Figure 3-8 illustrates the ^-directed stress distribution, T\ (z), acting on the avfaces 
of the beam (T\ (z) is the reduced matrix notation for T xx (z)) where the beam has been 
subjected to bending with a radius of curvature, p K . Along the unstrained neutral axis of 
length ds, 

p K d(j) = ds 

For small deflections, 

d(f) d(f) 1 d 2 w 
ds dx pk dx 2 

Here, dx is the length of the neutral axis in the initial undeformed state and the curva¬ 
ture, -j-, is equal to the second derivative of the deflection profile, w(x), of the beam. 

A longitudinal differential line element at a distance z — g from the neutral axis location, 
g, exhibits extensional strain of the form: 


where 61 (z) is the ^-directed strain component acting on the face normal to the avaxis 
(Si (z) = S xx (z)). 

Moment balance requires that the stress distribution, T\ (z), multiplied by the moment 
arm and integrated over the cross-sectional area of the beam be zero, as there are no 
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Figure 3-8: The radius of curvature, p K , due to bending moments caused by contraction of the 
piezoceramic can be determined from the resulting stresses and strains in the composite beam. 

externally applied moments or forces acting on this structure: 


M = 0 = b j (z — g)T\ (z) dz 

Z 

Next, we use the stress-strain constitutive law to restate this moment balance expression 
in terms of the curvature. Stress, strain and applied field in the monomorpli are related by 
the piezoelectric strain relations (described in Appendix B): 

S = s s T + d t E 

D = dT + e r E 


where 


S = 

Strain tensor 

T = 

Stress tensor 

D = 

Electric flux density vector 

E = 

Electric field intensity vector 

s B ' = 

Elastic compliance tensor at zero E 

d = 

Piezoelectric strain tensor 

dt = 

Piezoelectric strain tensor transposed 

c T = 

Dielectric permittivity tensor at zero T 


Alternatively, the constitutive relationships between stress, strain and applied field can 
also be expressed in terms of the piezoelectric stress relations: 


T = c E S - e t E 
D = eS + e s E 
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where 


c E = Elastic stiffness tensor at zero E 
e = Piezoelectric stress tensor 

et = Piezoelectric stress tensor transposed 

e s = Dielectric permittivity tensor at zero S 

In the case of the monomorph presented here, the faces normal to the 2- and 3-directions 
are stress-free surfaces, and so: 

t 2 = t 3 = 0 

The stress distributions on the 1-faces due to an applied electric held in the 3-direction are: 

T\ (z) = c (z) Si (z) - e 31 (z) E 3 

The piezoelectric stress tensor, e, is related to the piezoelectric strain tensor and the 
stiffness tensor by: 

e = dc 


giving: 


e 3 i (z) = d 3 1 (z)c(z) 


This results in a stress distribution of the form: 


T\ (z) = c (z) Si (z) - d 3 1 (z) c (z) E 3 
where c (z) and d 3 i (z) are step functions at the materials interface: 

f Cm hp ^ Z fop T fob 


C z = 


0 < z < h v 


d 3 i (z) = 


0 


fop < z < h p + fob 


[ d 3 i 0 < z < hp 
Substituting the expression for (z) into Ti (z) gives: 


Ti (z) = c (z) ( ^ ) - d 3 1 (z) c (z) E 3 
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Plugging this stress expression into the moment balance equation yields: 


M = 0 = 


b J 

Z 


c (z 


( z ~ 9? 

Pn 


(* 


g) d 3 1 (z) c (z) E 3 


dz 


Solving for the curvature and noting that the second term in the integral exists only in 
the piezoceramic layer, we find: 

\ 

hp 

I (z- g )dz 

_ 0 _ 

hp h p -\-hb 

f c p (z - g) 2 dz+ f c m (z - gf dz 

0 hp 

Integration produces an expression for the curvature in terms of thickness and compli¬ 
ance ratios: 

1 _ 3 d 3 \E 3 1 
Pk 2 hp r K 

where ■p- is a dimensionless curvature parameter, 

1 = _ (1 ~ 2 fp) _ 

r K 1 - 3 f 0 + 3 / 0 2 + m(3p + 3 p 2 + p 3 - 6 pf 0 - 3 p 2 f 0 + 3p/ 0 2 ) 


— — d 3 \c p E 3 < 

Pn 


and 


m = 


h b 

P = T 

hp 

g 1 + mp 2 + ‘Imp 
0 hp 2(1 + mp) 

We can solve for the deflection profile, w(x), from the curvature-displacement relation¬ 
ship, 

1 d 2 w(x) 3 d 3 \E 3 1 
p K dx 2 2 h p r n 

by taking into account the boundary conditions of a simply-supported beam: 

/ x x d 2 w . . d 2 w . 

w (x = 0) = 0 w (x = L) = —— (x = 0) = —~pr (x = L) = 0 

dx 1 dx 1 
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to determine the constants of integration. The displacement profile becomes: 


3 d 3i E 3 1 
wix) = — 


4 h p r K 


^x 2 — LxJ 


and the resulting deflection at the midpoint is: 


W r , 


X = 


3 d 3 \E 3 L 2 
16 h p r n 


As the electric held is equal to the applied voltage, V 3 , divided by the thickness of the 
piezoceramic, we now have an expression for the static deflection of the beam in terms of 
the applied voltage: 

f _L\ _ 3 d 31 V 3 L 2 

Wmax { X 2) 16 h 2 p r K 

Note that d 3 \ is a negative quantity and so the resulting deflection is positive, as ex¬ 
pected. What we are really interested in however, is the magnitude of the dynamic deflection 
of the stator at resonance. At resonance, the system is damping controlled and dependent 
on the losses in the materials. The resulting dynamic deflection is equal to the static de¬ 
flection of the beam times the dynamic amplification factor, or quality factor, Q. Damping 
characteristics of a system are often also expressed in terms of the loss factor, r/, where: 


Q 


1 


For a composite structure such as the monomorph presented here, the loss factor of the 
entire structure can be expressed in terms of the loss factor of the piezoceramic material 
times the ratio of the strain energy of the piezoceramic layer to the strain energy of the 
entire monomorph [Nashif et al. 85]: 


— h'pi 


u. 


piezo,strain 


u. 


mono,strain 


where the strain energy stored in the flexural mode, U s t ra im is defined as the integral over 
the volume of the strain energy density: 


lb f 

u stram = — J S 1 (z) T\ (z) dz 
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Substituting in the expressions for stress and strain for flexural energy alone, gives: 


U strain = f C (z) (- -- ) dz 

^ J \ Pn 


Integration over the appropriate layers results in a stored elastic energy (in Joules) in the 
piezoceramic layer of: 


U, 


= -Lbh-ndliEnC, 


piezo,strain — 0 - LjU,t p Uj 3\ J - J 3 { ^p~ 
o 


(1 - 3 f 0 + 3/ 0 2 ) 


and in the entire monomorph of: 


U, 


_ 3 rrr ;2 r 2„ (1 - 2/ 0 ) 

mono,strain — L/Oilpd^ -C /3 Cp 

o T K 


The loss factor for the device is then calculated to be: 


Pmono 


Ppiezo 


1-3 f 0 + 3/p 
r K (l-2f 0 ) 


and the dynamic amplihcation factor: 


Qmono — Qt 


(1 ~ 2 f 0 ) 
1-3 f 0 + 3/ 0 2 


For example, we can calculate the expected deflection for our motors of Figure 3-5, 
ignoring the teeth and using the dimensions and elasticities given in Section 3.1. We find 
m = 2.86, p = 1.60 and f 0 = 1.56, resulting in a value for the dimensionless curvature 
parameter of: 

— = -0.291 

r K 

One experiment was performed earlier using an interferometer [Zhang, Pan and Cross 
88 ] to measure the deflection of the stator at the fourth flexure mode resonance under an 
applied voltage of 10 V. (This is a rather low held, = 5.32 X 10 4 11, whereas the maximum 
depoling held for the ceramic is E c = 2.36 X lO 6 ^ [Piezo Systems 85]). For such a drive, 
we would expect a dehection magnitude of: 


l^max 


X 


L 

2 


= 1.87 X 10 8 m 
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The actual measured value at resonance was 0.1 /im, which means there was a dynamic 
amplification factor of the monomorph of 5.4. 

For a 10 V drive, the total stored elastic energy in the piezo layer is: 


U, 


piezo,strain 


= 1.24 x io -9 / 


For the entire monomorph, the total strain energy turns out to be: 


U, 


mono,strain 


= 6.47 X 10“ 9 / 


The quality factor for the monomorph device then is: 


Q 


mono 


6.47 

T24 


Qr 


5.20 Q p 


The loss factor for the PTS-1195 piezoceramic is given as i] p i ezo = 0.015 [Piezo Systems 85], 
yielding a Q p of 67, and a consequent amplification factor of: 


Q 


mono 


350 
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Chapter 4 


Electrical Modeling of the Stator 


We can model the stator both mechanically, as a beam vibrating, and electrically, as an 
equivalent circuit looking into the electrical terminals. Electrode groups A and B of Fig¬ 
ure 3-6 each act to produce standing waves when excited near their resonant frequencies. 
At these frequencies, the stator can be modeled by the electrical equivalent circuits shown 
in Figure 4-1. The series -RLC branches are the mechanical arms. The resistive, inductive 
and capacitive lumped parameters are the mechanical damping, mass and compliance, re¬ 
spectively, transformed into electrical equivalents due to the piezoelectric effect. C 0 is the 
electrical capacitance in the absence of the mechanical motion which is induced at reso¬ 
nance. At resonance, the inductive and capacitive components of the mechanical arm, juL 
and become equal in magnitude and opposite in phase, leaving the system damping 
controlled. This resonant frequency, f r , shown in Figure 4-2, corresponds to an impedance 
zero. Above resonance, the series RLC branch becomes inductive and the equivalent circuit 
reduces to a parallel resonant circuit, displaying an antiresonance at the frequency f a , cor¬ 
responding to an impedance maximum. At even larger frequencies, the system again acts 
capacitively. 

The resonance and anti-resonance frequencies of a single electroded slab of piezoceramic 
(as opposed to the monomorph device used in the ultrasonic motor here) can be related 
to material properties and the electromechanical coupling coefficient (see Appendix B). 
Laboratory measurements of the impedance characteristics of a piezoelectric material can 
then be used to characterize and verify the electromechanical coupling factor, the damping 
coefficient and the quality factor of the device. 
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Sector A Sector B 


Figure 4-1: The electrical circuit model of an ultrasonic motor consists of two series-resonant, 
circuits, one for each standing wave. 
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Figure 4-2: Near a resonant frequency, a piezoelectric material has an impedance characteristic 
displaying both a resonance and antiresonance. This small blip can be seen on the \Z\ curve near 
26 kHz. The phase, 9 and electrical damping factor, Q, also change significantly at this point. 
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The electromechanical coupling factor for a piezoelectric material where a held is applied 
in the 3-direction and the output strain is taken in the 1-direction, is k^i. [Buchanan 86 ] 
gives &31 as: 


&31 = 


2 P an {2-7; 


d- 


31 


\ l + f 


fa 


tan tt 


. fg — fr 


f T 

fc 33*ll 


fr V 2 fr 

However, for a monomorph, the 1-directed strain must work against the stiffness of the 
non-piezoelectric steel layer, with the result being bending instead of lateral displacement. 
The effective electromechanical coupling factor for a monomorph then would be lower than 
for a plain slab of piezoceramic. For instance, [Smits and Choi 92] showed that for bimorphs 
used as cantilever beams (two electroded slabs of piezoceramic bonded together rather 
than one layer of piezoceramic bonded onto one passive layer), the squares of the effective 
electromechanical coupling factors, A e jj, when the tip works against a vertical moment, 
force, or distributed load, are respectively: 


Am = 3 (^ 


9 ( ^31 

4 U -kl 


A p 


5 / ^1 \ 

3 1 4 — fe^y 


For our ultrasonic motor, we do not have cantilever beams. Our monomorph structure 
can, however, be modeled as a simply-supported beam of length L = 

The effective electromechanical coupling factor for our simply-supported monomorph, 
which we denote as T e jj, will be different than the effective electromechanical coupling 
factor for a cantilevered bimorph device. For a simply-supported monomorph working 
against moments applied at the ends (as in the case of a traveling wave motor), we can 
expect: 


T M < 3 


ft 31 


4 fell. 


since a single-layer monomorph structure will be less effective than a bimorph in delivering 
output moments and forces. 

The electromechanical coupling coefficient, & 31 , for our PZT thin films was measured to 
be 0.22 [Udayakumar 91], and is given as 0.34 for the PTS-1195 ceramic in the motors of 
Figure 3-1 [Piezo Systems 85] which leads to the result that the percentage of mechanical 
energy converted from input electrical energy for our bulk-ceramic motors will be on the 
order of 11.6%. This does not imply efficiency however, as the remaining 88.4% is not 
dissipated but stored dielectrically. 
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Whereas energy conversion is described by k%i in a slab of piezoelectric material, or by 
A e ff and T e ff for devices such as bimorphs and monomorphs, energy loss in a piezoelectric 
material is denoted by the quality factors Qm (mechanical) and Qe (electrical, or dielectric). 

The equivalent circuit of Figure 4-1 can lead to determination of the mechanical quality 
factor of a piezoceramic, Qm , which is the ratio of average energy stored per cycle to average 
energy dissipated per cycle. [Jaffe 58] gives an expression in terms of electrical equivalent 
circuit parameters: 

— = 47tA/ | Z m m | (C 0 + Ca ) 

VM 

The electrical quality factor is related to the dielectric dissipation or loss tangent, tanS , 
the form usually specified with the material. tanS is the ratio of the imaginary part of the 
complex permittivity to the real part: 


—— = tanS 

Qe 


Im {e} 
Re {e} 


This electrical model of each stator section is applicable only in the vicinity of an isolated 
resonance. Unwanted excitation of nearby modes can result in mode coupling and unwanted 
forms of deformation. 
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Chapter 5 


Energy Conversion 


To this point, we have developed a model of the mechanical and electrical characteristics 
of the stator portion of an ultrasonic motor, from the input electrical terminals to the 
output displacements on the stator surface. It would be useful to compare such energy 
conversion processes across a spectrum of actuator technologies in order to expose relative 
merits among the possible choices. 

[Hollerbach, Hunter and Ballantyne 91], in their comparative analysis of actuator tech¬ 
nologies for robotic applications, point to stall torque density as a figure of merit amongst 
disparate technologies ranging from electromagnetic motors through pneumatics, hydraulics, 
magnetostrictors, shape memory alloys and piezoelectric inch worm motors. They claim that 
stall torque density is a better figure of merit than power density, as power density calcula¬ 
tions require the proficiency of the drive electronics to be taken into account. [Bart et al. 88] 
alternatively compare electromagnetic and electrostatic motors based on energy densities 
storable in the gap between the stator and the rotor. Here, we extend the latter analysis to 
energy densities storable in piezoelectric ultrasonic motors. 

In an electromagnetic motor, the storable energy density, u mag , in the air gap between 
the stator and the rotor is: 

_ IB 2 
Umag ~ 2 

where B is the magnetic flux density and is the permeability of free space. 

In an electrostatic motor, the storable energy density in the air gap is: 


■ — — F 2 

^electrostatic — ^ 
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where E is the electric held intensity and e 0 is the permittivity of free space. 

[Bart et al. 88] note that the maximum flux density available in an electromagnetic 
motor at saturation is approximately 1.5 T, yielding a maximum energy density of: 


^mag ,max 


= 9 x 10 5 


J 

m 3 


Similarly, for an electrostatic motor, the maximum electric held before breakdown is on 
the order of 10 6 —, which gives the maximum energy density attainable as: 


^electro static,max 



For small gaps on the order of 1 fj,m, achievable with silicon lithographic processes, [Bart 
et al. 88] also point out that electric helds on the order of 10 8 ^ are sustainable, leading to 
energy densities commensurate with those of electromagnetic motors: 


l^thin electrostatic,max 


= 4.4 x io 4 


m 3 


Torques however, are computed from the spatial derivatives (the changes in rotor po¬ 
sitions as the rotors slide past the stators) of the total coenergies, U' mag and U' electrostatic . 
For a linear material (i.e. an air gap), the energy, U, is equal to the coenergy, U'. Total 
coenergy in the air gap of an electromagnetic motor then is the integral over the volume of 
the energy density: 


U' 

mag 


vol 


-—dv = -Li 

2 jJ,Q 


where dv is a differential volume element, L is the inductance and i is the current. 
Total coenergy in an electrostatic motor becomes: 


^electrostatic = j \^E 2 dv = ^ CV 2 
vol 

where C is the capacitance and V is the voltage. 

The torque for an electromagnetic motor can then be calculated by taking the derivative 
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of the coenergy with respect to the rotor position, 6 , keeping i constant: 


^~mag — 


^^mag i 

86 1 


1 dL , 2 

2le l 


Similarly, for an electrostatic motor, the torque produced is calculated with V held 
constant: 


electrostatic 


8U' 

^ electrostatic 


86 


v = 


I^T/2 
2 d6 


If we look now at piezoelectric motors, we find that the storable energy density is 
related to the electric held applied across the piezoceramic, “the gap”. In these capacitive- 
type motors, the gap is no longer air, but the high dielectric material, PZT. The storable 
electrical energy density input to a piezoelectric motor becomes: 


u. 


piezo,input 




pzt^o 




where e pzt is the relative dielectric constant of PZT (1950 for the PTS-1195 ceramic in the 
motors of Figure 3-1 [Piezo Systems 85] and 1300 for our PZT thin films [Udayakumar 
91]). For bulk-ceramic PZT, the electric held breakdown strength is on the order of 10 6 ^-, 
similar to the breakdown held of macroscopic electrostatic motors. Thus we see that the 
storable electrical energy densities of bulk-ceramic piezoelectric motors are three orders of 
magnitude larger than those of electrostatic motors, due to the large dielectric constant. 

The total electrical energy input to a piezoelectric motor is also the integral over the 
volume of the energy density: 


U, 


piezo,input 


— f -tpzt^oE dv 


ol 


In a piezoelectric motor however, torque is not generated as the spatial rate of change of 
total energy, but rather is hrst converted to strain energy in bending and then converted to 
output forces through frictional coupling. The amount of input electrical energy converted 
to strain energy in a slab of piezoceramic is characterized as [Lucas 75] as is described 
in Appendix B. For a monomorph, the ratio is 

For example, we can calculate that the maximum allowable energy density input to our 
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bulk-PZT motors will occur at the breakdown field, 10 6 —: 

’ m, 


3 J 

b^bulk piezo,max input — 8.6 X 10 T 


The output strain energy density is the total strain energy, U p i ez0tS t r ain, divided by the 
volume, Lbh p : 


_ ^ i2i jp 2 

ILpiezo, strain — ~Q^13^- J 3^i 


(1 - 3 f 0 + 3/ 0 2 ) 


which, using the maximum breakdown field strength, calculates to: 


b^bulk piezo,max strain — 6.65 


J 


mr 


The ratio of the latter, strain energy density, to the former, input electrical energy 
density, is T 2 e jp. 


T 


2 _ 
eff ~ 


6 . 65^3 

_ m^ 

8.6 X 10 3 ^3 

m^ 


= 7.7 X 10“ 4 


or 

T eff = 0.028 


This calculated value is roughly a factor of 12 lower than that calculated from the given 
(measured) value of ko,\ = 0.34: 

Tm < 3 (j&r) = 0.090 

When macroscopic electrostatic motors are shrunk to microscopic sizes with gap widths 
on the order of 1 /rm, electric field breakdown strength increases by two orders of magnitude. 
This is due to the fact that as the gap gets smaller or the pressure in the gap decreases, 
there is a decreased chance of air molecules colliding, leading to breakdown. 

For thin film dielectrics, it is also the case that electric field breakdown strengths see 
a two-order of magnitude improvement from the roughly 10 6 ^) available for PTS-1195 ce¬ 
ramic [Piezo Systems 85] which is 188 /im thick, to the 10 8 ^ possible in PZT thin films 
[Udayakumar 91] which are 0.3 /im thick. In thin films, the reason for improved break¬ 
down strength is not a decrease in the pressure-gap product, but rather a decrease in the 
likelihood of point defects due to thin layers. 
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Figure 5-1: The breakdown strength of PZT increases dramatically as the material moves from 
bulk to thin-film form. The round dots indicate data from [Gerson and Marshall 59] while the square, 
block signifies new data on thin films [Udayakumar 91]. 


Studies of breakdown strength vs. thickness in bulk PZT samples were first conducted by 
Gerson and Marshall [Gerson and Marshall 59]. They found that the relationship empirically 
fit the curve which is illustrated in Figure 5-1 where: 


E bd = 27.2/7- 0 ' 39 


Here, E bd is the electric field intensity at breakdown, given in and h v is the thickness 
of the piezoceramic layer, in centimeters. 

Although, at the time, Gerson and Marshall had no tliin-film PZT, our recent fabrica¬ 
tions runs have shown that thin film breakdown measurements do fit precisely with their 
earlier work. 

We can see then, that while electrostatic micromotors display improved energy densi¬ 
ties over electrostatic macromotors, microfabricated piezoelectric motors exhibit the same 
improvement in breakdown strength, but again contribute three orders of magnitude im¬ 
provement in energy density due to the large dielectric constant. The maximum energy 
density then, for piezoelectric micromotors is: 


^ Ih in piezo,max input 


— , ^pzt^oEf j — 5.8 X 10" „ 

2 


Plugging in the above expression for E bd , we derive an expression for energy density 
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Figure 5-2: Energy density and total energy, for a given active electrode area, depend on piezoce¬ 
ramic thickness. 

as a function of piezoelectric thickness which is plotted in the top graph in Figure 5-2. 
Converting to mks units for h p in meters: 

Uptezopnput = (4-26 X 10 4 ) ((100fc p p 0 ’ 78 ) 

Multiplying by the area and thickness, we find the total energy as: 

U ptezo , tnp ut = (4.26 X 10 4 ) ((lOOM" 0 - 78 ) U){h p )J 

This expression of total energy as a function of piezo-layer thickness is illustrated in the 
lower graph in Figure 5-2 for an area equal to the 8 mm outer diameter, 5 inner diameter 
motor of Figure 3-1. For a him thickness of 0.3 /mi, we find the total input electrical energy 
for a thin-film PZT version of the motor in Figure 3-1 to be: 


Uthi 


n piezo , input 


= 5.3 x 10“ 3 J 


For the same size motor, but using bulk PZT 188 /mi thick, we find 


Ubulk piezo,input — 2.1 X 10 J 
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Although the total electrical energy storable in this bulk PZT motor is four times greater 
than for a thin-hlm motor of equivalent area, we are interested in the possibilities of lower 
manufacturing costs which microfabrication holds. If motors can be printed rather than 
machined, sanded, grinded and soldered, then we can hope to develop small, compact, cheap 
actuators. In addition, moving away from machine tools and towards microfabrication 
and lithographic manufacturing processes provides a means for scaling to even smaller 
dimensions or patterning finer features. 
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Chapter 6 


Contact Mechanics 


We have described the energy conversion process which takes place in the stator of an ultra¬ 
sonic motor. That is, we have discussed the first stage of the energy transduction mechanism 
of a vibromotor which we outlined in Figure 1-1. We have seen how input electrical energy 
is transformed into output strain energy in the form of mechanical displacements, where 
the percentage of input electrical energy converted to strain energy is characterized by T 2 e ^ 
for a monomorph. Losses in this resonant motor were described by the mechanical quality 
factor, Qm- 

What we are really interested in, however, is how much power can be delivered to a load 
and what form the speed-torque curve will take. We would also like to know the overall 
system efficiency and the stall torque density. To answer these questions, we focus on the 
second stage of the energy conversion process, the vibration rectification achieved through 
frictional coupling between the rotor and the stator. 

Actuators utilizing these mechanisms form a much wider class than the various embodi¬ 
ments of ultrasonic traveling wave motors. Vibratory feeders transporting granular particles 
on a bed in process-control applications impart a flight phase or sliding phase to the par¬ 
ticles to be moved by inducing an oscillatory diagonal motion in the bed [Gaberson 72]. 
Early piezoelectric actuators used ratchet-style mechanisms, wherein a single-point contact 
would oscillate and impact against the rotor at an oblique angle, creating a uni-directional 
drive [Barth 73]. Some modeling of these types of systems is described in [Sashida 82] 
and [Kobrinskii and Lennox-Napier 69]. In all of these actuators, the stator is in contact 
with the rotor for one portion of the cycle pushing it in one direction, and is retracted 
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from contact in the remaining portion of the cycle when the stator point is moving in the 
opposite direction. A primary disadvantage of this type of actuator is that the stator must 
be adjusted and held fixed to within a few microns of the rotor when at rest. 

One subset of vibromotors, in which the stator and rotor always maintain frictional 
contact are the hammer-type actuators described in the work of [Higuchi et al. 89], which 
quasistatically inflict an impulse to a mass which produces a reaction force larger than the 
frictional sticking force. 

Resonant traveling-wave ultrasonic motors which were first invented in the Soviet Union 
in the sixties [Ragulskis et al. 88] and again in Japan in the eighties [Sashida 83] contributed 
a means for creating larger amplitude displacements, maintaining continuous rotor-stator 
contact and allowing for electronic control of bi-directional motion. While much work has 
been done on modeling the stator energy conversion processes, far less has been reported 
on modeling the frictional coupling at the rotor-stator interface. 

On the friction topic, [Hosoe 89] and [Okumura and Mukohjima 87] have reported on 
wear problems and efforts at improving efficiency in Canon ultrasonic motors used in auto¬ 
focus lenses. [Maeno, Tsukimoto and Miyake 92] have simulated the slip-stick interaction of 
the Canon motors using finite element analysis and [Maeno and Bogy 92] have investigated 
the hydrodynamic effects on the coupling. Other studies have been undertaken in Germany 
to understand the tribological properties of ultrasonic motors [Rehbein and Heinz 92] and 
experiments with linear ultrasonic motors using a variety of frictional interface materials 
have been reported by [Endo and Sasaki 87]. While many papers cite output torques and 
speeds of prototype motors with efficiencies in the range of 35% to 45%, [Kumada 85] claims 
the highest efficiencies, on the order of 80%. 

[Minotti and Lallement 93] have modeled the rotor-stator interaction parameterizing the 
contact conditions of a non-rigid rotor in terms of angle of contact. [Kurosawa and Ueha 
88 ] have performed closed-form analyses of interface losses, claiming theoretical limits of 
70% efficiency for traveling wave motors and [Hirata and Ueha 93] have described lumped 
parameter equivalent circuits of the frictional coupling. 

Our analysis of this process can be divided into several steps. First, we focus on the 
contact mechanics of a number of simple models of the interface conditions between the 
rotor and the stator, ignoring the tribological issues associated with surface asperities, 
hardness, wear and other surface effects. Later, we will make some comments on experiments 
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Figure 6-1: (a) We need to determine the form of the speed-torque curve. For each applied normal 
force, Fjv, and applied voltage, V r , there will be one curve characterizing the relationship of speed 
to torque. This curve will shift and possibly change shape as V or Fjv is increased, (b) For a given 
applied voltage, we are interested in determining how the curve shifts with F^. Parameterizing 
each curve by ce no _/ oa( j, we can find a relationship between ce no _/ oa( j and F^. (c) Alternatively, the 
performance of the motors can be described by parameterizing each curve by r sta ;; and plotting 
Tstaii vs. Fn- By calculating these curves, we will be able to apply an optimal normal force to these 
motors such that the power output will be maximized for a given applied voltage. 
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performed with rotors exhibiting different asperity magnitudes. 

The goal is to find a figure of merit for ultrasonic motors in order to compare them 
to other types of actuators. We have chosen stall-torque density as the metric for judging 
disparate technologies, and Figure 6-1 illustrates the approach. The torque produced by 
an ultrasonic motor is determined by the friction force at the rotor-stator interface. For 
a given excitation voltage and applied normal force, there will be a set of speed-torque 
operating points at which the motor will run, depending on the load it must drive. With no 
load, the motor will spin at its no-load speed and when the load is increased to the point 
that the rotor no longer moves, the motor is delivering its stall torque. The no-load speed 
and stall torque provide two operating points on the speed-torque curve as illustrated in 
Figure 6-l(a), but we would also like to know the relationship of speed to torque for the set 
of points determining the form of drop-off in between. 

As the voltage or normal force is varied, the characteristic speed-torque curve will shift 
in some manner. By setting the drive voltage to the maximum value before depoling occurs 
and varying the applied normal force, we can determine how the speed-torque curve shifts 
in relation to a new normal force. Figure 6-l(b) illustrates a succession of characteristic 
no-load speeds for a set of applied normal forces and a given excitation voltage. We would 
imagine that if there is no normal force (the rotor does not touch the stator), there will be 
no rotor speed. When a normal force presses the rotor against the stator, without a side 
load, the rotor will spin at a certain no-load speed. As the normal force is increased, the 
no-load speed wiU change and for some very large normal force, the rotor will not spin at 
all. For some value of normal force, Fjv, there will be a maximum no-load speed. Of course, 
for every speed-torque curve characterized by a no-load speed, there will be a related stall 
torque, and the variation of stall torque with normal force can then also be determined as 
shown in Figure 6-l(c). Once we find the maximum stall torque, we can divide by the mass 
or volume of the motor to determine our figure of merit, the maximum stall-torque density. 

The form the speed-torque curves take will depend on two things: 1) our model of the 
deformation, and 2) the friction law we choose. 

In terms of the friction law, a coefficient of friction describes the relationship between 
two surfaces and how forces are transmitted across their area of contact. To speak of “the 
coefficient of friction of a material” is meaningless, as a coefficient of friction describes not 
a property of a lone material, but the relationship between two. In general, the coefficient 
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Figure 6-2: If the rotor contacts the stator only at the apex of the traveling wave, the rotor velocity 
will be equal to the horizontal velocity, Vh,max, of the stator surface as long as there is no slipping. 

of friction of two materials, /<, is not a constant. In Coulomb friction, p, is a constant and 
is multivalued for dynamic and static scenarios. In rolling friction, such as a plate on a 
bed of rollers, the coefficient of friction is proportional to the difference in velocities of the 
plate and the rollers [Bliushan and Gupta 91]. Other, more complicated, friction laws have 
been proposed to include time histories and hysteresis for phenomena such as earthquakes 
[Linker and Dietericlr 92] and robot joints [Dupont and Bapna 92]. 

In pursuit of characteristic speed-torque curves, our approach will be to analyze a series 
of progressively more sophisticated contact mechanisms, along with several different friction 
laws. We will examine three cases: 

• A rigid rotor and a rigid stator contacting along a line, assuming Coulomb friction; 

F t = fxFff, 

• A compliant rotor and a rigid stator contacting over an area, also assuming Coulomb 
friction; Ft = (xFn- 

• A compliant rotor and rigid stator contacting over an area, but assuming viscous 
friction; F T = -£ [v sta tor i x ) ~ v rotor ] %• 

6.1 Line Contact, Coulomb Friction 
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Figure 6-3: Free-body diagram of the rotor and stator with the stator modeled as having only a 
line contact by the rotor. 


The first case is illustrated in Figure 6-2, where we assume the rotor is rigid and contacts 
the stator only at the peaks of the traveling wave, as the stator bending profile does not 
distort. The contact between the rotor and the stator then is along a fine of length b into 
the beam. Due to our assumption of the beam having a uniform cross-section, we will use 
two-dimensional free-body diagrams to depict forces acting on the system. These forces 
are, however, actually forces per unit length. This ideal case of line contact will lead us to 
upper bounds for the no-load speed and the stall torque of an ultrasonic motor. 

If we assume Coulomb friction and no slipping, then at the line of contact the velocity 
of the stator will be equal to the velocity of the rotor, as the stator contact point pulls the 
rotor along. 

Figure 6-3 illustrates the forces acting on the rotor and the stator. For line contact 
between the rotor and stator, the normal force acting on the rotor is transmitted as a force, 
Ra, acting solely on the apex of the stator (Ra = Fn)- While in the condition of no slip, 
when a load or braking force, Fb, is applied to the motor at a level less than that required 
to stall the motor, then at steady-state, the frictional force, Ft, delivered by the rotor will 
increase to match it, being equal and opposite to Fb, and the velocity of the rotor will be 
equal to the velocity of the stator. 

There can be slipping, however, when larger loads are applied. The largest load that 
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^rotor 



Trotor 


Figure 6-4: For point contact and Coulomb friction, rigid rotor and undistorted stator, the speed- 
torque curve shows constant rotor speed for any torque in the region of no slipping. For larger loads, 
slipping occurs and the speed-torque curve is vertical, displaying a constant, maximum torque for 
any rotor speed. 


can be applied without slipping is equal to: 


F t = n a F N 

where /.i s is the static coefficient of friction. 

After sliding is initiated, the friction force becomes: 

Ft = l-id.FN 

where ^ is the dynamic coefficient of friction. Figure 6-4 illustrates the speed-torque curve 
for this ideal case of line contact between a rigid rotor and rigid stator. 

The velocity of the stator contact point is determined by examining Figure 6-2. In this 
figure, the £-( coordinate system is centered at a point P on the surface of the originally 
undeformed stator. The thin-dashed horizontal lines in the figure represent the top and 
bottom surfaces of the undeformed stator. From Appendix A, we saw that when a traveling 
wave of bending is induced, a stator surface point, P, moves with a transverse displacement 
of: 

£ = w 0 cos (kx — ujt) 
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Figure 6-5: (a) The velocity profile for a counter-clockwise elliptical trajectory of a point on 

the surface of the stator is everywhere at right angles to the displacement, (b) The longitudinal 
components of velocity point in opposite directions for ^ > 0 and ^ < 0. 


and with a longitudinal displacement of: 


c 


2k h/w 0 
A 


sin (kx — 


LOt) 


tracing out an elliptical trajectory which satisfies the equation: 

c 2 e 2 

—--- + — = 1 

^ 2whw n ^^ Wq 

The displacement and velocity profiles of a stator surface point undergoing this elliptical 
trajectory are shown in Figure 6-5 (a). The longitudinal speed of a surface point is the 
derivative of the horizontal displacement with respect to time: 



'2iujjhw 0 

A 


cos(kx — ujt) 


The horizontal components, ly,, of each velocity vector in Figure 6-5(aj are shown in Fig¬ 
ure 6-5(b). 

For the case of a rigid rotor and rigid stator, the rotor contacts the apex of the traveling 
wave, giving a rotor speed of the maximum horizontal velocity of the surface point: 

27 Tujh/w 0 

Uh,max — 7 


We see from the speed-torque curve of Figure 6-4 that the rotor speed will be constant 
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Figure 6-6: (a) For a rigid rotor and rigid stator, the no-load speed is independent of applied 
normal force, (b) For this case of line contact, the stall torque can be increased arbitrarily by 
applying larger normal forces. 


at the value: 


Vdrotor 


Vh ,max 

r 


for all load torques within the range of no slipping. Once slipping occurs, the torque 
delivered by the motor remains constant at: 


Trot or — I 1 d l \ !' 

This describes the motoring action of quadrant I of Figure 6-4. As Fb is increased beyond 
this point, the motor slows down while still delivering torque equal to /.i^F^r, until finally 
the rotor stops moving. Increasing Fb further causes the rotor to spin in the reverse 
direction as illustrated by quadrant IV behavior. When Fb is directly oppositely to the 
way it is defined in Figure 6-3, the friction force acts to oppose it and so the rotor and 
stator velocities remain equal and in the same direction. The rotor torque is thus negative 
in quadrant II and remains equal and opposite to the pulling force, —Tg?’, until slipping 
occurs, at which point the rotor velocity increases and the rotor delivers a constant torque 
of: 

7~rotor — “h d.FNV 

Once the speed-torque curve has been characterized, we can investigate how it shifts 
with Fjv for a fixed excitation voltage. As the rotor speed is independent of iqy and the 
torque is directly proportional to iqy, increasing the applied normal force will increase the 
stall torque of the motor without affecting its speed before stall. Figure 6-6 shows the 
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no-load speed and stall torque dependencies on F Xr. For this case, it seems that the stall 
torque can be increased by arbitrarily increasing Fjv- However it is surely the case that the 
rotor and stator will deform at some point and the case of point contact will no longer hold. 
This will lead us to the next case, that of a compliant rotor and rigid stator. 

However, continuing with the point contact condition, we note that the mechanical 
power output by a motor is the product of output speed and drive torque: 

Pmech — Fotor^rotor 

Mechanical output power is plotted as a function of rotor torque in Figure 6-7 (a). Since 
the mechanical output power is P mec h = T rotor^rotors and urotor is constant for torques from 
zero through stall, the output power grows linearly with torque until stall. 

In this case, the maximum output power would be delivered at operating point A in 
Figure 6-4: 

„ {2iruhw 0 \ r ^ , 

p mech,max — l ^ J [J^d^N) 

where it appears that output power is directly proportional to normal force, Fjv- 

The efficiency of this second stage of energy conversion is calculated as the ratio of 
mechanical output power to input power: 

Pmech 

Vm ~ — p 

p in 

The power loss in the system is due to frictional forces at the interface causing heating. 
The power loss is the product of tangential force produced by the rotor times the relative 
velocity of the two surfaces: 


Ploss — Ft (vstator "V rotor') 


or in angular coordinates: 


Plo 


= T. 


rotor 


stator ^rotor ) 


The motor efficiency then, is: 


him 


P 

± rr 


sch 


P 


loss 


p 


mech 


^rotor 


^stator 
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Figure 6-7: (a) The mechanical output power delivered by the motor grows linearly in the region 
of no slipping, from zero torque up through stall, (b) For the case of line contact, the efficiency of 
an ultrasonic motor is 100%) while there is no slipping. 


For line contact, u> stator is a constant and efficiency is linearly related to rotor speed. The. 
relationship of efficiency to rotor torque is illustrated in Figure 6-7(b) where for operating 
points that incur no slipping, that is, up until stall, the rotor speed is equal to the stator 
speed and the motor is perfectly efficient. 


6.2 Line Contact, Viscous Friction 

We can modify the model further by changing the friction law to viscous friction, where 
the friction force is proportional to the relative velocity between the two surfaces times the 
normal force: 

Ft — - \_ Vstator (- r ) Vrotor\ F]\f 

where e is a coefficient of viscous friction with units of seconds per meter. 

For the case of line contact as in Figure 6-2, where the rotor contacts the stator only along 
a line through the width of the motor at the apex of the wave, then v stator ( x ) = v stator,max 
and the rotor no-load speed is: 

_ _ V h ,max 

^no-load. — ^stator,max — 

r 

The net friction force that the rotor applies to the load is equal and opposite to the friction 
force acting on the stator, as in Figure 6-3: 


FT ne t 


—Ft 
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F't — stator,max '^rotorj-^N 

^rotor 



Figure 6-8: For line contact and viscous friction, the speed-torque curve falls off linearly. 


and so the rotor torque is, converting to angular coordinates: 


7~rotor — £ ( [^stator,r 


— U, 


rotor 


\ Fjyr 


Solving for w TO tor-> we see that for viscous friction, the rotor speed falls off linearly with 
increasing torque: 


Ur, 


otor — ^stator ,r 


1 rotor 


eFjsir 2 

The resulting speed-torque curve is shown in Figure 6-8 where the stall torque occurs at: 


7~stall — stator,max 


F N r 2 


The slope of the curve is determined by iqv- Figure 6-9 illustrates how the slopand 
stall torque change as iqy is increased. The rotor no-load speed is unaffected, as we have 
assumed that the stator does not distort and the rotor contacts the stator only along a line 
at the apex. Consequently, the LJ no -i oa d vs. iqy curve and the T sta u vs. iqy curve look very 
similar to those of Figure 6-6. In Figure 6-10(a) we see that uj no _i oa d is constant for all 
values of iqv- Figure 6-10(b) illustrates how stall torque increases linearly with iqy with a 
slope of £LL> s tat 0 r,maxr 2 instead of as in the case of Coulomb friction with line contact. 
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Figure 6-9: For line contact and viscous friction, as the normal force, Fjv, is increased, the stall 
torque shifts out. 



Figure 6-10: (a) For line contact and viscous friction, rotor no-load speed remains constant with 
Fjv. (b) Stall torque increases without bound with F^. 


74 



The output power is again the product of torque and speed: 


Pmech — ^TotoT^Totor 


For this case of viscous friction, the power vs. rotor torque curve is an inverted parabola: 


P, 


mech — 


t-, 9 T ^ stator .max^rotor 

Sr TV' 


where the maximum power point can be found from: 


dP, 


mech 


dr,, 


= 0 


rotor 


and occurs at: 


or: 


- 1 2 
h~rotor Q pt — Sy£^stator,max-T 


7~rotor opt — install 

The rotor speed at r T0i0Topt is: 

_ 1 

^rotor Q pt — stator,max 

The output power vs. rotor torque curve is plotted iu Figure 6 -ll(a). The resulting maxi¬ 
mum power output is: 

Pmech,max — ~^ J ^no — load,h~stall 

The mechanical efficiency of the friction coupling can be found by again calculating the 
power loss due to frictional heating: 


hjm — 


Pmech 


Pmech + Po 


where Pi 0S s is the product of rotor torque and the relative speeds of the sliding surfaces: 


Ploss — h~rotor (pJ stator,max ^rotor ) — 


7~rotor 

sFnP 
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Figure 6-11: (a) Mechanical output power reaches a maximum when the rotor speed is L 0 rotor = 
stator,max and the rotor torque is T rot0 r = b T staii- (b) The efficiency falls off linearly with rotor 
torque from 100% at zero torque to 0% at stall torque. 


yielding a mechanical efficiency of: 


T~rotor 

rim = 1- 

Tstall 


Figure 6-ll(b) displays the efficiency vs. rotor torque curve which falls off linearly with 
rotor torque. 


6.3 Area Contact, Coulomb Friction 

While the line contact assumption was useful in generating theoretical bounds on perfor¬ 
mance of an ultrasonic motor, line contact is an ideal case. We can make our model more 
realistic by relaxing the constraints of a rigid rotor and rigid stator. Here, we analyze the 
contact condition of a rigid stator with a compliant rotor, maintaining the assumption that 
the friction law acts as Coulomb friction. 

6.3.1 Physical Interpretation 

By assuming that the stator is rigid (i.e. its displacement profile is undistorted after contact 
with the rotor), we can say that the rotor conforms to the stator over some finite area and 
the normal pressure at the interface will be proportional to that displacement. Since we 
know the displacement of the stator, we then need to determine the area of contact. 

This contact condition is illustrated in Figure 6-12. We change the coordinate system 
here to make the following discussion more straightforward. Whereas earlier diagrams of 
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Figure 6-12: For the case of a compliant rotor and rigid stator, the rotor undergoes deformation 
such that the contact is over a finite area of the stator. The coordinate system here is centered on 
a peak at the original top surface of the stator. 


the deforming stator were drawn for a transverse displacement of £(.t,/) = w 0 cos(kx — ujt) 
for ujt = here we take ujt = 0 and place the ^-axis on a peak of the stator. The .r-axis is 
placed on the top surface of the undeformed stator. 

From Figure 6-12 we see that a given normal force will cause the rotor to comply and 
interpenetrate the stator surface over an area parameterized by x 0 . While the rotor moves 
with velocity v rotor , the points on the stator surface from —x 0 to x 0 all have different 
horizontal velocities. As was illustrated in Figure 6-5(b), for a traveling wave moving in the 
positive .r-direction, the horizontal components of stator surface point velocities point in the 
negative .r-direction for £ > 0 and in the positive .r-direction for £ < 0. Due to symmetry 
of this problem about the s-axis, we can focus on just the positive .r-axis in describing 
the system. For the moment, we consider the case where the applied normal force is such 
that the interpenetration of the rotor into the stator is such that 0 < x 0 < ^ (i.e. the 
interpenetration depth is less than w 0 ). In this case, stator surface points in contact with 
the rotor all have horizontal velocity components in the negative x- direction, although of 
different magnitudes. 

The operating point of the motor depends on the load applied, Fb- We would like to 
find the characteristic speed-torque curve as in Figure 6-l(a) for this contact condition and 
friction law. If the load, Fb, is zero, we can find the no-load speed of the rotor. Figure 6- 
13(a) illustrates how the stator surface points’ horizontal velocities will vary from Vh^max 
at the apex, x = 0, to Vh t min at the edge of the contact region, x = x 0 . The rotor, which 
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(b) 


Figure 6-13: (a) The rotor moves at one velocity while individual stator points have different 
horizontal velocity components, (b) The points of slip and no slip corresponding to the points 
shown in (a) are marked. In the region 0 < x < x r , points on the stator surface have horizontal 
velocity components larger than the rotor velocity. In the region.®,. < x < the stator horizontal 
velocity components are smaller than the rotor velocity, while at the point x r , they move at the 
same speed. 

moves at one velocity, will travel at a velocity equal to some stator surface point’s horizontal 
velocity, v ho . 

That is, for Fg = 0 and a given normal force, iyv, the no-load speed of the rotor is: 

_^ho 
^no—load — 

r 

For all other stator surface points, there will be a relative velocity between the rotor and the 
stator, which is a slip condition. These slip regions are illustrated in Figure 6-13(b) where 
for the stator surface point moving at velocity c/ lo , there is no slip. In region 1, all stator 
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surface points have horizontal velocities larger than the rotor velocity and in regions 2 and 
3, all stator surface points have horizontal velocities smaller than the rotor velocity. When 
a load Fg is applied, the friction force acts to oppose it, bringing the rotor to a steady-state 
velocity. 

The forces acting on the rotor and stator are shown in Figure 6-14. In the free-body 
diagram of the rotor, the dotted line represents the deformation the rotor undergoes when 
contacting the stator under an applied normal force, Fjv- This contact condition creates 
a pressure distribution at the interface. We consider only the normal component of the 
pressure distribution, p(x), here, as the curvature of the stator turns out to be small and 
the stator is nearly flat over the region of contact. 

The normal force is related to the pressure distribution by the integral of the pressure 
distribution over the contact area: 

X 0 

F N = Jp{x)ix 

— X 0 

For Coulomb friction, the resulting friction force is: 

f Av(x) 

F t = ~ I a ( ,F dP(x)dx 
J | An (A)| 

X o 

where the friction force is directed so as to resist motion. The relative velocity, Av(x), is: 

Av (A) — V stator (®) 'O rot or 

where v s t a tor (x) signifies the horizontal component of the stator velocity. 

For instance, in region 1, where the velocity of the stator is greater than the velocity 
of the rotor, Av(x) is pointed in the negative n-direction. The friction force acting on 
the stator in that region, Fx : , is positively n-directed to resist the motion. The force 
acting on the rotor is equal and opposite to Fr : of the stator and contributes a force 
component propelling the rotor in the negative ^-direction. Exactly the opposite effect 
occurs in regions 2 and 3 where the stator velocity is less than the rotor velocity, contributing 
force components Fj 2 and Ft 3 of the rotor, resisting rotor motion in the negative ^-direction. 
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Figure 6-14: These free-body diagrams of the rotor and stator illustrate distributed contact and the 
ensuing frictional forces which act in the direction opposite to the direction of the relative velocity of 
the two sliding surfaces. For this region of contact where the rotor and all stator horizontal velocity 
components are moving to the left, in region 1 the stator horizontal velocity components are larger 
than the rotor velocity, and in regions 2 and 3 they are smaller. Consequently, the frictional forces 
in each region are directed as shown. 

The resulting force acting on the rotor is: 

F b = ~F Tnet = - (F Tl - F T2 - F Ts ) 

which is zero, the no-load condition, if: 

Ft i = - (Ft 2 + F Ts ) 

When Fb is increased from zero, the point of no-slip must move down the stator profile, 
as region 1 must grow and regions 2 and 3 must shrink for Fx net to match Fb in steady- 
state. The rotor then moves at the new no-slip velocity which is lower than Vh Q . The point 
of no-slip is designated as x r in Figure 6-13(b) and moves down the velocity profile as Fb 
is increased. This drop-off in speed with increased torque continues until Fb is increased 
to where the no-slip point is at the edge of the contact region, propelling the rotor with a 
velocity of: 

_ v h ,min 
^stator .min — 

r 
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Figure 6-15: For area contact and Coulomb friction, the speed-torque curve falls off with torque 
from a given no-load speed up until the no-slip point reaches the edge of the contact region, which 
is moving at u>stator,min- Larger loads, Fb, cause slipping over the entire area of contact and the 
rotor reaches a limiting maximum torque. The exact shape of the curve in the region marked by the 
dotted line is yet: to be determined and will depend upon the pressure distribution and the model 
of the contact mechanics at the interface*,, 


For any larger loads, all points of contact will be sliding and the rotor will deliver its 
maximum torque: 


1~maxi — [Id I X t 

Maximum torque occurs at the point where x r = x 0 . That is, when the braking force, 
Fg , is increased just to the point where the rotor begins to slip everywhere, the motor is 
delivering its maximum net frictional force: 


Fb 


max 


— Ft 

-L net,max 


~(F t 1 - Ft 2 - F Ts ), 


Fx net is maximized when Tp, and Ft 3 are zero, which is the case when x r = x 0 . 

Figure 6-15 depicts the resulting speed-torque curve. The velocity of the stator at its 
apex, oj s t a tor,max, i s fixed by stator geometry and excitation voltage. The stator velocity at 
the edge of the contact region, oj s t a t 0 r,min, is determined by the interpenetration depth due 
to iqv- For an interpenetration depth less than w 0 , the no-load speed will be some point 
between testator,max and u>stator,min'-, ■ As Fb is increased above zero, the speed of the 

motor falls off from to ^stator,min■ We show a dotted line in the figure as at this stage 
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^rotor 



7 "m<i:r j ^~max2 ^~max 3 


Figure 6-16: Increasing the normal force, Fjv, creates a larger area of contact and therefore a lower 
u) s t at or, min- The speed-torque curve then shifts down as shown. 

the shape of this curve is unknown, and will depend on the deformation at the interface 
and the pressure distribution. 

Once the no-slip point lias moved to the edge of the contact region, larger torques cannot 
be sustained, and the motor reaches its torque limit. Loads larger than this cannot be met 
with opposing friction forces and so the rotor decelerates. Equilibrium points are shown as 
the vertical solid line, where increasingly larger loads, Fb, cause the motor to slow down, 
but the motor cannot push back with any force larger than r maxi . Finally, the speed of the 
rotor becomes zero and eventually reverses direction. 

Note that the point where the speed of the rotor is zero is defined as the stall torque. 
For contact conditions where 0 < x 0 < the maximum torque is the same as the stall 
torque and occurs when x r = xo, as shown by the vertical line in Figure 6-15. 

The speed-torque curve of Figure 6-15 depicts a set of operating points for one given 
normal force and excitation voltage. Figure 6-16 shows liow the speed-torque curve shifts as 
iqv is increased from the initial value of to iqv 3 , assuming iqv 3 causes an interpenetration 
depth just equal to w 0 . An increase in normal force from F^ 1 to iqv 2 causes the edge of the 
contact region to move down the stator displacement profile, resulting in a lower testator,min- 
The maximum stator speed, testator,max, is not affected by the normal force in this model 
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which assumes a rigid stator, and so for Fb = 0 and Fn = Fjv 2 , the no-slip point will start 
at a lower position on the the stator profile than for Fn = -Fjv t in order to maintain the 
no-load equilibrium condition: 

F Tl = - ( F T2 + F Ts ) 

As Fg is increased from zero, the no-slip point, x r , continues to move lower and the 
rotor velocity falls off with rotor torque until it reaches ustator,min 2 ? a t the new maximum 
torque: 

T~ r max 2 = HdFN 2 r 

When the normal force is increased to F^ 3 where the interpenetration depth is w 0 , then 
x 0 = j and the edges of the contact region are at x = —^ and x = At these points, the 
horizontal component of the stator velocity is zero and therefore u s tator,min 3 is zero. The 
no-load speed will be somewhere between testator,max and ^stator,min 3 and the velocity of the 
rotor will drop off with the load torque until the maximum torque is reached at: 

^~max 3 = HdFN 3 r 

For larger loads, Fb, the rotor begins to decelerate and spin in the other direction as the 
rotor cannot supply enough torque to match the load. 

The scenario just related describes quadrants I and IV behavior. Loads Fb applied in 
the negative ^-direction yield quadrant II behavior where rotor velocity increases from the 
no-load speed to testator,max U P until a limiting negative torque is reached at —HdF^r where 
i = 1 to 3. Again, we draw dotted lines in the fall-off region as the exact form of these 
curves has not yet been determined. 

While Figure 6-16 illustrates how a family of speed-torque curves evolves, for quantifying 
a figure of merit for classifying the performance limits of this class of actuator technology, 
we are especially interested in how the stall torque for each curve changes as a function of 
F n . 

Figure 6-17(a) illustrates the condition of r maxi for a normal force of F^. A pressure 
distribution, p(x), is generated over the area of contact. At the edge of the contact region, 
p(x) = 0. When a braking load of magnitude Fs maxl is applied such that the motor is 


83 



supplying its maximum net torque, then: 


-Ft. 


net .max 


= -Ft , 


and all points on the surface of the stator in the region of contact have horizontal velocity 
components greater than the rotor velocity, and the no-slip point has reached the edge of 
the contact region; x r = x 0 . 

Figure 6-17(b) illustrates the maximum torque condition when the normal force has 
been increased to Fjv 3 , where x 0 = The pressure distribution now extends from x 0 = —^ 
to x 0 = ^ and has a magnitude such that: 


F ]v = p (x ) dx 


For a pressure distribution which is proportional to displacement, most of the total 
pressure derives from the region near the apex. When a braking load of Fs max?j is applied, 
ff ]ietmai3 must again be composed solely of contributions from Ft i: and so x r = x 0 = 

As FT net is related to the normal force by: 


F T net = PdF N = Pd J p(x)dx 


we realize: 


Ft, 


> Ft 

J-n 


net,max 3 -‘-net,maxi 


since: 

Fn 3 > F n i 

Figure 6-18 depicts the geometry of the rotor and stator when the normal force is large 
enough that the interpenetration of the rotor is greater than w 0 . Figure 6-19 illustrates 
the resulting pressure distribution and maximum torque condition for this case of Fjv 4 
where the contact region has moved beyond x 0 = ^ and encompasses points on the velocity 
profile where £ < 0 and stator horizontal velocity components are pointed in the positive 
^-direction. 
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Figure 6-17: (a) Normal force creates a pressure distribution over the region of contact. As 
the braking load is made larger and larger, Fx nei increases to match it as the rotor slows down and 
the no-slip point,. x r , moves down the velocity profile to x r = x 0 . (b) For a larger normal force Fn 3 , 
which squishes the rotor and causes it to comply such that x 0 = a larger braking load can be 
applied before Fx nei is reached, since Fx nei is directly proportional to F^. In this case, the onset 
of this max condition occurs when x r = x 0 = At x r = u>rotor = 0. 
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Figure 6-18: If the rotor deforms more than the amplitude of stator displacement, w e , horizontal 
velocity components below £ = 0 will work against propulsion of the rotor. 



Figure 6-19: As the normal force is increased such that the area of contact extends beyond 
x 0 = j, contributions to Fy net are picked up which work against the driving force. However, these 
contributions are small, as the pressure is low at the edges of the contact region. The increase in 
the magnitude of p(x) near the apex more than offsets these contributions and can lead to a larger 
stall torque, T sta u 4 , than 

^stalls • 
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In these regions, the stator has horizontal velocity components opposing the direction of 
the velocity of the rotor. However, due to the larger normal force, F^ i , a larger load, i^B max4 , 
can be applied than FB max3 . This condition can occur at a negative velocity where x r = x 0 
is beyond As shown in Figure 6-20, the onset of r maXi occurs at the negative velocity 
^ stator, min 4 , the velocity corresponding to the no-slip point, x r , being at the edge of the 
region of contact. While T maXi occurs at a negative velocity, the stall torque corresponding 
to F Ni is still defined as the torque at which Lv roior = 0. That is, for cases where ^ < x 0 < 
the stall torque occurs at the point where x r = 

Note that in this scenerio, r sta ;; 4 can be larger than T sta u 3 , even though the contact 
area has extended to regions on the displacement profile where stator horizontal velocity 
components oppose the direction of the rotor velocity. The reason this is possible is that 
while the normal force, Fjv 4 , is larger (and hence Te max4 is larger) than for the case of 
F N3 , the pressure distribution remains primarily over the apex where stator horizontal 
velocity components are moving in the same direction as the rotor. Since p(x ) has grown 
in magnitude, relatively small contributions to the friction force due to the relation: 

X 0 

F T„et = dd J p(x) dx 

— X 0 

come from the additional p(x) at the edges of the contact region covered by this new pressure 
distribution. 

Would it be possible for an even larger normal force, Fjv 5 , to shift the speed-torque curve 
in such a way that the new stall torque, T sta n 5 , would be less than r sta ;; 4 ? In this model, 
where we assume the stator profile to be undistorted and the pressure distribution to be 
proportional to that profile, larger normal forces will only contribute to a larger Fj net max 
and pick up very little opposition from the edges of the contact region where the pressure 
is low. For our model here, we would expect the stall torque to asymptotically reach a 
maximum as Fjv is increased to the point where there is complete interpenetration of the 
rotor into the stator, that is, where x 0 = This case is also illustrated in Figure 6-20. 

The variations of stall torque and no-load speed with Fjv are shown in Figure 6-21. The 
fall-off trend of no-load speed with normal force for first quadrant operation is plotted in 
Figure 6-21(a) where we see that Lv no _i oa( i decreases as normal force is increased. Comparing 
this graph to the corresponding curve for line contact, Figure 6-6 (a), we see that the no- 
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^rotor 



Figure 6-20: When the normal force is increased from Fn 3 to F^, the maximum torque, T maXi , is 
increased and the onset of this maximum condition occurs at a negative velocity, u> s tator,min 4 - Note 
that Tgtaiu f an be larger than T sta u 3 , even though x 0 > For this model of a non-compliant stator, 
further increases in Fjv never coerce the stall torque to decrease. 


load speed for area contact predicts a drop-off in LJ no -i oa d with iqv as opposed to remaining 
constant as in the case for line contact. 

Similarly, the stall torque variation with iqv is plotted in Figure 6-21(b), where stall 
torque reaches a maximum at x 0 = Comparing to Figure 6-6(b), the T sta ii vs. iqy 
curve for line contact, we see that the area contact condition asymptotically approaches the 
prediction for the line contact case for small applied normal forces, but falls off for large 
F n . 

A t s t a u vs. iqv relationship like that originally sought in Figure 6-l(a), where T sta ii 
falls to zero for very large values of iqy would not be realized here. Extensions to this 
model of area contact with Coulomb friction between a rigid stator and compliant rotor to 
incorporate the condition of a compliant stator would lead towards such expectations, but 
for now we examine the model here more closely to investigate exactly how the performance 
curves develop for various assumptions about the contact mechanics at the interface. 
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no—load2 
^no-load 3 
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^no-load 5 


t t t t h Fn 

Fn } F N2 F/v 3 F n 4 F N5 

(a) 

Fstall 



Figure 6-21: (a) Assuming area contact, Coulomb friction, a rigid stator and compliant rotor, we 
can predict that, the no-load speed will decrease with normal force from a maximum of w no _; oa( j = 
Vh '™ ax . (b) Stall torque will increase with Fn, asymptotically reaching a maximum at i 0 = 
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Figure 6-22: The coordinate system for the model of a cylinder contacting an elastic half space is 
centered on a peak at the original top surface of the undeformed stator and the area of contact is 
taken as a. The edge of the contact region then, is x 0 = In the region 0 < x < x r , the horizontal 
velocity components of the stator surface points are larger than the rotor velocity. In the region 
x r < x < x 0 , the stator horizontal velocity components are moving slower than the rotor. At x r , 
the rotor and stator move at the same speed. 

6.3.2 Hertzian Contact Model 

Arguing physically, we have found bounds for motor torques and identified contact mech¬ 
anisms for describing liow the speed-torque curves shift with normal force. However, the 
exact form of the speed-torque drop-off, marked by the dotted lines in the previous figures, 
needs to be determined by examining the pressure distribution and contact area at the 
interface. 

In this section, we investigate the effects of using a Herztian model of contact [Hertz 
1882] to predict the speed-torque characteristics. Figure 6-22 illustrates the coordinate 
system used here along with the definition of a, the width of the area of the contact region 
used in Hertzian theory. 

We note that the friction force is related to the rotor speed by: 


F t = -/.id 


/ 


— X 0 


Av (x) 

I Ac (x )| 


p ( x) dx 


which is, equivalently: 


•L o 

Ft = —'2pd J s 9 n [Ac (x )]p (x ) dx 
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Figure 6-23: One model for a contact pressure distribution is a cylinder contacting an elastic 
half-space, Here, the stator is modeled as a cylinder with an equivalent radius of curvature and i lie 
rotor is modeled as an elastic half-space contacting the cylinder over an area of width a. 

where the sgn [Ai'(r)] function is defined as: 

1 Vstator (- r ) Vrotor 

sgn [An (x )] = 0 v stator (x ) = ^ rotor 

- 1 ^stator (* e ) ^ ^rotor 

arid: 

V stator (x) — l 'j, ,,,, a.r f ' f> S (kx Ojt) 

where for this coordinate system, the figures are again drawn for ujt = 0. 

Finding the area of contact as parameterized by x 0 will lead us to the limits of inte¬ 
gration, and our model of the pressure distribution will lead us to p(x). This friction law 
will then give us a set of operating points along the speed-torque curve for rotor speeds of 
v h,min < v rotor < v h,max- That is, in this section we will be calculating the speed-torque 
curves valid for the regions drawn by dotted lines in Figure 6-16 and Figure 6-20. 

While there are innumerable models of contact mechanics, and [Johnson 85] gives com¬ 
prehensive coverage in this area, one simple model of contact applicable to the problem 
here, is that of Hertzian contact of a cylinder contacting an elastic half-space as shown 
in Figure 6-23, where we model the deformed stator as a cylinder of equivalent radius of 
curvature. The rotor is modeled as an elastic half-space contacting the cylinder over an 
area of width, a, when a force per unit length, iqv, is applied. [Young 89] gives the width 
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of the area of contact, a, for this model as: 


a = 1.60\/ FnDC'e 


where Fn is the applied normal force per unit length, D is the diameter of the cylinder and 
C'e takes into account materials parameters: 


C'e 


1 ~ v\ , 1 ~ v l 

E x + E 2 


Here, E\ and E 2 are the moduli of elasticity for the two materials at the interface, and 
v\ and v 2 are their corresponding Poisson’s ratios. 

This cylinder model will be adequate so long as the undulating stator can be modeled 
by the curvature of a cylinder over the region where the stator contacts the rotor. That 
is, for values of Tjv where the rotor interpenetrates only as far as the top half of the 
stator’s sinusoidally-shaped surface, a cylinder on an elastic half-space model should suffice. 
Certainly for Tjv such that the interpenetration depth is larger than w 0 , where the stator 
curvature has an inflection point, we would not expect this model, where a is assumed to 
grow as the y/F^r, to apply. 

The constants D and C'e are determined by w 0 , the stator deflection (which in turn 
is a function of the applied drive voltage), and the material properties of the frictional 
interface. For example, for the case of an aluminum rotor contacting a steel stator, if we 
take the modulus of elasticity for aluminum as E a i = 70 GPa, the modulus of rigidity of 
aluminum as G a i = 26 GPa, the modulus of elasticity for steel as E st i = 200 GPa and the 
modulus of rigidity of steel as G s ti = 77 GPa, Poisson’s ratio can be determined from: 


G = — -- 

2(1 + //) 

resulting in v a i = 0.35, v st i = 0.38 and Ce = 1-68 X 10 -11 ^. 

The diameter of the equivalent cylinder, D, is twice the radius of curvature of the 
deformed stator, 2 p K , which is related to the stator deflection by: 


1 

Pn 


d 2 t{xG) 

dx 2 


k 2 w 0 


2ir 

T 


2 


W 


o 
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For our prototype motor of Figure 3-5, A = 5.1 mm and w 0 = 0.1 p, giving D = 13.2 m. 
The contact area for a unit force per length of Fn = 1 ^ works out to be a = 24 pm. That 
is, the percentage of contact area of a half-wavelength of the motor is approximately 1% for 
unit normal force per length. In terms of the limits of integration of our friction law, this 
contact area width, a, is related to x 0 by: 


a 

X o —— 

2 

After calculating the width of the contact area, the pressure distribution remains to 
be found. The pressure distribution at the interface can also be modeled to various levels 
of sophistication [Johnson 85]. Here we assume p(x) is proportional to the transverse 
displacement of the stator, £ (x,t) = w 0 cos (kx — ut ), where, again, the figure of the stator 
has now been drawn for ut = 0. Consequently we take: 

p (x) = A [cos (kx) + C] 

where A is a function of Fn (and therefore x 0 ) with units of ^ which satisfies: 

X 0 

f n = / „(,)<(, 

— Xo 

Here, C is a constant which accounts for the boundary condition that p(x) must equal 
zero at x = — x 0 and x = x 0 . We find C to be: 


C = —cos (kx 0 ) 

As the pressure distribution is symmetric about x = 0, the expression for the normal force 
simplifies to: 

Xo 

Fjv = 2 A j [cos (kx) — cos (kx 0 )\ dx 
o 

yielding an amplitude of the pressure distribution of: 


A = 


F N k 


2 [sin (kx 0 ) — kx 0 cos (kx 0 )\ 
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Now that we have determined p(x), we can substitute it into the equation for the friction 
force, Ft- The limit of integration, x 0 , is calculated from the width of the contact region, 
a, and is a function of F Xr: 

x 0 = - = 0 . 80 \JFnDCe 

The velocity of the stator is: 

Vstator (®) — ^h,maxCOS (kx~) 


giving the friction force acting on the stator as: 


fT„,., = -2 W A/, 9 M^co,(fa)-w][co,(U)-co,(U 0 )]^ 

0 

The net friction force per unit length that the rotor applies to the load is equal and 
opposite to the friction forces per unit length acting on the stator, as was shown in Figure 6- 
14 : 


F T„ e t = -Ft,, 


F T net , then, is: 


F Tnet = 2g d A j sgn [v h , max cos (kx) - v rotor } [cos (kx) - cos (kx 0 )\ dx 
o 

Breaking the sgn[ \ function into two integrals, where the stator velocity is greater than the 
rotor velocity for 0 < x < x r and less than the rotor velocity for x r < x < x 0 , the net force 
per unit length that the rotor applies to the load becomes: 

X r X 0 

F T ,„ = 2w .4/[«»(fa) - coHkX')]ix-2 M Aj[coMkx) - coHkx,)]ix 

0 Xr 

which, after integration, simplifies to: 


FT ne t 


VdFN 


[sin (kx 0 ) — kx 0 cos (kx 0 )\ 


[2 sin (kx r ) — 2kx r cos (kx 0 ) — sin (kx 0 ) + kx 0 cos (kx 0 )\ 
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Plotting T rotor vs. tjJ ro tor gives us the desired speed-torque relationship that we were 
looking for, as outlined earlier in Figure 6-16. We can use the values taken before, ^ = 0.3 
for the coefficient of dynamic friction, k = 1.23 X 10 3 m -1 for the wavenumber, D = 13.2 m 
for the radius of the cylinder modeling the stator, r = 3.25 mm for the radius of the motor, 
A = 5.1 mm for the spatial wavelength of the motor, C'e = 1-68 X 10 -11 ^ for the materials 
interface parameter, and plot the speed-torque curves for various values of Fjv- As Fjv is 
increased, the contact area becomes wider and x 0 therefore becomes larger. Figure 6-24 
illustrates speed-torque curves for Fjv values yielding deformation where x 0 moves down 
the sinusoidal stator profile by ^ where n ranges from to 1 to 4. At x 0 = the rotor and 


95 




Figure 6-24: These are the speed-torque curves for Hertzian contact and Coulomb friction for 
different values of Fjv such that the edge of the contact region, x 0 , moves down the stator profile 
from xq = yjr to xq = y|. 


stator are contacting only along a line at the peak of the wave. At x 0 = the rotor lias 
deformed over the stator from x = — y to = y. That is, the rotor lias deformed over the 
top half of the cylinder which is depicted in Figure 6-23. 

These curves represent the dotted-line regions predicted in Figure 6-16 for quadrants I 
and II behavior. The curves here are graphed out to the point of maximum positive torque 
in quadrant I, where they then fall off vertically at positive stall, and out to maximum 
negative torque in quadrant II, where they increase vertically. Note, that as predicted in 
Figure 6-16, for large negative torques, the rotor speeds approach testator,max- Also, we see 
that for x 0 = y, the onset of maximum torque, T maX3 , occurs at u> ro t or = 0 as predicted in 
Figure 6-16. 

We can try to continue utilizing this model of a cylinder contacting an elastic half-space 
beyond the condition that x 0 = y. That is, if we continue with the model of a growing as 
\/Fn and: 

1 — 1 ( ^rotor \ 

x r = —COS - 

™ stator,max J 

then we will find speed-torque curves such as those illustrated in Figure 6-25. There we 
see that as x 0 is increased to x 0 = and so on up to x 0 = the stall torque continues to 
increase. 

After characterizing the speed-torque curves, we would like to find ^ no -i oa d vs. iqy 
and t s t a u vs. iqv- The no-load speed, u> no -i oa d, can be found by setting r rotor = 0 in the 
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Figure 6-25: For Hertzian contact of a cylinder and an elastic half-space, as x 0 is increased beyond 
x 0 = the speed-torque Curves shift such that the no-load speeds decrease while the corresponding 
stall torques increase. Here, x 0 moves down the stator profile from *0 = to xo = 4- 

speed-torque equation: 

° = ^ ( ' Xr ^ ^ )] 

which reduces to: 

<j)(x r ) = ^(x 0 ) 

We solve for x r and convert to LJ no -i oa d with: 

^no-load. — ^statoi^max^^'-i^^r) 

Since each value of xo is related to the applied normal force per unit length from: 

Fk = W^(oJo ) 2 

the 0 J no -i oa( i vs. Fn performance curve for the motor as outlined in Figure 6-21(a) can be 
found. ijJ no -i oa< i is plotted vs. iqy in Figure 6-26(a) for the same constants used earlier and 
for values of x 0 ranging from ^ to We see that uJ no -i oa d falls off with Fn as expected. 
Turning to stall torque, we find that the maximum torque per unit length is developed 
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(b) 


Figure 6-26: (a) u> n o-ioad is plotted fprsus Fn for values of Fn creating regions of contact param¬ 
eterized by x 0 ranging from ^ to (b) Th#?:r sta // vs. Fjv curve is plotted for the same values of 
Fn- Tstaii is seen to be linear with Fjv in this model. 
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when x r = x 0 in the speed-torque equation, leaving: 

T m(l x = VdF N r = T st all f OX 0 < X 0 < j 


Again, this is consistent with our physical interpretation in Figure 6-16, 
directly proportional to Fn up to For contact beyond x 0 = stall 

when x r = 


Tstall 


fJ-dFNr 

<t>{x 0 ) 


2 <f> 


A 

4 


4>{Xo) 


for 


A 

4 


A 

< x 0 < - 


where T sta ii is 
torque occurs 


Simplifying, T sta ii reduces to: 


Tstaii 


fJ-dFNr 

4>{x 0 ) 


[2 - 4>{x 0 )\ 


t A A 

for - < x 0 < - 


where for Hertzian contact, Fjy and x 0 are related by: 


x 0 = 0.80 V FnDC'e 


This T sta ii vs - Tjv curve is plotted in Figure 6-26(b) for normal forces all the way up to 
deformations of x 0 = We see that for A < Xo < T s t a ii increases at a lower rate than 
for 0 < x 0 < A. 

Note that T s t a lh T rotor an( i have been given throughout as torques and forces per 
unit length. For total torque or force magnitudes, these values should be multiplied by b, 
the stator width. For the motors of Figure 3-1, b is 3 mm. 

The mechanical output power can be calculated as the product of rotor torque and rotor 
speed in quadrant I operation (i.e. where rotor torque and rotor speed are both positive: 

Fmech — ^rotor Aotor 


Graphs of output power in this region of operation are shown in Figure 6-27. Peak 
output power is shown to occur in the graph for x 0 = with a value of approximately 
31 — and a torque per unit width of 15 j? or the motors of Figure 3-1 where the width, 

m l Jr m o 7 

b, of the stator is 3 mm, and the number of wavelengths in the motor is 4, we find maximum 
total output power to be 370 mW at a total output torque of 0.18 Nm. In reality, our model 
undoubtedly breaks down in regimes of such large normal forces as it does not take into 
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account damping and compliance in the stator. 

Perhaps more useful figures to keep in mind would be power and torque at a normal 
force corresponding to x 0 = There we find that maximum output power per width for 
one wavelength is approximately 11 ^ at a torque per unit width of 2.5 and rotor 
speed of 4.4 (i.e. 42 rpm). For a 3 mm width and 4 wavelengths, total output power is 

130 mW and total output torque is 0.030 Nm (i.e. 3.0 kgf-cm or 4.2 oz.-in.). The normal 
force here is approximately 8.6 N (i.e. 0.87 kgf or 1.9 lbs.). 

The efficiency for this Hertzian model of area contact can also be calculated. The 
efficiency is the ratio of mechanical power output to power input: 


P 


Vm — 


mech 

~P 

J- m. 


P 


mech 


Pi 


loss 


P 


mech 


where Pi oss is the frictional heating loss which arises from the product of the rotor torque 
and the relative velocity of the two surfaces. We formulate Pi oss using the same integral as 
was taken earlier to find the rotor torque, but now multiplying the pressure distribution by 
the relative velocity between the rotor and stator in the integrand: 

X 0 

Pioss = 2/i d r J sgn[Au(x)]p(x)Au(x)dx 
o 


Here, A^(r) is the relative velocity between the rotor and stator, testator ~ ^ rotor , and 
sgn [Au; (A)] is the direction of the relative velocity of the rotor and stator as before: 


sgn [Au; (x)] = < 


1 ^stator ^ ^rotor 

0 ^stator (^) — ^rotor 

^ 1 ^stator ^ ^rotor 


Splitting the integral into two parts over the separate regions where the relative velocity 
changes sign, Pi 0S s expands to: 


x r 

Pioss = 2 Apdr f [cos ( kx ) — cos ( kx 0 )\ [w sm cos ( kx ) — u r \ dx 
0 


Xrj 

— 2 Apdr f [cos (kx) — cos ( kx 0 )\ [oJ S mCos (kx) — ui r ] dx 

X r 
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(a) 



/ Nm\ 

\ m J 
(b) 

Figure 6-27: (a) Mechanical power output per width, 6, of the motor over one wavelength is shown 
vs. T ro t or for values of x 0 from yj? to y for quadrant I operation where T ro tor and u) rotor are both 
positive and the motor is doing work on the external system, (b) The corresponding curves for x 0 
from y to 4 show that mechanical output power reaches a maximum for a normal force which causes 
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Vm (1 = 100 %) 



Figure 6-28: (a) Mechanical efficiencies for quadrant I behavior are shown for speed-torque curves 

where the normal force is such that x 0 varies from yjr to (b) Mechanical efficiency decreases with 

increasing F^. Graphs here correspond to interpenetration of the rotor into the stator where x 0 

varies from 77 to 4 . 

16 2 
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where we have used uj sm for u stator,max and w r for uj ro t or - After integration and algebraic 
manipulation, Pi oss can be expressed in the form: 


Ploss = ^[2$(j: r )-$( j£ 
<p yXo) 


where we have substituted for: 


(x) = —kxu) srn + —uj sm sin (2 kx) — [oj sm cos (kx 0 ) + ui r \ sin (kx) — kxuj r cos ( kx c 


and (f>(x) is, as before: 


cf) (x) = sin (kx) — kxcos ( kx 0 ) 


Plugging Pio SS and P mec h into °ur expression for mechanical efficiency: 


P 


Vm — 


mech 

~P 

J- in. 


P 


mech 


Pi 


loss 


P 


mech 


rj m becomes: 

_ 7~rotor^rotor 

[2$ (X r ) - $ (X 0 )] + T rotor U rotor 

The mechanical efficiencies for quadrant I operation of the speed-torque curves of Fig¬ 
ure 6-24 and Figure 6-25 are plotted in Figure 6-28. We see from Figure 6-28(a) that while 
maximum efficiencies decrease with increasing Fjv, for very light normal forces (i.e. Fjv such 
that x 0 = the maximum efficiency approaches that for the ideal case of line contact, as 
was illustrated in Figure 6-7(b). 

6.3.3 Linear Spring Model 

In the Hertzian contact model, we used an expression for the contact area that Hertz 
derived based on the stress distribution in two contacting cylinders, and we assumed a 
pressure distribution that was proportional to the displacement of the rigid stator. Taking 
one cylinder to have an infinite radius produced a contact between a cylinder and an elastic 
half-space. The relationship between Fjv and x 0 in the Hertzian formula effectively provided 
a stiffness of the system. 

We can see this via a simpler model which considers the rotor to be merely a distributed 
linear spring. While perhaps less accurate, this viewpoint can be helpful in building in- 
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tuition. Figure 6-29(a) depicts this model where the rotor has a thickness of h r and a 
distributed stiffness of K r . The rotor is modeled as a series of springs placed in parallel 
along the rotor. That is, the displacement of the rotor at any position, x, is unaffected by 
the piece next to it at x + Ax. If the area of contact for a given normal force per unit 
length, Fjv, extends from x = — x 0 to x = x 0 , we take the normal force per unit length to 
be related to the pressure distribution by: 



p(x ) dx 


— Xrj 


and: 


p (x) = K r [vj 0 cos (kx) — w 0 cos (kx 0 )\ 


for 


0 < x 0 < 


A 

2 


as w 0 cos (kx) — w 0 cos (kx 0 ) is the displacement of each spring at any location x in the region 
of contact. K r here has units of since Tjv is taken as force per length. For a rotor of 
thickness h r , K r is simply the modulus of elasticity of the rotor divided by its thickness: 


K r 


E r 

h r 


Carrying out the integral for the normal force, we find that Fn is related to x 0 by: 


F n 


2 K r w 0 
k 


[sin (kx 0 ) — kx 0 cos (kx 0 )\ 


Figure 6-29(a) illustrates a contact case where a normal force, F^, creates a region of 
contact where 0 < x 0 < For some normal force, Fjv 2 , there will be a contact condition 
where x 0 just equals ^ as in Figure 6-29(b). For the scenerio where x 0 = and Fn is 
increased further, the pressure distribution becomes dependent upon the degree of com¬ 
pression of the rotor. For this case, we take a coordinate system as shown in Figure 6-29(b) 
where the ,?-axis is pointing down and the origin is centered at the top surface of the rotor 
just as the bottom of the rotor is touching the valley of the undulating stator. As the 
normal force is increased from Fn 2 to Fn 3 as shown in Figure 6-29(c), the displacement of 
the spring at x 0 = then becomes h r — d and the resulting pressure distribution is: 

p(x) = I( r [vj 0 cos (kx) — (h r — d)] for x 0 = — 
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Figure 6-29: (a) This is a model of the rotor as a linear elastic spring with a distributed stiffness 
of A" r . (b) At complete interpenetration, for a given normal force, Ajv 2 , x 0 equals 4- ( c ) If Fn is 
increased beyond Ajv 2 , the rotor becomes compressed at all locations x. 

Going back to the case where 0 < x 0 < we see that the pressure distribution for a 
linear spring: 


P ( X ) LinSpr = K r [ W 0 COS (kx ) - W 0 COS (kx 0 )] f 01' K r 


E r 

h r 


is similar in form to the pressure distribution for Hertzian contact: 


P ( X )Hertz = ^ [ C0S ( kx ) “ COS ( kx 0 )] f 01' A = 


Fn k 

2(f) (x 0 


We can find a rotor thickness, that signifies the equivalence between these two models 
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by setting: 


P(%)LinSpr — P(.3'') Hertz 


The cosine terms drop out, leaving: 

A = K r w 0 

Substituting in the expressions for A and K r , h r is found to be: 


h r 


2 E r w 0 <f>(x 0 ) 
kF N 


This can be simplified further by noting that the Hertzian formulation gave a relationship 
between Tjv and x 0 of: 

where D was related to w 0 by the radius of curvature, p K \ 


D = 2 p n and — = k 2 w 0 

Pn 


We hnd the ratio of w 0 to Tjv to be independent of D: 

w 0 _ 2 C'e (0.80) 2 
Fn k 2 x 2 0 

Plugging y 2 - into the expression for h r , h r can be seen to be a function of x 0 : 


h r = 2.56 


E r C E <t> (% c 

1,3 ™ 2 

lb dj _ 


Thus at a given x 0 , the Hertzian formulation is equivalent to a linear spring model with a 
distributed stiffness of: 

K r = —— = 0.39- 


E r k 2 x 2 


h r C E 4* (^o) 

The speed-torque equation for this linear spring model is found by solving the friction 
integral: 

Ft „, = 2„ d K r w/j [cos (*x) - cos (*x„)] c/x - 2 W A >*/ [cos (*x) - cos { k x „)\d x 
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Figure 6-30: These are the speed-torque curves for a linear elastic spring model of the rotor with 
Coulomb friction for various values of x 0 . Here, x 0 moves down the stator profile from xq = A to 

= f- 


After carrying out the integral, Fx net is: 


Fx nPt = ‘ ^ d ' —- X [2sin (kx r ) — 2kx r cos (kx 0 ) — sin ( kx 0 ) + kx 0 cos ( kx 0 )] 


If we again take cf>(x) such that: 


(f>(x) = sin (kx) — kxcos (kx 0 ) 


the expression for Fx net reduces to : 


F T„et = 2fldI ^ rWu [2(/>{x r ) - (t>{x 0 )] 

As the speed-torque equation is simply T rotor = Fj net r , we find: 

•2n d K r w 0 r 

Trotor = -£- [2(f) (x r ) ~ <f) (x 0 )J 

We are still assuming that the stator is rigid, and so the velocity profile remains a cosine 
function and x r is, as before: 


x r 


—cos 
k 


_-i 


^rotor 


stator, r 


The resulting speed-torque curves are plotted in Figure 6-30 and Figure 6-31. We have 
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Figure 6-31: These are the speed-torque curves for larger values of Fjv such that x 0 moves down 
the stator profile from *0 = y| to *0 = 4- 

used the same values for the constants, w 0 and r, as in the speed-torque plots for the 
Hertzian model of Figure 6-24 and Figure 6-25. For K r = we have taken the rotor as 
aluminum, where E r = E a j = 70 GPa as before and we have chosen so as to make K r be 
the equivalent stiffness of the Hertzian case for x 0 = That is, for the k and Ce values 
used earlier, works out to be 0.31 mm and K r to be 2.3 X 10 14 

In comparing the speed-torque curves of the linear spring model to the speed-torque 
curves of the Hertzian model, we see that while we have set the curves to be the same for 
x 0 = the linear spring model predicts the same no-load speeds for the family of curves, 
but predicts larger stall torques for a given interpenetration, x 0 . 

The equivalence of the no-load speed as a function of x 0 between the two models can 
be seen analytically. The no-load speed is calculated by setting T rotor = 0, giving the same 
relationship for the Hertzian case: 

(f){x r ) = ^4>{x 0 ) 

However, the relationship of x 0 to iqv is different here than in the Hertzian case and we 
see from comparing Figure 6-32(a) and Figure 6-26(a) where LJ no -i oa d is plotted vs. iqv for 
each model, that larger normal forces are required here for each subsequent value of x 0 after 

x - A. 

•Co — 16 - 

The stall torque is also similar to the Hertzian case. For a normal force where 0 < x 0 < 
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Figure 6-33: (a) These curves represent the mechanical output power of the motor as a function 
of torque for normal forces such that x 0 varies from x 0 = yjr to x 0 = (b) Here, output power is 

graphed for interpenetration of x 0 from x 0 = y| to x 0 = 4- 
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we again have that r sia u occurs when x r = x 0 . 


Tstall 


2/i, d K r w 0 r 

k 


4>{x 0 ) 


0 < x 0 



In this linear spring model where x 0 and F Xr are related by: 


F n 


2 K r w 0 

k 


4>{x 0 ) 


the stall torque reduces to: 


Tstall = VdF N r 0 < x 0 < - 


where T sta ii is again linearly related to F Xr. 

We have a similar relationship for T sia u in the regime of ^ < x 0 < where r sia u occurs 
when x r = x 0 : 


Tstall 


2 fj, d K r w 0 r 

k 


[2 -4>{x 0 )\ 


A 

4 


A 

< x 0 < - 


Substituting the relationship between x 0 and _Fjv gives the same expression for stall torque 
as in the Hertzian case: 

_ n d F N r 
Tstall ~ <j>(x 0 ) [ 

Comparing Figure 6-32(b) with Figure 6-26(b), we see that the main difference between 
the two models is the stiffness of each system, where for the linear spring model, which is 
stiffer, it takes more normal force to achieve the same x 0 . 

Mechanical output power for this linear spring model is again taken as: 


Pmech — Fotor^rotor 


Figure 6-33 graphs this function vs. T rotor for every operating point along each speed-torque 
curve. Comparing to Figure 6-27, we see that for a given x 0 , power output is uniformly 
larger than for Hertzian contact as the ability to apply larger normal forces for a given x 0 
leads to larger torque outputs and hence larger power outputs. 

The expressions for Pi oss and r] m are also the same as for the Hertzian case and graphs 
of efficiency for the linear spring model are shown in Figure 6-34. 
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Figure 6-34: (a) Efficiency vs. torque curves for x 0 ranging from ^ to are similar in shape to 
those of Hertzian contact, except that the corresponding stall torques are larger, (b) Efficiency vs. 
torque curves here range from x 0 = y| to 4- 
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6.4 Extensions 


The models described to this point are useful in building intuition about which assumptions 
and constraints impose which kinds of performance behavior. The closed-form solutions de¬ 
veloped to this point do not take into account squashing of the stator deformation. The 
normal force pushing the rotor against the stator will act against the piezoelectric actuation 
of the stator and create a reduced amplitude of deformation. Approximate solution tech¬ 
niques, such as Rayleigh-Ritz methods, can be used to formulate a model which combines 
the various models of contact described in this section with stator squashing. [Hagood and 
McFarland 95] have recently implemented such a model using a linear spring contact mech¬ 
anism. Appendix B gives background on these variational techniques and shows simulation 
results for the 8 mm bulk motors discussed here. 
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Chapter 7 


Experiments with 8 mm Bulk 
Motors 


7.1 Initial Tests With Early Prototypes 


Measuring speeds and torques of microfabricated ultrasonic micromotors such as the one 
shown in Figure 2-7, is rather difficult as the torques are quite small, bearings are not 
incorporated into the fragile membrane structures and the fabrication iterations require 
significant turn-around time. 

Consequently, we have fabricated a number of ultrasonic motors from bulk PZT and 
steel parts cut in a machine shop at a scale which is easier to handle than the micromotors 
of Figure 2-7 yet still smaller than any commercially available electromagnetic motors. One 



Figure 7-1: This is the stator of an 8 mm bulk PZT ultrasonic motor. A rotor would be pressed 
against the stator. The stator structure shown here has an inner diameter of 5 mm with a ring of 
piezoceramics bonded to the bottom in a four-wavelength electrode pattern. 
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u z (Microns) u z vs. Freq 



Freq (kHz) 

Figure 7-2: Interferometer readings of one place on the stator as the frequency is swept. The 
amplitude of out of plane deflections as the traveling wave moves through the stator is under 1 pm 
for this 8 mm motor. 

of these 8 mm diameter ring-type motors is shown in Figure 7-1. The idea is to study these 
devices, understand the phenomena at the rotor-stator interface, discern design tradeoffs, 
and then fabricate a new set of motors using improved tliin-film PZT developed, in parallel 
with these efforts, at the Pennsylvania State University’s Materials Research Laboratory. 
Achieving large-area electrode coverage of tliin-hlm PZT which is free of pinholes is difficult, 
lienee the decision to keep these motors rather small. An 8 mm diameter motor was chosen 
as a valid tradeoff between ease of handling and testing, and the difficulty of fabricating 
large-area tliin-hlm PZT devices. 

With these bulk 8 mm motors, we have been able to study and measure a number 
of important characteristics of ultrasonic motors, such as impedance profiles, out-of-plane 
deflections and speed-torque curves for varying normal forces. 

In order to determine which mode shapes correspond to which frequencies, we have 
run finite element simulations and also used an interferometer at the Penn State Materi¬ 
als Research Laboratory [Zliang, Pan and Cross 88] for measuring deflections. Figure 7-2 
illustrates the results of interferometer measurements at one position on the stator as the 
frequency was swept from 50 kHz through 63 kHz. A resonance occurs near 53 kHz and the 
maximum deflection is under 1 /mi. This graph shows clearly that on either side of reso¬ 
nance, the displacement falls off. Furthermore, this data emphasizes the need for smooth 
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Magnitude of Displacement Around the Stator 



Figure 7-3: Interferometer readings around the stator at one frequency show non-uniform ampli¬ 
tudes of displacements. For traveling wave excitation, all points should see the same amplitude of 
deflection. 

interface surfaces. We have found that glass lenses spin much better than brass rotors 
turned on a lathe, even after polishing. Subsequent Dektak measurements of surface rough¬ 
ness gave 300 A average surface roughness for the glass rotor, as opposed to 7000 A average 
surface roughness for the brass rotor. For out-of-plane displacements of just under 1 /mi, 
the 7000 A asperities of the brass rotor nearly touch the valleys of the traveling wave. 

Interferometric measurements were also performed at one frequency but at equally 
spaced points around the perimeter of the stator. Figure 7-3 graphs the resulting mea¬ 
surements of the amplitude of out-of-plane deflections. For an ultrasonic motor working 
on a traveling wave principle, all points along the stator should see equal amplitudes of 
displacement. Here, we see non-uniform displacements, a clear area for improvement. This 
anomaly could be due to the electrode pattern, which is non-uniform, or to an imperfect 
bond. 

For measuring speeds and torques under varying normal forces, we contracted to C'ete- 
lior, the French watch industry’s center for research and development in horology, to build 
a dynamometer to fit these 8 mm motors. Shown in Figure 7-4, the dynamometer is capable 
of measuring torque from 0 to 60 X 10 -4 Nm, and normal forces from 0 to 20 N. An optical 
encoder delivers 6 pulses per motor turn. 

When running the motor in the dynamometer with a glass rotor, as the side loads are 
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■ CETEHOR M 


Figure 7-4: This dynamometer, built by Cetehor, measures speeds and torques along a range of 
operating points and for varying normal forces. The steel stator, with wires attached, sits just below 
a glass rotor. Above the rotor, is an optical encoder with a wrapped string acting as a friction brake. 
The string wraps around a leaf spring at the left and attaches behind to an LVDT displacement 
transducer. The load on the motor is changed by the knob on the left, which pushes on the leaf 
spring. The normal force is adjusted by the knob at the top of the picture which acts on a spring 
pushing the rotor against the stator. 


increased and the motor delivers more torque, the glass rotor tends to wear and get eaten by 
the steel stator. It is important then, not only to have very smooth surfaces, but also hard, 
wear-resistant surfaces. After switching to sapphire rotors, we have been able to measure 
speed-torque curves for five operating points as shown in Figure 7-5. The no-load speed has 
been found to be approximately 230rpm and the stall torque was measured to be 3.0gf-cm 
(3.0 X 10 -4 Nm). 

All of these experiments yielded insight to the workings of these ultrasonic motors, but 
this was only a preliminary pass and an exercise in putting together the appropriate instru¬ 
mentation for characterizing such small motors. Unfortunately, just as the experimental 
setup came together, the motors all broke. They delaminated. Having been only assembled 
with superglue, the extensive testing on the interferometer while driven at 70 V peak caused 
the motors to heat up such that the initial stall torques first measured became so small as 
to be unmeasurable. 
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Figure 7-5: Speed and power as functions of torque are shown here for measured operating points 
of an 8 mm steel stator with sapphire rotor using the Cetehor dynamometer. 


7.2 Second Generation 8 mm Bulk Motors 


After the initial batch of motors passed away, we undertook fabrication of a second gener¬ 
ation of motors incorporating improvements from lessons learned in the first experiments. 
The goal was to produce as complete a data-set as possible, both to verify proposed models 
and to provide useful open-literature real data, to other researchers in the held. 

The second generation of motors was a set whose parameters were chosen using design- 
of-experiments techniques. Design-of-experiments methods are another approach to cre¬ 
ating and fitting a model: a statistical-experimental investigation. The large number of 
parameters involved in the design of an ultrasonic motor creates a huge design space. Size, 
geometry, number of teeth, tooth height, base height, type of ceramic, epoxy composition, 
poling held, drive voltage, drive frequency, electrode pattern, stator material, surface finish, 
rotor material and normal force are just a few of the parameters to be varied. 

Consequently, an experiment consisting of 32 trials, or combinations of stators and 
rotors, was set up with stall torque and no-load speed taken as the quality parameters to be 
measured. All of the stators maintained the 5 mm inner diameter and 8 mm outer diameter 
dimensions, but the stator designs varied 4 parameters in 8 combinations. Three different 
rotor parameters were combined with these basic stator types. The way in which these 7 
parameters were varied in the 32 trials was based on design-of-experiments methods which 
enable problems of '2 n configurations to be studied with approximately 2n experiments, 
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Figure 7-6: This is the second generation electrode pattern used on the set of 8 mm bulk PZT 
stators. The pattern at left is laser cut, then poled and reconnected to produce the pattern at right. 
In this manner, the entire electroded surface is utilized and the final number of wires remains three. 

where n is the number of parameters. 

7.2.1 Bonds and Electrodes 

From the preliminary tests run on the first batch of 8 mm motors, we learned of a few 
problems which we addressed in all motors in the second generation motors. The most 
obvious realization was that a higher temperature epoxy was needed. After a number 
of trials, we found Masterbond EPTS-10 conductive epoxy to produce a strong laminate. 
The second problem addressed was the 11011 -uniformity of out-of-plane displacements. The 
electrode pattern used for the first set of 8111 m motors was a simple pattern that could be 
etched by hand with a surgeon’s scalpel, poled and reconnected with silver paint, as shown 
earlier in Figure 3-6. However in that design, the two phases of standing waves were set up 
011 opposing halves of the stator and spacer segments of one wavelength were required. For 
a four-wavelength motor, fully 25% of the available electrode area was undriven. 

I 11 the second generation motors, we switched to a more sophisticated electrode design 
which leaves 110 undriven ceramic and which interlaces the two phases of standing waves. 
The new electrode pattern used is shown in Figure 7-6 and was contracted for from Piezo 
Systems, Cambridge, MA, and uses PTS-5A ceramic. The pattern 011 the left is first pat¬ 
terned by laser machining with line widths of 1.95mm (7.5mils). The same width spaces 
are left at the inner and outer edges. Pairs of neighboring segments are then poled and are 
reconnected after poling to produce the final pattern as shown on the right in Figure 7-6. 
The inner set of segments is then driven 90 degrees out of phase from the outer set of 
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segments to produce a traveling wave. Each segment is one quarter wavelength long and so 
this pattern of 12 segments produces a three-wavelength ultrasonic motor. First tests with 
this electrode pattern bonded onto a steel stator using Masterbond conductive epoxy gave 
promising results. Stall torque was measured to be 4.6gf-cm, 50% higher than any previous 
measurement. 

7.2.2 Design of Experiments 

Due to the large number of parameters involved in ultrasonic motors, we used statistical de¬ 
sign of experiments methods in our second batch of motors. Design-of-experiments methods 
are the planning of experiments so that appropriate data is collected which can be analyzed 
by statistical methods for formulating objective conclusions [Montgomery 91]. 

The term “experiment” relates to the entire set of 32 trials, or runs. There are 7 
parameters, or factors, chosen to be varied: 4 pertaining to the stator and 3 pertaining to 
the rotor. The 4 stator parameters are the tooth height, the base height, the number of 
teeth and the stator material. For the rotor, we vary the rotor material, the rotor liner and 
the lubrication. Voltage and phase are kept constant throughout all the runs. Each motor 
is run at its third mode resonant frequency and the normal force is adjusted to extract the 
maximum stall torque or no-load speed. 

We use what is known as factorial design, specifically a factor of two design. That is, 
each parameter is assigned one of two binary possibilities, denoted as ‘ + ’ and The stator 
and rotor parameters are assigned the values shown in Figure 7-7. 

The stator tooth height is either 30 mils (0.76 mm) or 50 mils (1.3 mm). The base 
height of the metal below the teeth is either 10 mils (0.25 mm) or 15 mils (0.38 mm). The 
number of teeth in the stator is either 24 teeth or 36 teeth, and the stator material can be 
either aluminum or stainless steel. 

The rotor parameters are the material composing the rotor (sapphire or 01 tool steel 
hardened to a Rockwell C scale hardness of 64), the rotor liner (no liner or kapton tape) 
and the lubrication (no oil or WD-40 oil) 

For seven parameters, 2 7 possible combinations result. A fractional factorial design 
selects a subset of these possibilities and aliases some of the effects. The idea is to alias 
the main effects with higher-order interaction effects which can be considered negligible. 
Interactions that can be assumed, a priori, to have significant effects may be left in the 
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Parameters 


- 

+ 

1 

Tooth Height 

30 mils (0.76 mm) 

50 mils (1.3 mm) 

2 

Base Height 

10 mils (0.25 mm) 

15 mils (0.38 mm) 

3 

#Teeth 

24 

36 

4 

Stator Material 

6061 Aluminum 

304 Stainless Steel 

5 

Rotor Liner 

None 

Kapton 

6 

Rotor Material 

Hardened 01 Tool Steel 

Sapphire 

7 

Lubricant 

None 

WD-40 Oil 


Figure 7-7: The stator and rotor parameters in the designed experiment are one of two possibilities 
each. 

design set unaliased and accounted for in the analysis. 

If we seek to perform a smaller number of runs than the 2' possible, say 32, we can 
formulate a one-quarter fraction factorial design, denoted as a ' 2 k ~ 2 design, where in this 
case k = 7. Figure 7-8 illustrates the appropriate mapping of the 2 5 combinations. The 
first three factors, tooth height, base height and the number of teeth (denoted rq, rq and 
i >3 respectively) are assigned all possible (2 3 ) combinations of levels. The fourth factor, the 
stator material, is assigned to the combination that is equivalent to the interaction between 
v -2 and 113 , which is assumed insignificant. Note that column V 23 is the product of columns 
iq and 113 . I 11 this way, four parameters can be studied with 8 runs. Figure 7-8 illustrates 
the designed experiment. 

If it is judged that some interactions will be significant, they are explicitly left in ad¬ 
ditional columns of the matrix, to be taken into account in the analysis of the measured 
response. For the response, we measure both stall torque and no-load speed. The average 
column is also included to account for common mode results in each trial. 

Each row of Figure 7-8 represents a specific motor trial. Row 1 signifies a motor trial 
where the stator has a tooth height of 50 mils, a base height of 15 mils, 36 teeth and is 
made out of steel. The rotor is sapphire with a kapton liner and has a WD-40 oil lubricant. 
Row 2 signifies the same stator and same rotor, but no lubrication, and so on. 

Note that the stator material parameter, S ma t or C 23 , is aliased with the interaction 
between rq and 113 . Thus we have assumed the interaction between the base height and the 
number of teeth is insignificant. 

For each trial, a response was measured: maximum stall torque and maximum no-load 
speed, as the normal force was varied. The set of motors fabricated is shown in Figure 7- 
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8 mm Bulk PZT Ultrasonic Motor Experiment Design 
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Figure 7-8: This one-quarter fraction factorial design selects 32 trials out of a possible 2 '. The 
Trow blocks, given names as the colors violet, fucia, tan, ffiblue, green, yellow, orange and grey, 
signify the 8 different stators that we have manufactured. 




Kfc, 

Figure 7-9: This set of eight different 8 mm diameter stators was fabricated for the designed 
experiment. Two different rotors, with and without liners and/or lubricants combined to make a 
32 -trial experiment. 

9. Design-of-experiments methods provide for a way of analyzing the data and computing 
both quantitative and qualitative conclusions. [Fiegutli, Spina and Staelin 94] describe a 
software tool for assisting in setting up the design matrix and in processing the results of 
the trials. A regression analysis is performed on the matrix of plus ones and minus ones, 
the vector of unknown factor effects and the vector of measured responses. Thus a number 
is assigned to each of the columns and is the weight in the resulting predictor polynomial 

for the optimum response. A qualitative result, a sequence of -|-1-b + values is also 

determined. This string refers to the optimum combination of variables for obtaining the 
minimum response, (negative signs are simply input with the measured stall torques and 
no-load speeds in order to find the maximum response) and this combination may not be 
one of the trials. 

7.2.3 Fabricating Bulk PZT 8 mm Motors 

In order to perform these experiments with the new electrode pattern, a more sophisticated 
process had to be developed for etching the electrode pads. Laser machining offered the 
potential of a fast, mask-free, process which is ideally suited for building prototypes and 
changing designs quickly. However, significant amounts of testing had to be done to find 
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Figure 7-10: PZT-5A material has been cut. and etched by Piezo Systems, Cambridge MA, with 
the electrode pattern required for poling. The dot signifies that pad was poled with positive voltage 
on that side. 

the appropriate laser parameters for both cutting through the thickness of the PZT and for 
etching the top electrode without burning the material. 

We did extensive tests with both YAG and excimer lasers through Laser Services, West- 
ford MA, with limited success. Both laser processes tended to chip and burn the ceramic 
at the edges of the cuts. Nickel-coated PZT-5A from Piezo Systems, Cambridge MA, was 
somewhat easier to machine than silver-coated PZT-4 from Transducer Products, Goshen 
C'T, as the silver electrode tended to melt at a lower temperature. 

We contracted to Piezo Systems to provide PZT-5A cut in the 8 mm rings and etched 
with the electrode pattern of Figure 7-6. They used a proprietary laser-machining process 
and the results are shown in Figure 7-10. These parts are free from burn marks and the 
ceramic has been successfully poled and used in the second generation of motors. 

After this pattern is poled, we reconnect segments by painting conductive epoxy over 
the appropriate etch. This was done by hand, and details of the procedure and other 
aspects of the manufacturing process are given in [Franck 95]. We also experimented with 
a laser direct-write system with the assistance of Dr. Daniel Ehrlich at the MIT Lincoln 
Laboratory Solid State Division. The laser-direct writing process was also able to perform 
the reconnect. 

Once the electrode interconnect was complete, wires were attached to each of the inner 
and outer sets of pads. A third wire was soldered to thffi side of the stator for a ground 
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contact. Stainless steel stators were not troublesome in this regard, but the aluminum 
stators were. In order to solder to the aluminum stator, we conductive-epoxied brass shims 
to one small patch of the stator and then successfully made solder contacts to the shims for 
ground leads. 

The stator teeth are polished on very fine sandpaper to make them all of uniform surface 
finish. We performed Dektak experiments to characterize the surface finishes of both the 
stators and the rotors and results are shown later. 

Once the stators were completed, the rotors were placed on top of the stators and the 
combination was then placed into the dynamometer shown in Figure 7-4. The metal head 
above the optical encoder can be moved up vertically to allow for positioning the motor. 
Then the head is lowered and normal forces are adjusted by the normal force knob at the 
top of the dynamometer. The dynamometer provides the bearing for the motor. That 
is, the rotor and stator are two separate pieces and do not work as a motor outside of the 
dynamometer. We have postponed developing bearings and building up discrete component 
motors until the surface interactions and normal force conditions are better characterized. 
This technique of using the dynamometer as the bearing allows easy switching between 
rotor and stator combinations for the 32 trials in the experiment. 

7.3 Results From Second Generation 8 mm Motors 

7.3.1 Design of Experiments Results 

The 32 trials in the experiment were run on the Cetehor dynamometer and the normal force 
and side loads were adjusted to find the maximum stall torque and no-load speeds. Each 
motor was run at its resonance frequency with a peak voltage of 60 V and the two phases 
of the drive signal were set 90 degrees apart. 

The results of the experiment are tabulated in Figure 7-11. It was found that the motors 
did not work at all when WD-40 lubricant was applied. Consequently, every other row in the 
table has a zero result. This means that the matrix can be recast as a 16-trial experiment 
where the oil parameter column is removed as shown in Figure 7-12, where the frequency 
and normal forces are also noted. 

Figure 7-12 shows that stator Green (0.76mm tooth height, 0.38mm base height, 36 
teeth, made out of stainless steel), in combination with an unlined sapphire rotor and no 
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Design of Experiments Results: Stall Torques and No-load Speeds 



Stator 

Liner 

Rmat 

Oil 

Max Stall Torque 
(gf-cm) 

Max No-load Speed 
(rpm) 

1 

Violet Exp 1 

+ 

+ 

+ 

0 

0 

2 

Violet Exp 2 

+ 

+ 

- 

4.7 

1660 

3 

Violet Exp 3 

+ 

- 

+ 

0 

0 

4 

Violet Exp 4 

+ 

- 

- 

4.8 

1600 

5 

Fucia Exp 1 

+ 

+ 

+ 

0 

0 

6 

Fucia Exp 2 

+ 

+ 

- 

5.5 

1750 

7 

Fucia Exp 3 

+ 

- 

+ 

0 

0 

8 

Fucia Exp 4 

+ 

- 

- 

4.5 

1610 

9 

Tan Exp 1 

- 

+ 

+ 

0 

0 

10 

Tan Exp 2 

- 

+ 

- 

2.7 

746 

11 

Tan Exp 3 

- 

- 

+ 

0 

0 

12 

Tan Exp 4 

- 

- 

- 

3.4 

893 

13 

Blue Exp 1 

- 

+ 

+ 

0 

0 

14 

Blue Exp 2 

- 

+ 

- 

4.7 

787 

15 

Blue Exp 3 

- 

- 

+ 

0 

0 

16 

Blue Exp 4 

- 

- 

- 

4.6 

813 

17 

Green Exp 1 

- 

+ 

+ 

0 

0 

18 

Green Exp 2 

- 

+ 

- 

10.0 

870 

19 

Green Exp 3 

- 

- 

+ 

0 

0 

20 

Green Exp 4 

- 

- 

- 

8.4 

781 

21 

Yellow Exp 1 

- 

+ 

+ 

0 

0 

22 

Yellow Exp 2 

- 

+ 

- 

5.3 

800 

23 

Yellow Exp 3 

- 

- 

+ 

0 

0 

24 

Yellow Exp 4 

- 

- 

- 

5.6 

840 

25 

Orange Exp 1 

+ 

+ 

+ 

0 

0 

26 

Orange Exp 2 

+ 

+ 

- 

4.8 

1060 

27 

Orange Exp 3 

+ 

- 

+ 

0 

0 

28 

Orange Exp 4 

+ 

- 

- 

4.6 

1100 

29 

Grey Exp 1 

+ 

+ 

+ 

0 

0 

30 

Grey Exp 2 

+ 

+ 

- 

4.3 

885 

31 

Grey Exp 3 

+ 

- 

+ 

0 

0 

32 

Grey Exp 4 

+ 

- 

- 

5.3 

901 


Figure 7-11: Resulting stall torques and no-load speeds for the designed experiment. 





Design of Experiments Results: Stall Torques and No-load Speeds 



Stator 

Liner 

Rmat 

Frequency 

(kHz) 

Max Stall Torque 
(gf-cm) 

@FN 

(gf) 

Max No-load Speed 
(rpm) 

@FN 

(gf) 

1 

Violet Exp 2 

+ 

+ 

44.5 

4.7 

62 

1660 

6.9 

2 

Violet Exp 4 

+ 

- 

44.5 

4.8 

66 

1600 

11 

3 

Fucia Exp 2 

+ 

+ 

55.5 

5.5 

81 

1750 

11 

4 

Fucia Exp 4 

+ 

- 

55.5 

4.5 

97 

1610 

4.8 

5 

Tan Exp 2 

- 

+ 

38.0 

2.7 

38 

746 

4.1 

6 

Tan Exp 4 

- 

- 

38.0 

3.4 

40 

893 

21 

7 

Blue Exp 2 

- 

+ 

35.5 

4.7 

76 

787 

35 

8 

Blue Exp 4 

- 

- 

35.5 

4.6 

83 

813 

14 

9 

Green Exp 2 

- 

+ 

57.7 

10.0 

110 

870 

6.9 

10 

Green Exp 4 

- 

- 

57.7 

8.4 

110 

781 

30 

11 

Yellow Exp 2 

- 

+ 

63.0 

5.3 

66 

800 

6.9 

12 

Yellow Exp 4 

- 

- 

63.0 

5.6 

59 

840 

2.8 

13 

Orange Exp 2 

+ 

+ 

47.4 

4.8 

97 

1060 

17 

14 

Orange Exp 4 

+ 

- 

47.4 

4.6 

97 

1100 

14 

15 

Grey Exp 2 

+ 

+ 

42.2 

4.3 

94 

885 

17 

16 

Grey Exp 4 

+ 

- 

42.2 

5.3 

110 

901 

14 


Figure 7-12: Resulting stall torques and no-load speeds for the designed experiment with the oil 
parameter removed. Stator Green, in combination with an unlined sapphire rotor and no lubricant, 
produced the largest stall torque: lOgf-cm. Stator Fucia, with a kapton-coated sapphire rotor and 
no lubricant, produced the highest no-load speed: 1750rpm. 

lubricant, produced the largest stall torque: lO.Ogf-cm. Stator Fucia (1.3mm tooth height, 
0.38 mm base height, 24 teeth, made out of aluminum), with a kapton-coated sapphire rotor 
and no lubricant, produced the highest no-load speed: 1750rpm. 

The factorial design model represents a quality measure of a system as a linearly weighted 
sum of the parameters of the system where the parameters are valued as either +1 or -1. 
For instance, the quality measure for each row of the experiment shown in Figure 7-8 would 
be represented as: 


Quality Measure=TT a „ 3 v avg + TFii’i + W 2 v 2 + W 3 v 3 
+W 23 v 23 + W 4 V 4 + W 5 v 5 + W 6 v 6 


where v n is +1 or -1, depending on the trial, and the weights, W n , are found from a 
regression analysis after performing all trials (measuring stall torques, for instance). 

The analysis was performed using a software tool [Fiegutli, Spina and Staelin 94] devel¬ 
oped in the MIT Leaders for Manufacturing Program which not only determines the weights, 
W n , for the predictor polynomial, but also performs a search through the parameter space 
for the optimal combination of parameters to minimize the quality factor measured. 
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For the maximum stall torque quality measurement, the analysis results are as follows 
(maximum stall torque data were input as negative values for maximizing the result): 


Experiment Analysis Program. 

Matrix Identification Code: tes!14 


Regression 

Coeff 

-5.200000 

0.837500 

-0.900000 

-0.225000 

-0.650000 

0.387500 

-0.050000 

0.687500 

0.312500 

0.012500 

-0.225000 

0.037500 

-0.075000 

0.312500 

-0.025000 

0.062500 


Matrix 

Column 

Average Quality 

varl ( var2-var4 ) 

var2 ( varl-var4 var3-var23 ) 

var3 ( var2-var23 ) 

var23 ( var2-var3 ) 

var4 ( varl-var2 ) 

var5 

(unassigned) ( varl-var3 var23-var4 ) 

(unassigned) ( varl-var23 var3-var4 ) 

(unassigned) ( varl-var5 ) 

(unassigned) ( var2-var5 ) 

(unassigned) ( var4-var5 ) 

(unassigned) ( var3-var5 ) 

(unassigned) 

(unassigned) ( var23-var5 ) 
(unassigned) 


Optimal Parameter Assignments: 
varl var2 var3 var23 var4 var5 

- + + + - + 


Optimized Quality: -8.250000 


Free Parameters: 9 

RMS Unconfidence: 0.284861 


The predictor polynomial for maximum stall torque then is: 


t stall,max = —5.200 + 0.8375tq — 0.9000r 2 — 0.2250r 3 — 0.6500u 23 + 0.3875r 4 — 0.0500r 5 

where the v n represent +1 or -1 values for tooth height, base height, the number of teeth, 
the stator material, the rotor liner and the rotor material, respectively, for each trial motor. 
Plugging in any other combination of parameters will yield a predicted value of stall torque 
based on the regression analysis of the prototypes fabricated. 

The analysis also performs a search through the space of combinations and find the 
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combination “—b + H—b”, which refers to stator Green, with an unlined sapphire rotor as 
the optimal configuration for maximizing stall torque. This combination happens to be 
one of the trials performed, but in general the predictor string may point to an untried 
combination. 

The maximum no-load speeds were also analyzed and the results were found to be: 


Experiment Analysis Program. 

Matrix Identification Code: tes!14 


Regression Matrix 
Coeff Column 


-1068.500000 

-163.875000 

-170.375000 

-20.250000 

31.375000 

-252.250000 

-1.250000 

27.875000 

-14.000000 

-2.125000 

-29.875000 

-16.750000 

6.000000 

19.125000 

-12.125000 

7.000000 


Average Quality 
varl ( var2-var4 ) 
var2 ( varl-var4 var3-var23 ) 
var3 ( var2-var23 ) 
var23 ( var2-var3 ) 
var4 ( varl-var2 ) 
var5 

(unassigned) ( varl-var3 var23-var4 

(unassigned) ( varl-var23 var3-var4 

(unassigned) ( varl-var5 ) 

(unassigned) ( var2-var5 ) 

(unassigned) ( var4-var5 ) 

(unassigned) ( var3-var5 ) 

(unassigned) 

(unassigned) ( var23-var5 ) 
(unassigned) 


) 

) 


Optimal Parameter Assignments: 
varl var2 var3 var23 var4 var5 

+ + + - + + 

Optimized Quality: -1707.875000 

Free Parameters: 9 

RMS Unconfidence: 17.475031 


Here, the predictor polynomial for maximizing no-load speed turns out to be: 


N no _i oa d,max = —1069 — 163.9fi — 170.4t> 2 — 20.25r 3 + 0.31.38r 2 3 — 252. 3v 4 — 1.25r 5 

In this case, the predictor variable string is “+ + -1—b+”, which is a different combination 
than any of the trial motors built. This suggests that the next experiment would be to 
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fabricate such a stator: tall teeth, large base height, large number of teeth, made out of 
aluminum, with a kapton-coated sapphire rotor, and check. Additionally, replicates should 
be made of Fucia, as Fucia Experiment 2 actuahy gave a higher no-load speed, 1750 rpm, 
than the predictor polynomial for the indicated optimal parameter string. 

7.3.2 Speed-Torque Curves 

Two complete replicate sets of the eight stators have been fabricated, but the second set 
has not yet been tested. Of the first replicate set, all motors were measured for maximum 
no-load speed and maximum stall torque, as discussed in the previous section. Here we 
refer to stators in the first replicate set as Violet 1, Fucial, Tanl, and so on. Stators in the 
second replicate set will be referred to as Violet2, Fucia2, Tan2, etc. 

In addition to measuring maximum no-load speeds and maximum stall torques, three 
trial motors of the first replicate set have been tested across a spectrum of operating points: 
Yellowl Experiment 2 (unlined sapphire rotor), Greyl with an unlined sapphire rotor and 
Fucial, also with an unlined sapphire rotor. These complete speed-torque curves are shown 
below in this section for a constant drive voltage of 60 V peak as the normal force is varied. 

7.3.3 Yellowl 

Figure 7-13(a) illustrates four speed-torque curves for Yellowl Experiment 2. The Yellowl 
stator is an aluminum stator of 24 teeth, with a tooth height of 0.76 mm and a base height 
of 0.38 mm. The rotor in this run was an unlined sapphire rotor. 

It can be seen in Figure 7-13(a), that as the normal force is increased, the no-load speeds 
fall off and the stall torques increase. For a normal force of 55 gf, the stall torque is 5.2 gf-cm. 
However, as the normal force is increased to 69 gf, as shown in Figure 7-13(b), the no-load 
speed continues to drop but the stall torque does not increase. It falls to 4.6 gf-cm. No- 
load speed versus normal force is plotted in Figure 7-14(a) and stall torque versus normal 
force is plotted in Figure 7-14(b). The no-load speeds falloff as predicted in Figure 6-21, 
but the stall torque was predicted in Figure 6-20 to reach an asymptotic value, under the 
assumption of a rigid stator. Here we see the effect of stator squashing come into play for 
large normal forces and the stall torque decrease with increasing normal force. 

By multiplying speed and torque we can find the power versus torque profiles. Figure 7- 
15 shows that the peak mechanical output power for Yellow Experiment 2 is 12 mW. 
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Figure 7-13: (a) The first, four speed-torque curves of Yellow Experiment 2 (stator Yellowl run 
with an unlined sapphire stator) are shown here with second order regression polynomial curves fit 
onto the scatter plot, (b) As the normal force is increased from 55 gf to 69 gf in Yellow Experiment 
2, the stall torques stops increasing and begins to decrease. 
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Figure 7-14: (a) The no-load speed versus normal force curve shows a dropoff as predicted, (b) 
Stall torques versus normal force are plotted here for Yellow Experiment 2. The decrease in stall 
torque for large normal forces indicates squashing of the stator. 
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Figure 7-15: Power versus torque is plotted with curve fits. The maximum mechanical output 
power is 12 mW and occurs at a normal force of 28 gf. 


7.3.4 Fucial 

The largest no-load speeds in the trial motors of the designed experiment were exhibited by 
stator Fucial which was made of aluminum, had 24 teeth, a tooth height of 1.3 mm and a 
base height of 0.38 mm. Although the designed experiment specified a kapton liner for both 
Fucial Experiment 2 and Fucial Experiment 4, we also experimented with unlined rotors 
and found that plain sapphire rotors produced very high speed motion. 

Figure 7-16 illustrates speed-torque curves for data taken at various normal forces for 
operating points closer to no-load speed than stall, in order to measure output power quickly. 

Figure 7-17 depicts the resulting power vs. torque curve for Fucial. Second order 
polynomial regression curve fits are also shown. The peak power is produced at a normal 
force of 28 gf and was measured to be 27 mW. 

To recognize the significance of this amount of output power, it is interesting to compare 
this ultrasonic motor to the smallest commercially available component electromagnetic 
DC motor. In the MIT Mobile Robot Group at the Artificial Intelligence Laboratory, we 
typically use Namiki DC motors for small mobile robot applications. These motors are 
7 mm in diameter and 16 mm tall. They are specified [Namiki 92] to have a stall torque 
of 1.8gf-cm and a no-load speed of 12,000rpm. Taking output power as one-quarter of 
the product of no-load speed and stall torque gives a peak output power of 57 mW. Thus 
this ultrasonic motor, Fucial with a plain sapphire rotor, produces approximately half the 
output power of a Namiki motor in 1/6 the volume. Even more importantly, this power 
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Figure 7-16: Speed vs. torque curves for Fucial run with an unlined sapphire rotor (not a trial in 
the designed experiment) show very large no-load speeds. The peak speed here is 2080rpm. 



Torque (gf-cm) 


Figure 7-17: The power vs. torque curves for Fucial running with an unlined sapphire rotor 
give the largest mechanical output powers of any motors fabricated here. Peak output power was 
measured to be 27 rnW at a normal force of 28 gf. 
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is delivered at a much lower speed and higher torque, a better impedance match for most 
real-world applications. 

7.3.5 Greyl 

The motor we have collected the most extensive data set on is stator Greyl with a plain 
sapphire rotor. Speed-torque curves were measured at seven different normal forces ranging 
from Tjv=28gf to 97 gf. The drive voltage was again 60 V peak and the frequency was 
42.2kHz as noted earlier in Figure 7-12. 

If we look back to the speed-torque curves of Figure 7-5, we can compare the number 
of operating points we were able to resolve on the early prototypes with the number of 
operating points and speed-torque curves on this second generation batch. The stall torques 
are approximately a factor of three larger and half a dozen speed-torque curves at different 
normal forces are resolvable before the stator is completely squashed down. This is likely due 
to larger displacements of the stator resulting from improved bonds and electrode patterns. 

The first four speed torque curves are shown in Figure 7-18(a). The no-load speeds 
decrease with normal force as the stall torques increase, similarly to Figure 7-13(a). As 
the normal force is increased from 69 gf to 83 gf, the stall torque grows from 8.4gf-cm to 
8.9gf-cm. However, as the normal force is increased further to 97gf-cm, the stall torque 
falls off to 7.9 gf-cm. 

The stall torques and no-load speeds are plotted versus normal force in Figure 7-19. The 
curves develop in much the same manner as for Yellowl Experiment 2 shown in Figure 7- 
14. Again, we see no-load speed drop off with increased normal force. Stall torques also 
increase with normal force up until a point where stator squashing causes the stall torque 
to decrease. 

The mechanical output power curves can then be calculated and are shown in Figure 7- 
20. The peak output power is 16 mW at a normal force of 69 gf. The speed-torque operating 
point is 368 rpm at a load torque of 4.2 gf-cm. Again, this low-speed, high-torque operating 
point at maximum power compares favorably with electromagnetic motors such as the 
Namiki DC motor which is commonly used as a vibration alarm in silent pagers. In the 
pager application, a geardown is not required, simply an eccentric cam to promote shaking. 
In fact, Namiki motors are not even available with gearboxes. Gearing them down to useful 
torques and speeds to drive a small robot requires building custom gearboxes. 
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Figure 7-18: (a) The first four speed-torque curves for motor Greyl with an unlined sapphire 
rotor, with regression curves, (b) Fifth, sixth and seventh speed-torque curves for motor Greyl 
show maximum stall torques of 8.9 gf-cm. 















Figure 7-20: (a) The first four power versus torque curves for motor Greyl. (b) Fifth, sixth and 
seventh power versus torque curves for motor Greyl. Peak output power of 16mW occurs at a load 
torque of 4.2gf-cm under a normal force of 69 gf. 














Figure 7-21: (a) The plot of the magnitude of impedance versus frequency shows a resonance at 
43.6kHz and an antiresonance at 44.4kHz for the inner electrode, (b) The phase of the impedance 
shows a resonance near 44 kHz and capacitive characteristics away from resonance. 

7.3.6 Magnitude and Phase of Impedance and Damping Measurements 

All the speed-torque curves just displayed for stator Greyl were measured at 42.2kHz, 
the frequency which produced the largest stall torque. Due to the piezoelectric effect, the 
mechanical resonance can be viewed electrically, as was discussed in Figure 4-1. 

A Hewlett Packard 4192A Low Frequency Impedance Analyzer, which was frequency 
swept via General Purpose Interface Bus control from an IBM 755C Tliinkpad laptop run¬ 
ning Labview instrumentation software, was used to measure the magnitude and phase of 
both the inner and outer electrodes on stator Greyl as shown in Figure 7-21. 

The magnitude of the impedance reaches a minimum at 43.6 kHz and a maximum at 
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Figure 7-22: The turn-off transient of stator Greyl from a GOV peak drive at 42.2 kHz displays an 
initial increase in voltage due to the piezoelectric effect, and than a damped response down to zero. 

44.4 kHz for the inner electrode. The outer electrode is closely matched. Initial prototypes 
with the first electrode pattern often did not display such uniformity, most likely due to 
the fact that the ceramics then were ground by hand and etched with a surgeon’s scalpel. 
These laser-etched ceramics, even though the two electrode areas are very slightly different, 
display a better matching of resonance conditions. 

The damping of the unloaded stator can be measured either from the frequency response 
or from the time-domain transient response. Figure 7-22 illustrates the turn-off transient 
from an initial excitation of GOV peak at 42.2 kHz. The quality factor, Q, can be measured 
from the turn-off transient. Q is defined as the number of radians in the signal in the time 
it takes to reach ^ of its initial value [Staelin, Morgentlialer and Kong 94]. 

Figure 7-22 shows that the two drive signals, sin and cos, are not perfectly matched 
in their damping characteristics. The channel driving the outer electrode is trace 1 at the 
top of the oscilloscope display. Using the zoom feature to examine the waveform, it was 
found that the signal decayed to ^ in 240 ^s. At 42.2kHz, or 265,000 the Q of this 
signal works out to be on the order of 63. The damping coefficient, £, defined by Q = -^ 
is then £ = 0.8% For the second trace, at the bottom of the oscilloscope display, the signal 
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Figure 7-23: The interferometer used to measure displacements on the 8 mm diameter stators. 
The stator is mounted on the vertical X-Y stage at left such that the laser spot hits the surface of 
a stator tooth for measuring out-of-plane deflection. 


decays to ^ of its value in 460 /.is, yielding a Q of 122 and a damping coefficient of £ = 0.4%. 
The value of Q is the dynamic amplification factor by which the static deflection will be 
multiplied at resonance. As the two electrode patterns yield different quality factors, we 
can expect that each electrode pattern will contribute different amplitudes of deflection to 
the traveling flexure wave. 


7.3.7 Interferometric Displacement Measurements 

To get further insight into the mechanisms underlying torque production in ultrasonic mo¬ 
tors, we put together a Michelson interferometer to measure out-of-plane displacements of 
the stators. The interferometer is illustrated in Figure 7-23. The motor can be seen at 
left mounted on a vertical X-Y stage. The laser is at far right, a polarized beamsplitter 
with two attached quarter-wave plates is shown at the center and an adjustable polarizer 
is shown off to center-right. The detector is off the photo beyond the adjustable polarizer. 
Also out of the picture is the reference-arm mirror of the interferometer which would be 
across from the adjustable polarizer. A microscope objective is usually placed between the 
beamsplitter and the stator to focus the beam onto the small surface area of a stator tooth, 
but is left out of the photo so that the stator can be seen. A spatial filter and collimating 
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Figure 7-24: Interferometer measurements of displacement on stator Greyl at 11 o’clock. The 
stator is driven at GOV peak and at 42.2kHz. The big blips in the lower trace mark the turn-around 
points from when the stator deflection changes from moving away from the laser beam to moving 
towards it. The fringes in between two big blips then yield the peak-t.o-peak displacement of the 
motor. Notice that the motor has a peak-t.o-peak displacement once per drive period. 


lens lias also been added between the laser and beamsplitter to filter out the structure from 
the laser light. The entire optical setup is mounted on a vibration-isolated optical table in 
a basement laboratory. Our original setup on the ninth floor of the Artificial Intelligence 
Laboratory, with subway trains running underneath, was far too shaky of an environment 
to achieve interference. 

The stator is mounted on the X-Y stage in a free-free manner, much as it is mounted in 
the dynamometer when speed-torque measurements are taken. In Figure 7-23 the motor can 
be seen mounted such that the surface of the stator, which would normally be contacting 
the rotor, is facing the laser spot. The stator is held in place merely by a piece of tape 
across the lead wires. A piece of ka.pton tape is placed below the stator on the X-Y stage 
to prevent electrical shorting by the stage. The only other attachment holding the stator 
fixed is a thin piece of ka.pton tape placed over one tooth of the stator at approximately six 
o’clock. This is also visible in Figure 7-23 as the long piece of tape hanging down off the 
stator. The stator is mounted in this manner so as to permit as much vibration as possible 
while preventing rigid body translation or rotation. 
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Figure 7-25: (a) Interferometer data showing the magnitude of the displacement at one point on 
the stator, at 42.2kHz, as the voltage is swept from 20 V peak to 70 V peak, (b) Interferometer 
data showing the magnitude of the out-of-plane displacement at one location on the motor as the 
frequency is swept from 1 kHz through 90kHz. A resonance is clearly visible. 
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The fringe pattern resulting from driving the unloaded stator is shown in Figure 7-24. 
The top trace is one phase of the stator drive signal and the bottom trace is the output 
of the interferometer’s detector. Fringes between turn-around blips on the trace can be 
counted to find the peak-to-peak out-of-plane displacement of the stator. A helium-neon 
laser of 632.8nm wavelength is used in the interferometer. Half the fringe count times the 
wavelength gives the amplitude of the displacement. Here, the stator is driven at 42.2kHz, 
60 V peak excitation. The laser spot is focused on a tooth positioned at approximately 11 
o’clock. While the outer diameter of the stator is 8 mm and the inner diameter is 5 mm, the 
teeth are canted at such an angle that the actual flat part of the tooth is 0.5 mm wide in 
the radial direction. That is, the inner edge of the top flat part of the tooth is at a radius 
of 3.0 mm from the center of the stator and extends to a radius of 3.5 mm. The spot size of 
the laser is small enough to take about two readings per tooth. 

Figure 7-25(a) illustrates how the out-of-plane deflection changes with applied voltage 
at a frequency of 42.2 kHz. The displacement grows from 1.6 p at 20 V to 4.5 p at 70 V. 

Figure 7-25(b) shows the displacement on the same tooth of the stator, at a drive voltage 
of 60 V peak, as the frequency is scanned. A peak displacement of 5.3p is achieved at 
41.0kHz. 

If we leave the stator drive voltage fixed at 60 V peak and the drive frequency fixed at 
42.2kHz and move the X-Y stage such that the laser beam sequentially measures displace¬ 
ments on each tooth of the stator, we see fringe patterns like those shown in Figure 7-26. 
Whereas Figure 7-24 was the result of focusing the beam on a tooth at approximately 11 
o’clock, Figure 7-26(a) and Figure 7-26(b) illustrate fringe patterns resulting from the beam 
being focused on teeth at 12 o’clock and 1 o’clock around the stator, respectively. By count¬ 
ing fringes, we find displacements of 4.9/um, 4.9/um and 3.9 p at 11-, 12- and 1 o’clock 
respectively. 

Notice that the turn-around blips move to the left, with respect to the reference drive 
signal in the top trace, in subsequent readings from Figure 7-24, Figure 7-26(a) and Figure 7- 
26(b). This phase change signifies where the maximum displacement points are with respect 
to each other as the circumference of the motor is scanned. That is, if we measure the 
phase with respect to the drive signal, for all points around the stator, we can mark off zero 
crossings in phase and visualize the mode of the traveling wave that is generated. 

Figure 7-27(a) plots the amplitude of displacement for points around the stator and 
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Figure 7-27: (a) Interferometer data showing the magnitude of the displacement, at resonance 
and 60 Vpeak drive, around the stator (halfway around and then the interferometer went out of 
alignment). Shows fairly uniform displacement, (b) Interferometer data showing the phase of the 
displacement with respect to the drive signal. This shows which mode is acting, namely the third. 

Figure 7-27(b) shows the phase relative to the drive signal. Every 180° phase shift is marked 
as a zero crossing and denotes half a wavelength of the flexure mode. Unfortunately, the 
interferometer went out of alignment halfway through this experiment, so only half of the 
readings around the stator were measured in Figure 7-27. 

Nevertheless, it can be clearly seen that the magnitudes of the displacements are much 
more uniform for the traveling wave generated here than for the initial prototype motor’s 
displacement shown earlier in Figure 7-3 where the magnitude of deflection sometimes fell 
below 0.1/im. 

Here, displacements average roughly 4.5 /mi, with some deflections as high as 5.3 /mi and 
some as low as 4.0 /mi. For traveling wave deformation, one would expect equal displace¬ 
ments at all points around the stator. However, this is clearly not the case as we can see a 
standing wave component of deflection superimposed on the traveling wave in Figure 7-27. 
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This is very likely due to the different resonance characteristics of the two electrode patterns 
as illustrated in the damping measurements of the last section. A number of reasons are 
possible for the differences in quality factor between the two electrodes. First, the outer 
electrode covers a slightly larger area and at a larger distance from the center of the stator 
than the inner electrode. Second, the ceramic segments may not be poled uniformly, thus 
granting higher deflections to some areas of the stator than others. Third, the bond may 
not be perfectly uniform, or weakened in spots. The latter was definitely a possibility in 
the initial prototypes, but probably is not the case here, where a higher quality epoxy has 
been used. 

However, the sheer magnitude of the improvement in the displacement measured in 
these stators as compared to the initial prototypes is clearly the reason for the drastic 
improvement in output performance. We are seeing almost 10 times as much displacement 
as in the initial prototypes. Again, this is primarily due to improvements incorporated into 
the second batch of stators of a higher temperature epoxy bond and a full-area coverage 
electrode pattern. 

If we look at Figure 7-27(b), the plot of the phase of the turn-around displacement point 
with respect to the drive signal, we can see points of zero crossings at 4 locations around 
the stator. By symmetry, we can see that if we had been able to complete the experiment, 
that there would be two more zero crossings for a total of 6. Six zero crossings represent 
three wavelengths, so we can see that a three-wavelength traveling wave has been generated 
- precisely what the electrode pattern was designed for. 

This interferometric method then, gives a nice way of visualizing traveling-wave modes 
and verifying that the appropriate drive conditions are achieved. 

7.3.8 Efficiency Measurements and Drive Electronics 

Efficiencies have not yet been measured on the second generation 8 mm bulk PZT ultrasonic 
motors. While the dynamometer is adequate for measuring mechanical output power, the 
input power is not easily measured at this point because the currents are rather small and 
fairly noisy. A little bit of work has to be done on the drive electronics to clean up the 
signals in order to facilitate ease of measurement of the phase angles offsets between the 
two current and voltage signals, when calculating average input current. 

The drive electronics all reside on two small printed circuit boards. A Motorola 68HC11 


147 



processor controls a Qualcomm 2334 2-channel digital frequency synthesizer and also six 
LCD thumbwheel switches from American Control Technology. The LCD thumbwheel 
switches have embedded 4-bit microcontrollers and can be used as simple smart input and 
display devices for changing the frequency, phase and amplitude of the drive signals. The 
digital frequency synthesizer is controlled by the microprocessor via user input from the 
LCD thumbwheel switches and outputs two channels of 12-bit data at the appropriate 
frequency and relative phase. The two 12-bit ports are connected to AD7845 multiplying 
analog-to-digital converters to produce sin and cos drive signals for the motors. An AD558 
analog-to-digital converter, directed from the microprocessor, controls the amplitude of 
the output of the AD7845 converters. These two resulting channels of drive signals are 
passed through Apex PA85 high-voltage high-bandwidth (expensive) operational amplifiers 
to boost the output voltage to 60 Vpeak typically. 

A low-pass filter needs to be added between the AD7845 analog-to-digital converters 
and the Apex high-voltage stage in order to reduce noise on the output waveforms. The 
motors run just fine with these inputs, but the current signals are somewhat noisy and make 
it hard to measure efficiencies. This problem is next on the agenda. 

7.3.9 Coefficients of Friction 

The coefficients of static friction between the various combinations of rotor and stator 
materials used in the experiments here and are shown in Figure 7-28. Most combinations 
gave coefficients on the order of 0.3, but sapphire on aluminum had an extraordinarily high 
value of 0.5, which makes sense since like materials have very high coefficients of friction and 
sapphire is aluminum oxide. So aluminum oxide on oxidized aluminum would be expected 
to have large friction coefficients. Sapphire on stainless steel had the lowest coefficient of 
friction at 0.22. 

7.3.10 Matching to a Model 

Once we have measured stator displacements and coefficients of friction between various 
rotors and stators, we can compare measured output performance to predictions from the 
sequence of models developed earlier. 

First, we look at the case of line contact with Coulomb friction. Examining the data 
taken on stator Greyl, if we approximate the amplitude of displacement from Figure 7- 
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Static Coefficients of Friction 


01 Hardened Steel 

Sapphire 

Kapton-on-steel 

Kapton-on-sapphire 

304 Stainless Steel 

0.32 

0.22 

0.38 

0.38 

Aluminum 

0.35 

0.49 

0.35 

0.35 


Figure 7-28: The coefficients of static friction for the rotor and stator materials used in the designed 
experiment are shown here. These coefficients were determined by finding the minimum angle at 
which a rotor would slide down an inclined piece of stator material. 


27(a) as 4.5/mi, the wavelength as A = where the radius of contact is 3.25mm, or 
A = 6.8 mm, and the coefficient of static friction between 304 stainless steel and sapphire 
as 0.22, we can compute the predicted no-load speed and stall torque. 

The no-load speed for line contact is: 


2'Kujhw 0 

^rotor — ^stator,max — T 

A r 

where uJ ro t or is the rotor speed (equal to the no-load speed for all loads up until stall when 
slipping occurs), oj is the drive frequency in ppf, h is the lialf-lieiglit of the stator, A is the 
wavelength, r is the radius of contact and w 0 is the out-of-plane displacement of the stator. 

For stator Greyl, which has a tall base and short teeth on top of the PZT ceramic, h is 
0.67mm. The rotor speed then works out to: 

rad 

uj rotor = 219 -= 2093 rpm 

secs 

Similarly, if we note the normal forces in Figure 7-18 that achieved the largest and 
smallest stall torques, we see they ranged from iqy = 28 gf to iqy = 97gf. For a coefficient 
of friction of 0.22, predicted stall torques are: 


TStall = f-lF N r 


which would predict stall torques of 2.0gf-cm and 6.9gf-cm, respectively. Figure 7-29 illus¬ 
trates these two models of line contact superimposed on the data sets from Figure 7-18 for 
normal forces of iqv = 28 gf and iqy = 97 gf. 

It was found that if one goes through all the calculations of the Hertzian contact models 
for the actual data found from these experiments with the second generation bulk motors, 
that the contact condition is actually equivalent So that of line contact. 
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Figure 7-29: Data from the first speed-torque curve for stator Greyl with a plain sapphire rotor, 
where the normal force was the lightest, at Fjv=28gf is marked by the diamonds. Data from the 
seventh speed-torque curve where the normal force was the highest, at F ^=97 gf is marked by the 
squares. Predictions based on line contact assuming Coulomb friction are shown for both cases. 
Hertzian contact assumptions give the same prediction as line contact for the light normal forces 
only sustained here. 
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Recall that Hertzian contact theory states that the length of the contact region grows 


as the square root of the normal force per unit length: 


x 0 = - = 0.8 \/FnDCe 


For 304 Stainless Steel, the Young’s modulus is E = 28 X 10 6 psi (190x 10 9 and Poisson’s 
ratio is v = —0.21 [Oberg et al. 88] For sapphire, or alumina, the Young’s modulus is 
E = 54 X 10 6 psi (37xl0 9 ^r) and Poisson’s ratio is v = —0.20 [Coors 95]. These materials 
parameters give a Ce of: 


C'e 


1 ~ V 1 , 1 ~ 
E x + E 2 


where Ce works out to be Ce = 1.14 X 10 -11 

Taking thedisplacement of stator Greyl as 4.5 p, the equivalent radius of curvature 
becomes pj~ = 0.26 m. Modeling the stator as a cylinder of diameter, D = 2 pj~ gives 
D = 0.52 m. Plugging in these materials parameters we find the length of the region of 
contact is 32 pm or ^ for a normal force of Fn = 28 gf of 61 pm and ^ for a normal force 
of Fjv = 97gf. Looking back at Figure 6-24, the speed-torque curves for Hertzian contact, 
we see that this condition approaches the case of line contact. Indeed, following through 
and calculating the speed-torque relations from: 


where 


7~rotor 


nF N r 

4>{x 0 ) 


[ 2 4>{x r ) - 4>{x 0 )\ 


cf) ( x ) = sin ( kx ) — kxcos ( kx 0 ) 


and 


^rotor — ^ stator, maxCOS 


(kx r ) 


and 

1 —1 / ^rotor \ 

x r = —cos - 

yX? stator,max / 

gives essentially the same lines as in Figure 7-29 for line contact. Note that the no-load 
speeds in reality are lower than the predicted values and the stall torques are higher. The 
coefficient of friction is the scale factor for the stall torques and is not known with complete 
accuracy. The no-load speeds are lower than those predicted for line contact, both because 
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of the assumption of a rigid, unsquashable stator in these closed-form models and because 
of the decrease in the amplitude of deflection as power is taken out of the system. This 
shows that it is important to model the squashing of the stator in addition to area contact. 
Appendix B illustrates variational techniques for producing models that relax the rigidity 
assumption on the stator. 

7.3.11 Surface Profilometry Measurements 

The second generation set of 8 mm motors saw significant improvement in performance 
over the initial prototypes whose speed, torque and power characteristics were plotted in 
Figure 7-5. With the initial prototypes, we found that glass rotors spun well, while brass 
rotors turned on a lathe did not work at all. This was understandable when later we found 
that the mean roughness of the glass rotor was 300 A while that of the brass rotor was 
7000 A. In addition, the amplitude of displacement on those first stators was under 1 fim 
meaning that the out-of-plane displacement of the stator was not much greater than the 
height of the asperities on the rotor. 

In the second generation of motors, we took great pains to achieve large out-of-plane 
deflections and indeed we have seen a factor of 5 to 10 improvement. We also paid more 
attention to the surface properties of the rotor and stator materials. We not only chose 
harder materials, but we polished each surface before testing in the dynamometer. The 
stators and rotors were lapped flat and then polished with 9 /im-, 5 /im-, 3 /im- and finally 
1 n m-grit lapping paper. 

After polishing, surface properties were examined with a Dektak prohlometer. Figure 7- 
30 illustrates surface profiles for the polished steel and sapphire rotors, respectively. The 
avaxis shows the length of material over which the stylus was run, measured in microns, 
while the y -axis denotes the asperity height, measured in angstroms or kilo-angstroms, as 
noted. The steel rotor, in Figure 7-30(a) had a mean roughness of 115 A, polished almost 
as smoothly as the sapphire surface which displayed a mean roughness of 110 A. 

The kapton-coated rotors were not so smooth, obviously. Figure 7-31 displays a mean 
roughness of 768Afor the kapton-coated steel rotor and a mean roughness of 411 Afor the 
kapton-coated sapphire rotor. 

It is important to examine not only the rotors, but also the surface of the stator teeth. 
Figure 7-32 depicts the surface profile across a tooth of stator Greyl just after polishing. 
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Figure 7-31: (a) The kapton-coated steel rotor measured after polishing displays an arithmetic 
mean roughness of 768 A. (b) The kapton-coated sapphire rotor has a mean roughness of 411 A. 
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Figure 7-32: This is the stator Grey 1 tooth surface.; The arithmetic mean roughness is 631 A. 

The measured mean roughness here is 631 A. These profilometer traces in conjuction with 
significant improvements in torque output, as compared with initial prototypes, point out 
the need for out-of-plane displacements to be much larger than surface asperities. 

Wear Tests 

Surface profile measurements were also performed on the rotors after running in the dy¬ 
namometer in order to study wear phenomenon. After the rotors are run on the stators 
for 30 minutes or so, the rotors and stators are removed from the dynamometer and tlieir 
surfaces inspected. Typically, a discolored ring is seen on the rotor at the radius of contact 
and the stator teeth have a radial section which is highly polished while the remaining areas 
are rather dull in color. 

Figure 7-33 illustrates wear profiles on plain steel and sapphire rotors after having run 
on stator Greenl, which was stainless steel. The stator makes a circular gouge mark on the 
rotor and the profilometer is run across the gouge in the radial direction from the center 
of the rotor out towards the edge. Figure 7-33(a) shows that for the 01 hardened tool steel 
rotor running on a stainless steel stator, material is deposited on the rotor to a thickness 
of 1942A. Similarly, Figure 7-33(b) shows that the sapphire rotor, after running on the 
stainless steel stator, also receives a deposit of 901 A. Thus the rotors are harder than the 
stator and the stator does not gouge the rotor, as was the case in the initial prototypes 
where the stators ate into the first glass rotors. Here, the rotors actually tend to polish the 
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Figure 7-33: (a) The steel rotor, measured across the discolored circular mark, shows a peak of 
1942A. This rotor was last run on stator Greenl (stainless steel), (b) The sapphire rotor measured 
across the discoloration, acquires a peak of material 901 A high. This rotor was also last run on 
stator Greenl. 
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Figure 7-34: (a) The kapton-coated steel rotor measured across the circular mark shows an inden¬ 
tation into the kapton of nearly 2pm . This stator was run on stator Greyl (stainless steel), (b) 
The kapton-coated sapphire rotor measured across the circular mark displays a valley of 121 A which 
is followed by a peak of 1956 A. 
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stators. 


Figure 7-34 portrays surface profiles of the kapton-coated rotors after having been run 
on stator Greyl, which was also stainless steel. Figure 7-34(a) shows what happens to the 
kapton-coated steel rotor after 30 minutes or so of running. The stainless steel stator eats 
into the kapton to a depth of nearly 2 fj,m. This is as would be expected since the polymer 
material is rather soft. Figure 7-34(b) shows a similar result for the kapton-coated sapphire 
rotor. The difference here, however, is that a valley of 121A is followed by a peak of 1956 A, 
so the kapton material is smushed to the side as the teeth dig into the polymer. 

7.4 Comparing Motor Figures of Merit 

7.4.1 Stall Torque and Output Power Densities 

By dividing the measured maximum stall torques and output powers by the weight or volume 
of the stator-rotor assembly, we can examine figures of merit for our motors. [Hollerbach, 
Hunter and Ballantyne 91] argue that the most useful figure of merit for comparing disparate 
motor technologies is stall torque density (^j 1 ) as power density figures really depend on 
the associated power electronics. 

Figure 7-35 lists the stall torque densities resulting from our 16-trial experiment. The 
maximum stall torque density of 2.9 is produced by stator Greenl with a sapphire rotor. 
Note that the motors here do not have bearings, but are held together in the dynamometer. 
True stall torque density numbers will be smaller after bearings and mounts are added in 
the future. We have put off working on bearings and mounts until the motors were better 
understood and useful torques could be shown. This has been achieved and in the near 
future we will mount a simple bearing to one of these motors. 

Figure 7-35 lists the corresponding output power densities for the second generation 
8 mm bulk motors. Stator Fucial with a sapphire rotor produced the highest power density 
of 108 g. 

7.4.2 Comparisons to Other Actuator Technologies 

It is interesting to compare piezoelectric ultrasonic motors and see how other actuator 
technologies measure up. Figure 7-37 [Hollerbach, Hunter and Ballantyne 91] shows how a 
wide variety of technologies, such as hydraulic systems, pneumatic actuators, shape memory 
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Design of Experiments Results: Stall Torque Densities 


Stator 

Max Stall Torque 
(gf-cm) 

Mass (g) 
St+Rot=Tot 

Volume 

(mm A 3) 

Max Stall Torque 
Density 
(Nm/kg) 

Max Stall Torque 
Density 
(gf-cm/mm A 3) 

i 

Violet Exp 2 

4.7 

.28+.13=.41 

120 

1.2 

0.039 

2 

Violet Exp 4 

4.8 

.28+.37=.65 

120 

0.74 

0.040 

3 

Fucia Exp 2 

5.5 

43+43=.26 

120 

2.1 

0.046 

4 

Fucia Exp 4 

4.5 

.13+.37=.50 

120 

0.90 

0.038 

5 

Tan Exp 2 

2.7 

.12+. 12=.24 

87 

1.1 

0.031 

6 

Tan Exp 4 

3.4 

.12+.36-48 

87 

0.71 

0.039 

7 

Blue Exp 2 

4.7 

.27+42=.39 

87 

1.2 

0.054 

8 

Blue Exp 4 

4.6 

.27+.36-60 

87 

0.77 

0.053 

9 

Green Exp 2 

10.0 

.22+42=.34 

68 

2.9 

0.15 

10 

Green Exp 4 

8.4 

.22+.36=.58 

68 

1.5 

0.12 

11 

Yellow Exp 2 

5.3 

41+42=.23 

68 

2.3 

0.078 

12 

Yellow Exp 4 

5.6 

.11+.36-47 

68 

1.2 

0.082 

13 

Orange Exp 2 

4.8 

40+43=.23 

61 

2.1 

0.079 

14 

Orange Exp 4 

4.6 

.10+.37=.47 

61 

1.0 

0.075 

15 

Grey Exp 2 

4.3 

.20+43=.33 

61 

1.3 

0.070 

16 

Grey Exp 4 

5.3 

.20+.37=.57 

61 

0.93 

0.087 


Figure 7-35: Stall torque densities resulting from the 16-trial experiment of 8 mm bulk PZT 
ultrasonic motors. Stator Greenl (short teeth, thick base, 36 teeth, stainless steel) with a plain 
sapphire rotor produced 2.9-=^- stall torque density. 


Power Density 

Stator 

Liner 

Rmat 

Mass (g) 
St+Rot=Tot 

Volume 

(mm A 3) 

Peak Power 
fmW) 

Peak Power 
Density 
(W/kg) 

Peak Power 
Density 
(mW/mm A 3) 

Fucia 

- 

+ 

43+42=.25 

120 

27 

108 

0.23 

Yellow Exp 2 

- 

+ 

41+42=.23 

68 

12 

52 

0.18 

Grey 

- 

+ 

.20+42=.32 

61 

16 

50 

0.26 


Figure 7-36: Output power densities for the second generation motors. Power densities as high 
as 108 have been achieved with stator Fucial (tall teeth, thick base, 24 teeth, aluminum) with a 
plain sapphire rotor. 
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Actuator 

Power 

Dens. 

(W/kg) 

Stall 

Torque 

Dens. 

(Nni/kg) 

McGill/MIT EM Motor 

15 

200 

Sarcos Dextrous Arm electro- 
hydraulic rotary actuator 

120 

600 

Utah/MIT Dextrous Hand 
electropneumatic servovalve 

20 

200 

NiTi SMA [Hirose 89] 

1 

6 

PVA-PAA polymeric actuator 
[Caldwell 90] 

17 

6 

Burleigh Instruments inchworm 
piezoelectric motor 

3 

0.1 

Magnetoelastic (magnetostrictive) 
wave motor [Kiesewetter 88] 

500 

5 

Human biceps muscle 

20 

50 

MIT 8 mm Fucia/Green 

108 

2.9 


Figure 7-37: Power density and stall torque density can be used as figures-of-nrerit. to compare 
actuator technologies. 

alloys, piezoelectric inch worms, magnetostrictive materials and human muscle compare in 
terms of both stall torque density and output power density. We have inserted our recent 
results into their table where the number for power density refers to the stator Fucial 
running on a plain sapphire rotor and the number for stall torque density refers to the 
stator Greenl also running with a plain sapphire rotor. 

Note that the piezoelectric ultrasonic motors vastly outperform the piezoelectric inch- 
worm drives. This is due the fact that resonance is used to advantage in the ultrasonic 
motors, whereas the inchworm drives are quasistatic. 

While the stall torque densities for the first three actuators in the table are significantly 
higher than those for our ultrasonic motors, it is useful to note that these are all fairly large 
actuators. The McGill/MIT electromagnetic motor is also water-cooled. If we compare our 
8 mm ultrasonic motors to small DC motors that would typically be the direct competition, 
our ultrasonic motors compare exceedingly well. 

Figure 7-38 gives figures of merit comparing ultrasonic motors reported in the literature 
with small electromagnetic DC motors typically used in mobile robot applications. The 
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Actuator 

Stall 

Torque 

(kgf-cm) 

No-load 

Speed 

(rpm) 

Power 

Density 

(W/kg) 

Stall 

Torque 

Density 

(Nm/kg) 

Peak 

Eff. 

[Micro Mo] 1319E003 

0.034 

13,500 

106 

0.29 

71% 

[Mabuchi] FK-280-2865 

0.16 

14,500 

NA 

0.42 

53% 

[Namiki 92] 7CL-1701 

0.0016 

9,100 

13 

0.06 

50% 

Aeroflex 

0.10 

4,000 

NA 

0.04 

20% 

Maxon 

0.13 

5,200 

NA 

1.13 

70% 

Astro 

0.76 

11,500 

NA 

0.21 

20% 

[Kumada 85] 

13.5 

120 

80 

8.8 

80% 

[Shinsei 89] USR-60 

6.0 

125 

23 

3.4 

35% 

[Panasonic 87] USM-40 

0.8 

800 

NA 

1.1 

40% 

[Hosoe 89] EF35-135 

1.2 

80 

NA 

2.3 

40% 

MIT 8 mm Fucia 

0.0055 

1750 

108 

2.1 

NA 

MIT 8mm Green 

0.01 

870 

NA 

2.9 

NA 


Figure 7-38: This representative sampling of electromagnetic DC motors (without gears), com¬ 
mercially available Japanese ultrasonic motors and our 8 mm ultrasonic motors are presented to give 
a feel of relative measures of stall torque density and efficiency. Stators Fucia and Green are run 
with plain sapphire rotors. 

first six motors listed in Figure 7-38 are electromagnetic motors, the next four are larger 
Japanese ultrasonic motors, on the order of 40 mm to 60 mm in diameter. 

The smallest commercially available electromagnetic DC motor is the Namiki 7C'L- 
1701. It is 7 mm in diameter and 18 mm long, roughly 5 times the volume of our 8 mm 
ultrasonic motors, yet stator Fucia running with a sapphire rotor lias 35 times the stall 
torque density and 8 times the power density of the Namiki motor. Stator Green with a 
sapphire rotor compares even more favorably to the Namiki in stall torque density. Its stall 
torque density is 48 times larger. Most importantly, however, is that the typical operating 
speeds of ultrasonic motors are much lower, granting a significantly better impedance match 
for nearly all loads. 
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Chapter 8 


Laser Etching of Thin-Film PZT 
Motors 


The original impetus for pursuing ultrasonic micromotors was the fact that thin-hlm PZT 
could now be deposited directly onto silicon and the resulting increase in breakdown strength 
provided extremely high energy densities. By depositing PZT in a sol-gel process onto sili¬ 
con, we do away with the need for an epoxy bond layer and we open the door to possibilities 
for integration with electronics and microsensors. Because of the high energy densities and 
because these ultrasonic motors inherently run at high torques and low speeds, we have a 
promising new technology for creating microfabricated motors which couple useful power 
to a load. 

To this point, we have shown process feasibility in fabricating thin-hlm actuators and we 
have modeled and tested larger motors from bulk PZT with very good results. Now we would 
like to join the two efforts and return to microfabricated ultrasonic motors. However, we are 
not necessarily interested in making these motors micro-sized, nor are we so constrained, as 
demonstrated by the larger motors we have built. While silicon electrostatic micromotors 
are constrained to be under a certain maximum radial dimension due to residual stresses in 
the rotor vanes causing them to touch down onto the surface below, friction-based ultrasonic 
motors can be made arbitrarily large in the plane of the stator. The constraint here is the 
ability to fabricate wide-area coverage films that are pinhole-free. 

To get motion out of the plain of the wafer, and to achieve true waveguide structures, 
it is desirable to cut the stator free from the wafer. Traditional machining is not amenable 
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Figure 8-1: The two rings at the bottom right are 8 mm composite rings consisting of layers of 
silicon, silicon oxide, titanium, platinum, sol-gel PZT and chrome-gold. The stator at left is 5 mm 
in diameter. The six-pin integrated circuit at top is a serial charge coupled imaging device, to give 
scale to the picture and to illustrate our desire to build actuators at the same size as our sensors. 

as silicon is fairly brittle, and etching processes create sidewalls at the angle of the lattice 
lines. Wliat would be best is a non-contact cutting tool. 

Figure 8-1 shows stators that we have machined using a laser-etching process. The laser 
etching was performed at a commercial job shop, Laser Services, in Westford, MA. The 
laser ablation process is used to cut rings out of wafers that have been prepped with all 
the layers of the necessary films: oxide, titanium, platinum, sol-gel PZT and chrome-gold. 
The laser process is also used to etch the electrode pattern in the chrome-gold and also the 
tooth structure on the backside to reduce stiffness. We also start with thinned wafers, and 
in this manner avoid the need for fragile membranes and silicon-rich nitride films. 

To connect to the titanium-platinum ground electrode, a chemical-assist laser-etching 
and deposition process was used. This step was done with the assistance of Dr. Daniel 
Ehrlich at MIT Lincoln Laboratory. A platinum gas is used with the laser to deposit a 
plug in a laser-punched via. This effectively brings up a contact from the bottom electrode 
to a top chrome-gold pad. In this manner, all the steps of stator microfabrication are 
accomplished without any masks or wet-etcliing. 


Chapter 9 


Conclusion 


The technology of micromotors has taken a new turn. Previously, micromotors were fabri¬ 
cated from predominantly silicon integrated-circuit processes. With the incorporation of an 
electromechanically active material, PZT, into a standard silicon micromachining process, 
we have been able to produce piezoelectric ultrasonic micromotors displaying significant 
improvements over previous incarnations of micromotors. 

Ultrasonic motors are new and not as well understood as conventional electromagnetic 
or electrostatic motors. A number of models with closed-form solutions for production of 
torque have been proposed here, where simplifying assumptions on the contact mechan¬ 
ics have been sequentially relaxed. This analysis has granted insight into the phenomena 
underlying the frictional coupling mechanism and provided useful upper bounds on perfor¬ 
mance. 

In addition to theoretical results, an experimental component of this research has yielded 
one of the most complete data sets on bulk piezoelectric ultrasonic motors in the open 
literature. Providing this data to other researchers in the held will not only contribute to 
elevating our understanding of these devices, but also hasten growth of this new technology. 

We have found that our second generation of bulk motors shows drastically improved 
performance over our initial prototypes, due to improvements in electrode design, bonding 
techniques and attention to detail in surface preparation. These new motors work primarily 
in the region of line contact. The experiments have shown that ultrasonic motors can deliver 
five times the stall torque of their electromagnetic brethren in | the size, with comparable 
power densities. Power densities of 108 and stall torque densities of 2.9-^^ have been 
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achieved. In addition, the impedance match for these motors to most real-world loads is 
significantly better than for electromagnetic motors. Peak power operating points occur at 
speeds on the order of a few hundreds of rpm rather than thousands of rpm. Careful design 
of the stator geometry can result in operating points tailored to a load such that the need 
for geardown can be circumvented. 

In the future, we plan to return to thin-hlm PZT-on-silicon technology and incorpo¬ 
rate these high-torque low-speed motors into complete microelectromechanical systems for 
microrobotic applications. A new laser-based etching process has been developed in this 
regard which allows for microfabrication of ultrasonic motors without fragile membranes 
and requires no masks and no wet-etching. This type of process is more amenable to the 
designer than batch fabrication processes and such a manufacturing capability will further 
aid the growth of the held by compressing the design, fabrication and testing cycles. 
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Appendix A 


Stator Motion 


[Inaba et al. 87] described elliptical motion on the surface of a beam due to 
wave as illustrated in Figure 3-4, where one solution to the beam equation is 
two standing wave oscillations offset by ^ in space and time: 

w(x,t ) = w 0 (sin kx smut + cos kx cos ut) 

which is simply the equation of a wave traveling along the neutral axis, 

w(x,t ) = w 0 cos (kx — ut ) 

For a beam of half-thickness h , however, the displacement of a point on the surface moving 
from P to P 0 has a horizontal displacement ( a and a vertical displacement £ a : 


a travehng 
the sum of 


£ a = —(h— h cos 0) + w 0 sin kx sin ujt 


£ a = — h sin 0 


To solve for 0 , we note that the slope of the normal to the neutral axis is the negative 
reciprocal of the slope of the tangent to the neutral axis at that point. 


tan(| + 6) 


1 

dw 

dx 


l _ l 
tan 6 ^ 
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Flexing Beam 


d i w(x,t ) pAd 2 w(x,t ) 

dx 4 + E/ df- = 

Traveling Wave Solution w(x,t ) = W 0 COs(kx — Lot ) 


ut = ^ 



Figure A-l: (a,) For a beam in flexure, a traveling wave is one solution to the beam equation. A 
vibration of amplitude w 0 will cause a point on the surface to displace from P to P Q . 


For small deformations, 9 is small and 


sin 9^9 


cos 6=1 


Therefore | or we can say, 

dx 


n ( l' W i , 

9 = —— = w Q k cos kx smwf 
dx 


= w 0 sin kx sin ujt 


Ca = —h9 = — hw 0 k cos kx sin u>t 


If a, second standing wave is induced which is spatially and temporally 90 degrees out of 


phase from the first, the displacements will be: 


7T . 7T 

^ = w 0 sinj/or + —) sin(cjf + —) = w a cos kx cos cot 


7T . 7T 

£& = —hkw 0 cos( kx + —) sin(ud + —) = —hk-w 0 sin kx cos 
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Superimposing these two waves gives the displacements due to traveling waves: 

£ = £ a + £& = w 0 (sin kx sin ut + cos kx cos ut) 

£ = £ a + £& = hkw 0 (cos kx sin ut — sin kx cos ut) 

From the trigonometric identities, 

cos(a — (3) = cos a cos (3 + sin ci sin (3 

sin(ci — (3) = sin ci cos (3 — cos a sin (3 
we see that traveling waves result, with: 


£ = w 0 cos (kx — ut) 


£ = hkw 0 sin (kx — ut) = 


2tt hw 0 
A 


sin(A;a: — ut) 


The relationship between the transverse and longitudinal displacements of a point on 
the surface can be seen to be elliptical if we square and add £ and £: 


1 = sin 2 (A;a: — ut) + cos 2 (kx — ut) = 


c 2 


£ 2 


' 2TThw„ \9. 1n 2 

- A ) 0 
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Appendix B 


Variational Methods in Modeling 
Active Structures 


Closed-form solutions cannot be formulated for the majority of structural systems which 
include active, or electromechanically coupled materials. For systems which exhibit complex 
geometry, boundary conditions and load configurations, approximation techniques based on 
energy methods and variational principles are employed. This appendix gives background 
for modeling active structures using variational methods such as Rayleigh-Ritz assumed¬ 
mode solution techniques. 

The aerospace community has a long history of using sophisticated analysis tools for 
modeling complex structures, such as composite load-bearing members which must be 
lightweight yet strong and subject to a wide range of forcing conditions, from mechani¬ 
cal to aerodynamic in origin. Active control of such structures is often employed to increase 
rigidity, dampen vibration or steer precision optical components. Recently, incorporation of 
active materials such as piezoelectrics, electrostrictors, magnetostrictors and shape-memory 
alloys into controlled structures has been a prolific area of research as a route to further 
reduce system size and weight by alleviating the need for discrete actuation mechanisms 
and linkages. 

Much of this recent work has helped to bridge the gap between the materials science and 
aerospace communities by developing engineering models of complex electromechanically 
coupled materials which can be used in the analysis tools employed for the design of active 
structures [Crawley and De Luis 87], [Crawley and Anderson 90] and [Rogers, Liang and 
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Jia 89]. While the aerospace devices being modeled are often large-scale structures such as 
space-station truss elements, heavy-lift launch vehicles, satellite optical systems or twistable 
helicopter rotor blades, many of the analysis techniques can be extended to the new held 
of microelectromechanical systems where the active materials may be thin films deposited 
directly onto silicon, the interest here. 

The remainder of this appendix gives background, fills in details, and equates notation 
from three papers which use these Rayleigh-Ritz methods to model piezoelectric structural 
elements such as static and dynamic beams and plates, [Hagood, Chung and Von Flutow 90] 
and [Crawley and Lazarus 91], and more recently, ultrasonic motors [Hagood and McFarland 
95]. At the end of this section we show a simulation of the latter model for the 8 mm 
ultrasonic motors characterized in this thesis work. It is hoped that this exposition will 
be helpful to the microelectromechanical systems community for modeling active materials 
incorporated into microstructures. 

B.l Hamilton’s Principle 

Hamilton’s Principle is the starting point from which the dynamic equations of motion are 
formed in [Hagood, Chung and Von Flutow 90]. The use of Hamilton’s Principle in modeling 
piezoelectric plates is also discussed in [Tiersten 67] and [Tiersten 69]. Hamilton’s Principle 
states that a system will evolve along a trajectory in such a manner as to minimize the 
action, where the action is defined as an integral over time of a variational indicator of the 
system. The deflections of the static system in [Crawley and Lazarus 91] are derived in a 
similar manner, but from the Principle of Minimum Potential Energy, which simply states 
that the variation of the variational indicator is zero. In a dynamic system, the action is 
minimized if the variation of the action is equal to zero: 

*2 

6 j ( V.I .) dt = 0 

ti 

Figure B-l illustrates the variation of a function. [Lanczos 70] describes the calculus of 
variations and its relationship to mechanical systems in which infinitesimal displacements 
and virtual work processes are applied. We assume here linear variations of a function f{x) 
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Figure B-l: A variation of a function is related to its dependent variable by 6f = §^-8x. 
by a function (f>(x) such that the variation of f(x) is: 

Sf (x) = €<j)(x) 

where e is infinitesimally small and the variation at the endpoints of the path are zero: 


s f( x ) x =a = S f( x ) x =b = 0 


A variation of a function is related to its dependent variables by the chain rule: 


,., df df df 

8 f — T. - 8x\ + T. - h x 2 + • • • + T. - 

UX\ UX 2 Ox n 


8x r 


The technique of applying an arbitrary variation in position to determine the virtual change 
in potential energy (a function of position) is the step which allows the equations of motion 
to be extracted from an accounting of energy. 

The particular form of the variational indicator chosen by [Hagood, Chung and Von 
Flutow 90] for inclusion into Hamilton’s Principle specifies stress and electric held as the 
independent variables. Stress and electric held are typically the control parameters available 
in driving active structures in actuator configurations. The interaction of these mechanical 
and electrical held variables is derived from thermodynamics. 

The first law of thermodynamics states that heat added or work done on a system 
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increases the total internal stored energy density of the system: 


dU tot = dQ + dW 

where U tot is the total internal energy density, Q is the heat per volume added and W is the 
work per volume acting on the system. Assuming a linear electro-elastic-thermal system, 
the change in total internal energy density due to infinitesimal changes in strain and electric 
displacement in the presence of uniform stress and electric held, becomes: 

dU tot = Ode7 + TijdSij + E k dD k 

where 0 is temperature, a is entropy, Tij is a component of the stress tensor as defined in 
Figure B-2, Sij is a component of the strain tensor as illustrated in Figure B-3, E k is a 
component of the electric held vector and D k is a component of the electric displacement 
vector. 

Legendre transformations can be applied to the energy expression to switch the inde¬ 
pendent variables from (S,D) to (S,E) [Nye 85]. The Electric Gibb’s Free Energy expression 
is one such transformation and is dehned as: 

G = U tot -da - E k D k 

where 

dG = dU tot — d(da) — d(E k D k ) 

and 

dG(0, S, E) = -add + T^dS^ - D k dE k 


By noting that: 


_ dG dG ia dG Jri 
dG (d, S, E) — + ~dE C ^ 

we hnd the relationships between the dependent variables and G: 


„ _ f rJG\ 

G ~ I, d6 jr 


T - — 

- 1 - 1. 1 - 


dG 


D k = — 


dG 


' t,e lJ \ dS u)e,E k \ dE k)e,T 

where the subscripts signify variables that are held constant. Expanding the change in 
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Figure B-2: A force per unit area applied to any given face of a cube (a traction vector) has 
stress components in three directions. The stress tensor is a 3 x 3 matrix of components X)j and is 
symmetric. The six independent components can then be written in reduced matrix notation as a 
1x6 vector. 

dependent variables as: 


dT(0,S,E) 

dD(0,S,E) 


dT DT DT lri 
— d0 + — dS + — dE 


dD 

~de 


dO 


dD 

~ds 


dS 


dD 

Je 


dE 


and assuming isothermal conditions, gives the constitutive relation coefficients for the lin¬ 
early electromeclianically coupled piezoelectric system: 


E,e = (mA = ( d 2 G \ 
i3kl \dSkiJ E ,e ) h . 

T ,e _ (dDk\ ( d 2 G \ 

£i3k \ dS Ej Tt0 \dS t 3 dE k ) Tg 

T ,e _ (dD k \ ( d 2 G \ 

€kl KdEj^g \dEkdEj T e 

where cf^ is the stiffness tensor at constant electric field and temperature, e Jjk is the 
piezoelectric stress tensor at constant stress and temperature and is the dielectric 
permittivity tensor at constant stress and temperature. 
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Because the stress tensor is symmetric: 


- 1 

t 12 

1 

CO 


1 

T 6 

t 5 - 

T 2 1 

T22 

CO 

= 

T 6 

t 2 

t 4 

_ 1 

T32 

1 

CO 


-T 5 

t 4 

t 3 . 


we can write Tin a reduced matrix form: 


-Tr 


'Tn " 


1 - 

b 

1_ 

t 2 


T22 


(Jy 

t 3 


CO 


°z 

t 4 


CO 


T yz 

t 5 


T 3 i 


T~zx 

-Te- 


-T21- 


L T xy - 1 


where < 7 , r is traditional notation in uncoupled systems for normal and shear stresses and 
is used in [Crawley and Lazarus 91]. The strain tensor is also symmetric and can similarly 
be written in reduced matrix notation: 



r 9 

u xx 

2 S xy 

2 S xz - 


-Si 

s 6 

s 5 - 

s = 

2 S xy 

°yy 

9 8 

^Uy Z 

= 

s 6 

S 2 

Si 


-2 S xz 

9 8 

^Uy Z 

Szz - 


-S 5 

Si 

-1 

^70 


The strain-displacement relations are defined in cartesian coordinates in Figure B-3 and 


given as: 



-sr 


' Sn - 


^X 


r -| 

dx 


s 2 


S22 


€ y 


dv 

dy 

s = 

s 3 


s 33 


€ z 


dw 

dz 

Si 


2S 23 


lyz 


dw _i_ dv 
dy ' dz 


s 3 


2S 3 i 


S zx 


dw 1 du 

dx ' dz 


-S 6 - 


-2 S21- 


-Sxy - 


du 1 dv 
dy ' dx 


where e, 7 is engineering notation for normal and shear strains, also used in [Crawley and 


Lazarus 91]. Assuming isothermal conditions, the piezoelectric constitutive relations can 
now be written in compact form [IEEE 87]: 
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(b) 

Figure B-3: (a) Strain produced by positive normal stresses, (b) Strain produced by positive shear 
stresses. 
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where the stiffness tensor at constant electric held, c E , becomes a 6 X 6 matrix in reduced 
form, the piezoelectric stress tensor transposed, e t , becomes a 3 X 6 matrix and the permit¬ 
tivity tensor at constant stress, e T , remains 3x3: 


-1 

_1 


1- 

to 

1_ 

T 2 


S 2 

t 3 

= c s 

CO 

t 4 


s 4 

t 5 


^5 

1- 

1_ 


-1 

to 
_1 


r^ii 


et 


E 2 

-E 3 _ 


r-s'ii 


■ d r 
d 2 

-D 3 _ 


S 2 

^3 

S 4 

^5 

-S 6 - 


-E x - 

E 2 

- E 3 _ 


For PZT, the material coefficient matrices take the following form: 


C PZT 


,E 

11 

r E 

c 12 

r E 

C 13 

0 

0 

0 

,E 

12 

r E 

c ll 

r E 

C 13 

0 

0 

0 

,E 

13 

r E 

C 13 

r E 

c 33 

0 

0 

0 

1 

0 

0 

r E 

c 44 

0 

0 

1 

0 

0 

0 

r E 

c 5-5 


1 

0 

0 

0 

0 

2( 



Ki 

0 

0 

1 


e PZT — 


( T 

e ll 


,E 

-'ll 


Lo o 


T 

e 33 J 
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epzT t 


'0 0 ei 3 ' 

0 0 ei3 

0 0 e 33 

0 ei5 0 

eis 0 0 

.0 0 0 . 


Alternate Legendre transformations of the energy expression can be used to state the 
constitutive relations in terms of other combinations of independent variables, such as (T,E) 
[Jaffe, Cook and Jaffe 71]: 


'S' 



+ dt' 


T 

D. 


. d 

e 5 . 


E. 


where d is the piezoelectric strain matrix, s E is the compliance matrix at constant electric 
held and e s is the dielectric permittivity matrix at constant strain. The compliance matrix 
is related to the stiffness matrix by: 



The piezoelectric strain matrix is related to the piezoelectric stress matrix by: 



and the constant strain and constant stress permittivity matrices are related by: 



With the constitutive equations defined, we return to the formulation of the variational 
indicator used in Hamilton’s Principle. In matrix notation and under isothermal conditions, 
the change in the Electric Gibb’s Free Energy, where G' is in units of energy rather than 
energy density, is: 

dG' (S, E) = J (T • dS — D • dE) dV 
v 

We take: 
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dU = f T ■ dSdV dW e i ec = / D • dEdV 

v v 

Integrating, we find expressions for U and W e [ ec : 

U = f |S t T dV W P j ec = f |D t E dV 

v v 

We can then set up an energy expression: 


U I V e { ec — IFmec/i 

where U is the stored mechanical energy of the system, W mec h is the work of mechanical 
origin done on the system and W e i ec is the work of electrical origin done on the system: 

The work of mechanical origin, W mec h , includes components involving accelerations. 
Rearranging so as to separate out this kinetic energy term, Tki n , we can define: 

fkmec/j, — Tfc i n T W 

where W is the work of mechanical origin not involving accelerations. Substituting the 
resulting variational principle back into Hamilton’s Equation, we find: 

*2 

s j (Tkin -U + Weiec + W) = 0 

ti 

which is the form of Hamilton’s Principle used in Equation (1) of [Hagood, Chung and Von 
Flutow 90] 

B.2 Laminated Plate Dynamics 

[Hagood, Chung and Von Flutow 90] continues by considering the dynamic excitation of 
a beam structure. Here, we show the derivation of the dynamic equations of motion of 
a laminated plate from Hamilton’s Principle using Rayleigh-Ritz assumed modes and the 
consistent plate hypothesis of [Crawley and Lazarus 91], and correlate the notation therein 
to that of [Hagood, Chung and Von Flutow 90]. We also describe an implementation of 
the model in both rectangular and cylindrical coordinates [Hagood and McFarland 95] and 
show deformation predictions of microstructures. 
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Hamilton’s Principle can be used with specific geometries, materials and loads to find 
the equations of motion of a structure. The constitutive equations and external work terms 
are substituted into the energy expressions and assumed forms of the deformation shapes 
can be used to approximate the solutions. 

For the thin plate shown in Figure B-4(a), plane-stress assumptions can be made which 
leave only the Ti, T 2 and Tq terms in the constitutive equations. By setting T 3 , T 4 and T 5 
to zero in: 


'S' 



d t 


T 

D. 


. d 

e s _ 


E. 


and converting the resulting reduced s E , d and e s matrices back into the form required in: 


T 


-c E 

— e t" 


'S' 

D. 


e 

e T . 


E. 


the plane-stress materials matrices can be found. The materials matrices for various layers 
of materials in the plate may have to be further modified if they are orthotropic and skewed 
with respect to the laminate axes as shown in Figure B-5. Appropriate tensor transfor¬ 
mations must be applied to find the equivalent coefficients in the plate coordinate system. 
Such a transformation can make the stiffness matrix fully populated, creating the effect of 
an anisotropic material, the characteristic which [Crawley and Lazarus 91] use to achieve 
twisting of their cantilevered plate. 

The strain-displacement relations for a plate can also be simplified by making approxi¬ 
mations as illustrated in Figure B-4(b) where the Bernoulli-Euler assumption of plane sec¬ 
tions perpendicular to the midline remaining perpendicular after bending, is shown. The 
displacements in the x,y and z directions become: 


u(x,y,z) = u 0 - z 

v(x,y,z) = v 0 - z 
w(x,y,z) = w 0 


dw 0 

dx 

dw 0 

dy 


where u 0 , v 0 and w 0 are the centerline displacements in the x, y and z directions, respec¬ 
tively. Plugging these deformation expressions into the strain-displacement relations, the 
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(b) 


Figure B-4: (a) For a thin plate, plane stress assumptions are used, where all but the in-plane 
stresses are assumed negligible, (b) The strain-displacement equations are found by observing a 
section of the plajfi and assuming that perpendicular sections remain perpendicular to the midline 
after bending. 
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X,1 X, 1 

(a) (b) 


Figure B-5: (a) Definition of direction cosines, for rotating material tensors, (b) An or¬ 

thotropic material, such as a fiber composite, may be oriented with its principal axes skewed from i lid 
plate axes. Rotating into plate coordinates creates a fully-populated stiffness matrix, an anisotropic 
characteristic. 


strain vector results: 


' Si' 

S-2 

-S 6 - 


du n _ - C) 2 Wo 

dx dx 2 

dv„ _ - 9 2 W n 
dy ~ dy 2 

du n I dv n _ O- 9 2 Wn 

dy ‘ 6 >x ' dxdy 


where centerline strains and ^-dependent terms can be separated: 


S = e 0 + zk 


Here, e 0 represents the vector of midplane strains and k represents the vector of curvatures. 
The e 0 -K vector can then be written as a differential matrix operator acting on the midplane 
displacements: 


K 


r JL 

dx 

o 

d_ 

dy 

0 

0 

. 0 


0 

d_ 

dy 

d_ 

dx 

o 

0 

0 


0 - 

0 

0 

9 2 
9a: 2 

d 2 
dy 2 

_2 92 

dxdy 


u 0 

V 0 

W„ 


L o4' v (x,y)p(t.) 


where Rayleigh-Ritz assumed solutions are inserted for the centerline displacements and the 
differential matrix operator is defined here as L 0 . In [Crawley and Lazarus 91] this matrix 
is called D but is mistyped, as transposed. [Hagood, Chung and Von Flutow 90] create a 
similar differential matrix operator, L„, the difference being that L„ incorporates the factor 
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z, 3 



Figure B-6: Cutting the plate and integrating over an area gives the stress resultants, N, M and 
Q for the deformed plate. 

of s multiplying the curvature terms to produce the S vector rather than the e 0 -K vector. 
For a laminated plate of many layers, it is convenient to separate out the ^-dependence 
explicitly in calculating the stress resultants. 

Before proceeding to Hamilton’s Principle, we note one more construction for the lami¬ 
nated plate problem, the stress resultants, as shown in Figure B-6. By taking a section of a 
plate and integrating the stresses through the thickness, the stress resultants can be found. 
N x , N y and N xy are the in-plane stress resultants, M x , M y and M xy are the moment stress 
resultants and Q x and Q y are the shear stress resultants (assumed small): 



Plugging in the constitutive relations and taking the coordinate system for the laminate as 
shown in Figure B-6, the stress resultants become: 
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of [Hagood, Chung and Von Flutow 90] to those of [Crawley and Lazarus 91]. Hamilton’s 
Principle can be written: 

*2 

j s (T km -U+ WpJec + W) dt = 0 

tl 

For a system driven only by electrical excitation with no external forcing, W = 0. The 
variation in electrical work done on the system is also zero: 


sw etec = J b Qe t d) dV = J (D • 6E) dV = 0 

V V 


as there is no variation in electric held if it is prescribed and where we have used the relation: 


c f df 

bf = —bx 
ox 


The variation in internal stored energy is given as: 


bU 


j b Qs t t) dV = j (T • tfS) dV 

v v 


which is equal to: 


bU = 


N+ M# 


x y 


be 0 

bn 


Sect 


bn t 


■ ]\je/ec ■ 

M elec 


dxdy 


where the t subscripts signify a transpose of the vector. In terms of the resultant stiffness 
matrices, A, B and D, the variation in internal stored energy is: 


bU 


x y 


A B' 


rr\ 

o 

_1 

B D. 


. \ 


^o t 


bn t 


■ j\je/ec ■ 

M elec 


dx dy 


Replacing the midline strains and curvatures by the differential matrix operator and the 
Rayleigh-Ritz assumed solutions gives: 


bU = j J |[L 0 ^] T p(t) 


x y 


A B 
B D 


[Lotpp] <5p - [L otppY 


j\je/ec 

M elec 


^p (t) > dxdy 
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The first term on the right can be written: 


Kp (t) Sp (/) 


where K is as used in [Crawley and Lazarus 91]: 


K = 


x y 



B' 

D 


[LoV’p] dxdy 


The second term on the right has a dependence on the electric field. As voltage is the 
negative of the gradient of the electric field, E can also be expressed as a differential operator 
acting on a matrix of assumed potential shapes and a generalized electrical coordinate v(/): 


E = L v dp v (x, y, z) v (/) 


The components of the stress resultants due to electromechanical coupling become: 


j\je/ec 

M elec 


= £ 

k =1 


e* (z k ~ z k - 1 ) 
e t ( - Zk+ * k ~ 1 ) (z k - z k -i) 


L v tp v v(t) 


The second term in the expression for the variation in internal stored energy then is: 

et (z k - z k - i) 


x y 


[Lot/’p] ^2 | (zy+zy^) ( y \ 

k=i Let k 11 (z k - z k - r)J 


L v if) v dxdyv (t) ^p (t) 


This statement can be rewritten as: 


#v (t) ^p (t) 


where 0 is the electromechanical coupling matrix formulated in [Hagood, Chung and Von 
Flutow 90]: 


0 = 


[L 0 vg T £ 


x y 


k = 1 


et (z k - z k - X ) 
- ( -^±f^(z k -z k - 1 ) 


L v %l) v dxdy 


If Tkin f° r a plate is zero and the plate is only deflected statically by an applied v, 
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Hamilton’s Principle reduces to: 


*2 


(Kp — Ow) Spdt = 0 


Since the variation, <5p, is arbitrary, it must be the case that: 


Kp - 0v = 0 


which is the equation of motion for the statically driven plate. The modal amplitudes can 
be found by inverting the stiffness matrix: 

p = K -1 #v 


The final deformations are then calculated by multiplying the mode shapes by the modal 
amplitudes and summing. 

For a dynamic plate: 

Tkin = J ^pu T udV 


V 


Ignoring rotary inertia and assuming u is the vector of centerline displacements, u 0 , and 
substituting Rayleigh-Ritz assumed solutions, gives the variation in kinetic energy as: 


ST km = J ptppPtp p SpdV 


v 


Integrating by parts: 


h 


h 


J STkindt = - J J pipp TpppSpdVdt 

h h v 

The mass matrix as defined in [Hagood, Chung and Von Flutow 90] is then seen to be for 
a plate: 

M = J JtppJ^Pk (zk ~ z k -\)i> p dxdy 


x y 


k =1 


Hamilton’s Principle now becomes: 


*2 


(—Mp — Kp + 9v) Spdt = 0 
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Because the variation in the modal amplitudes, <5p, is arbitrary, the dynamic equations 


of an electromechanically coupled plate result: 


Mp + Kp - dv = 0 


By calculating the M, K and 0 matrices for a laminated plate as derived, the time evolution 
of the dynamic system can be computed. 

B.2.1 Simulation Results 

The model of [Hagood and McFarland 95] was implemented for an asymmetric free-free 
composite piezoelectric plate in cylindrical coordinates and compared to measured results 
of a free-standing stator, Greyl. 

In cylindrical coordinates, plane-stress assumptions lead to displacements of the form: 

u(r,0,t) = u 0 - z-q- 

( a +\ 1 dw ° 

v(r,0,t) = v 0 - z —— 

r oti 

w(r,9,t) = w 0 

where u, v and w are displacements in the r, 0 and z directions, respectively. Strain- 
displacement relations in cylindrical coordinates, [Timoshenko and Goodier 87], are derived 
from cutting a cylindrical section of a thin plate as in Figure B-8. Figure B-8(a) illustrates 
the radial and tangential deformations while Figure B-8(b) depicts in-plane shear defor¬ 
mations. The deformation of section abed to a'b'c'd' involves a rotation, which must be 
subtracted from the total angular change, yff + fy- The plane-stress strain-displacement 
relations become: 


'Sr 


r i 

dr 

s 2 

= 

1 dv _i_ u 

r dO ' r 

.Se. 


1 du \ dv v 

~ r dO ' dr r " 
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dr 


Figure B-8: In cylindrical coordinates, the strain-displacement relations take on a new form, (a) 
The in-plane extensional strain is a function of radial displacement, u, while the in-plane tangential 
strain is a function of both the radial displacement, u, and the tangential displacement, v. (b) 
In-plane shear strains formulations must make sure to subtract off rotations. 
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which results in an e 0 -n vector of: 


£<o 

K 


1 1 d_ 

r r dO 

o o 

o o 

o o 

. o o 


o 

o 

o 


_j9i 

dr 2 

1 _d_ _ 1 d 2 N 

r dr r 2 dO 2 , 

2 d 2 , 2__d_ 

r dOdr ' r 2 dO 


V 0 

W n 


If we are interested in the out-of-plane deflection of the neutral axis of the plate, ignoring 
small in-plane extensions, we solve only for w 0 . We take the Rayleigh-Ritz assumed solutions 
for w 0 for the annular plate as a sum of products of mode shapes (independent of time) and 
modal amplitudes (independent of geometry): 


w 0 (r,e,t)='^2'ip Pi (r,e)pi(t) 

i 

To simulate our motor, we choose two mode shapes to meet the geometric free-free boundary 
conditions at the inner and outer circumferences, one with no radial variation and one 
linearly proportional to r: 


d P = [cos{U) (fry) cos ( 30)] 

For an unsymmetric (in the ^-direction) plate, the e 0 vector is not the vector of midplane 
strains, but more precisely, the neutral axis strains. The neutral axis location, g, from the 
bottom of the plate is taken as a ratio of weighted stiffnesses: 

/ zcy^dz 

f4?dz 

where Z(k-i) °f the 0 th layer is then equal to —g. The L 0 matrix can now be found and the 
M and K matrices calculated. 

The electromechanical coupling term, ffv, determines the forcing in this system. The 
piezoelectric layer is segmented and poled alternately in such a manner as to induce a three- 
wavelength bending mode, as shown in Figure 7-6. The generalized electrical coordinates, 
v(t), are taken as the applied voltages, V, while the mode shapes are either +1 or -1 
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depending on the poling direction of the segment. As the electric held is the negative of the 
gradient of the potential, the differential operator, L„, can be taken as the reciprocal of the 
layer thickness. The 0 matrix in then computed using: 


E 3 = ± 


V 


Zk Zk—l 


For a static plate deformation, the modal amplitudes are found from: 


p = K _1 0v 


and combined with the mode shapes to find the neutral axis deflection in the ^-direction, 

W 0 - 

w 0 {r , 0, t) = tj} pi (r, 9)pi (/) + ip P2 (r, 0)p 2 (/) 

Dynamical properties of the system can also be calculated. System characteristics are 
found by examining the undriven response of the system. For no forcing, the equation of 
motion is simply: 

Mp + Kp = 0 


When we take: 

Pi (t) = P t e 3UJt 

we have an eigenvalue problem and can solve for the eigenfrequencies and mode shapes by 
taking derivatives: 

[-w 2 M + K] P = 0 

Natural frequencies are computed by calculating the eigenvalues of the K -a; 2 M matrix: 

[K -uj 2 M] = 0 

The solution produces repeated roots which correspond to orthogonal modes of equal 
amplitude variation, R(r). A new mode shape matrix is created to model the appropriate 
sin and cos modes for generating the traveling wave: 

ip p (r, 6) = [ R (r) cos (3 6) R (r) sin (3 6) ] 
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The radial dependence of each mode, R(r ), is composed by using the eigenvector, P, to 
weight the two radial dependence terms of the initial assumed-modes matrix, ip p : 

* (i , = p 1 + p 2 (^) 

The M, K and 9 matrices are then recalculated using the new if) mode-shape matrix. 
The new M and K matrices then become diagonal and we assume a proportional damping 
matrix C: 

Ca = 2 

where £ is the damping coefficient related to the quality factor by Q = The complete 
dynamical system can then be expressed as: 

Mp + Cp + Kp = 0v 

This equation can be re-expressed in terms of the modal amplitudes as: 

p = M -1 Cp - M -1 Kp + M _1 0v 

where this second order differential equation can be reduced to two first order differential 
equations by taking: 

n = p 

The equation of motion can then be written in state-space form as: 


n" 


i 

s 

1 

o 

i 

s 

1 


" n" 


' M -1 0v' 


= 




+ 


P. 


I o 


-P- 


0 


Once the modal amplitudes, p, are found, they are combined with the modes of ij)' p to 
calculate w 0 \ 

w 0 (r, 6, t) = R(r)cos (30) pi (/) + R(r)sin (39)p 2 (/) 

Figure B-9 shows simulation results for the driven response of a neutral axis out-of- 
plane displacement, w 0 , of a composite plate similar to stator Greyl, consisting of PZT-5A 
material, product number PSI-5A from Piezo Systems, Cambridge, MA, [Piezo Systems 95], 
bonded onto a steel substrate. The composite plate has an inner diameter of 5 mm and an 
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Figure B-9: A cylindrical plate with a piezoelectric layer bonded onto a steel substrate. The 
piezoelectric segments are driven dynamically at 60 Vpeak, 90 degrees out of phase at 18 kHz. The 
damping is 0.2%. 


outer diameter of 8 mm. The PZT layer is 0.191mm thick and the base height of the stator 
is 0.381 mm thick. The density of PSI-5A used in the simulation is 7750 and the density 
of stainless steel is 7860^-. The measured value of stator Greyl’s damping coefficient, 
£ = 2% was used for the calculation of the modal damping matrix, C. The initial conditions 
are taken at rest and the plot in Figure B-9 is drawn at 0.01 seconds after turn-on. 

The materials matrices characterizing the PZT layer are the plane-stress piezoelectric 
matrix, e* pzT - 
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■ 0 
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and the plane-stress stiffness matrix, Cp* ZT 
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For stainless steel, the Young’s modulus is 190 X 10 9 giving the isotropic material a 
stiffness matrix of: 
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The simulation predicted a resonant frequency of 18.0kHz and a deflection of approx¬ 
imately 20 n m at resonance. The actual measured natural frequency at the third bending 
mode was 42.2kHz and the amplitude of deflection was on the order of 4.5 /im. The dis¬ 
crepancy in natural frequency can be attributed to the fact that the materials parameters 
used in the simulation may not be accurate, the teeth of the stator are not modeled, the 
interconnects on the electrode patterns are not modeled, the electrodes are not perfectly 
matched in terms of resonance characteristics and/or the epoxy bond layer is ignored. As 
the deflection amplitudes are dependent on the drive frequency, we expect the predicted 
amplitudes to be off if the predicted natural frequencies are not correct. We also know that 
stator Greyl does not exhibit complete uniformity in traveling wave amplitude at all points 
around the stator, so it is not acting in quite the manner that it is modeled. 

Rotor coupling can be appended to this model by calculating modal forces in the direc¬ 
tions tangential and normal to the stator surface. [Hagood and McFarland 95] illustrates 
this process for a symmetric ring-type motor clamped at the inner circumference where the 
equations of motion for the stator are written as: 


-M-P T Cp T Kp — 0V T F modal,N T ~F m odal,T 

and the rotor is modeled by equations of motion in the rotary direction, parameterized by 
the angle, a: 

Irotor^ T F a QL — 7~int F B T 

and in the ^-direction, parameterized by the flexure height, Wf , the distance between the 
undeformed stator height and the rotor lower surface: 


M ro f or Wf T C z Wf — Fint F/v 


where Ti nt and Fi nt are the interface torques and forces coupling the stator and rotor models. 
The rotor equations of motion can be written in state-space form as: 
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and combined with the stator model to model the entire assembly. 

This method of [Hagood and McFarland 95] can be extended here for an asymmetric 
free-free ring and the contact models discussed earlier can be superimposed to compare their 
predictions and yet take into account the decrease in stator deflection due to the applied 
normal force and the extraction of mechanical work. 
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